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Yield and suppression of electrons from open heavy-flavor decays in heavy-ion 

collisions 
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Measurements by the STAR and PHENIX collaborations indicate that a quark- 
gluon plasma, a hot and dense state of matter in which quarks and gluons are not 
confined inside hadrons, is formed in heavy-ion collisions at the Relativistic Heavy Ion 
Collider. Charm and bottom quarks have been predicted to interact with the medium 
differently than the light quarks; a study of heavy quark interactions with the medium 
provides an important test of theoretical models of the quark-gluon plasma. 

The spectrum of non-photonic electrons and positrons (e^ is dominated by e ± 
from the semileptonic decays of D and B mesons. Therefore, non-photonic e ± serve as 
proxies for heavy quarks. A measurement of the modification of the non-photonic e ± 
spectrum in nucleus-nucleus collisions relative to p+p collisions allows the interactions 
of heavy quarks with the medium to be studied. Previous measurements indicate that 
high-transverse-momentum non-photonic e ± are suppressed in Au + Au collisions at 
\/~s~nn = 200 GeV relative to p + p collisions at the same energy. The magnitude 
of that suppression is larger than was anticipated and it has been a challenge for 
theoretical models to predict the in-medium energy loss of light and heavy quarks 
simultaneously. 

This dissertation presents the first measurement of the yield of non-photonic e ± 
from open heavy-flavor decays in Cu + Cu collisions at v /s/v r /v = 200 GeV and the 
suppression of that yield relative to p + p collisions. A comparison of this result to 
similar results for Au + Au collisions provides some indication that the geometry of a 
heavy-ion collision affects the average amount of energy loss by heavy quarks passing 
through the quark-gluon plasma. 
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Chapter 1 
Introduction 



1.1 The Standard Model 

Throughout the early and mid-20th century, as the energies of particle accelerators 
increased and cosmic-ray measurements became more refined, a large number of par- 
ticles were discovered. Questions arose as to whether all observed particles were 
truly fundamental. In 1964, Zweig and Gell-Mann independently suggested that all 
hadrons could be understood as being combinations of fundamental particles, which 
Gell-Mann called "quarks." [3) It was hypothesized that protons, neutrons, and the 
other baryons are made up of three quarks, while mesons are made up of a quark and 
an antiquark. Initially, the existence of three flavors of quarks - down (d), up (u), 
and strange (s) - was proposed as an explanation for the observed properties of the 
baryons and mesons then known. A serious theoretical objection raised against the 
early quark model was that the Pauli exclusion principle should preclude the exis- 
tence of particles like the A ++ and A~. The quark model proposed that those spin-| 
particles were composed of three identical fermions (uuu for the A ++ and ddd for 
the A - ) with their spins parallel, i.e., three identical fermions in the same quantum 
state, a situation not allowed by the Pauli exclusion principle. In 1964, Greenberg 
proposed the existence of a new quantum number, color charge, to resolve this issue. 
The existence of three color charges, called "red," "green," and "blue," was proposed; 
each of the three u quarks in the A ++ would have a different color, and the Pauli ex- 
clusion principle would not be violated. The first direct experimental support for the 
quark model came in the late 1960s: deep inelastic scattering experiments, studies of 
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electron-proton scattering with large momentum transfer, indicated that the proton 
is not a point particle, but contains three point-like scattering centers. [3}|4] 

The quark model has been expanded to include six flavors of quarks and has 
become an integral part of the Standard Model of particle physics, in which the inter- 
actions among the various flavors of quarks and leptons are mediated by the exchange 
of vector bosons. The photon mediates the electromagnetic interaction between elec- 
trically charged particles, the W and Z° bosons mediate the weak interaction, and 
gluons mediate the strong interaction between quarks. Gluons themselves carry color 
charge: specifically one unit of color and one unit of anticolor. The Standard Model 
includes three of the four fundamental interactions; only gravity is excluded. The 
only particle in the Standard Model that has not been discovered is the Higgs boson, 
which has been the subject of an intense search effort at the Tevatron and now at the 
LHC. 

1.2 Quantum Chromodynamics 

The interactions of quarks and gluons is described within the framework of quan- 
tum field theory by Quantum Chromodynamics (QCD). QCD is a non-abelian gauge 
theory with each flavor of quark assigned to the fundamental representation of the 
577(3) symmetry group. The QCD Lagrangian |5) is the Yang-Mills Lagrangian for 
the SU (3) group: 

£ = -iV™ + $(i^D„ - m)V>. (1.1) 
Here, the covariant derivative is 

D, = d,-igAlt\ (1.2) 

where g is the coupling constant, A a ^ represents the gauge field (gluons), and the t a 
are the generators for SU(3). In Equation 1.1 F^ v is the field strength. 



[D„ D v \ = -igFp* =► F; v = d,A a u - d v A« + gf abc 4Al, (1.3) 
where the f abc are the 577(3) structure constants: [£ a ,£ 6 ] = if abc t c . When the La- 



grangian in Equation l.l|is expanded, the last term in Equation 1.3 gives rise to three- 



and four-gluon couplings. The coupling of the gauge field to itself, which does not 
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occur in Quantum Electrodynamics (QED), results in the strong coupling constant 
being asymptotically free. 
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Figure 1.1: Measured values of the strong coupling constant a s as a function of 
the momentum transfer Q. [6j Open symbols indicate Next-to-Leading-Order (NLO) 
calculations, while closed symbols indicate Next-to-Next-to-Leading-Order (NNLO) 
calculations. 



The QED vacuum contains virtual particle-antiparticle pairs which arise from 
the vacuum, live a short time (consistent with the Heisenberg uncertainty principle), 
then annihilate. If the virtual particle and antiparticle carry electric charge, the pair 
constitutes an electric dipole. In the presence of a real electric charge, the virtual 
dipoles are reoriented to be antiparallel to the field of the real charge, resulting in 
the partial screening of that charge. The vacuum functions as a dielectric and the 
effective electric charge seen at large distances (or small momentum transfers) is less 
than the "bare" charge seen at small distances (or large momentum transfers). [3] 
A similar situation arises in QCD: virtual qq pairs in the vacuum will couple to the 
gluon field produced by a real quark and partially screen that quark's color charge. 
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However, QCD also includes couplings of gluons to gluons. Those virtual gluon loops 
are found to give rise to an "antiscreening" effect, which is larger than the screening 
effect. The QCD coupling constant therefore decreases with decreasing distance and 
increasing momentum transfer q. The coupling constant a s = g 2 /^ as a function of 
momentum transfer has been calculated [5l|6l to lowest order to be 



2\ 2tt 
aM ] = { l in c -ln f )Hq/A) = (11 - \n f ) lnfa/A) ' (L4) 
where n c = 3 is the number of colors, n/ is the number of flavors, and A 200 MeV is 
the QCD scale. Figure 1.1 showg^j measurements of the strong coupling constant (6j; 



a s has been observed to decrease with increasing momentum transfer in a manner 
consistent with theoretical predictions. When q is significantly larger than A, QCD 
can be treated perturbatively; this is the case for the production of c and b quarks in 
heavy-ion collisions. 

When q < A, perturbation theory breaks down and other approaches must be 
used to perform calculations in QCD. The most productive approach to performing 
non-perturbative QCD calculations has been Lattice QCD. Continuous spacetime is 
replaced with a set of discrete points on a four-dimensional Euclidean lattice. Calcu- 
lations are performed numerically by evaluating path integrals between neighboring 
points on this lattice. [5] Lattice QCD calculations indicate that for large distances 
(or small momentum transfers) the potential between a quark and an antiquark in- 
creases linearly with the distance separating them (cf. the potential in the weak- 
coupling limit, which is a coulomb potential with a running coupling constant). As 
the separation distance between a quark and antiquark increases, the potential energy 
increases without limit. At some point it becomes energetically favorable for a second 
quark- antiquark pair to be created. If a quark-antiquark pair is separated, the result 
is not an isolated (anti) quark, but two mesons. Observing an isolated quark seems 
to be impossible and all particles that have been directly observed are color singlets. 

While quarks are confined in the states of matter normally observed, asymptotic 
freedom implies that in nuclear matter at high temperature, hadrons should "melt" 
and the new degrees of freedom of the system should be deconfined quarks and gluons, 
rather than hadrons. Such a state of matter is called a "Quark-Gluon Plasma" (QGP). 
The transition from a hadronic state to a QGP is accompanied by an increase in the 
number of degrees of freedom, implying an increase in the entropy density and pressure 



1 Color versions of all Figures are available in the electronic version of this dissertation. 
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as the temperature increases through the transition region [7j. For an ideal gas of 
quarks and gluons (in which the quarks are massless and the deconfined quarks and 
gluons do not interact with each other), the equation of state (81 is 
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where P is the pressure and T is the temperature. Lattice QCD calculations [8] show 
a rapid increase in P/T 4 with increasing T in the vicinity of a critical temperature 
T c ss 160 MeV. For temperatures > 2T C , P/T 4 saturates at values somewhat less than 



the ideal (Stefan-Boltzmann) limit given in Equation 1.5 Because the QGP does not 
reach the limit of an ideal partonic gas at these temperatures, it is sometimes referred 



to as a "strongly interacting QGP" (sQGP). These findings are shown in Figure 1.2 
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Figure 1.2: The ratio of pressure to T 4 calculated [8] using lattice QCD as a function 
of temperature. The different curves indicate the number of active quark flavors; "2+1 
flavor" indicates that the s quark is heavier than the d and u, while "pure gauge" 
indicates active flavors. The arrows indicate the Stefan-Boltzmann predictions of 
P/T 4 for an ideal partonic gas for each set of active quark flavors. 



Figure 1.3 is a schematic phase diagram for strongly interacting matter [9], with 
the temperature on the vertical axis and the baryo-chemical potential (ps) on the hor- 
izontal axis. Many details, including the numerical values and orders of the proposed 
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Figure 1.3: A schematic phase diagram for (thermally equilibrated) strongly inter- 
acting matter [9] with temperature on the vertical axis and baryo-chemical potential 
potential on the horizontal axis. The matter produced in heavy-ion collisions is ex- 
pected to lie at low baryo-chemical potentials and high temperature (in the vicinity 
of the "crossover" label). 



phase transitions, are the subjects of debate and active research in the theoretical and 
experimental communities. Nevertheless, the qualitative form of the phase diagram 



shown in Figure L3 is generally accepted. For low temperatures and baryo-chemical 
potentials, strongly interacting matter is in a hadronic state, either as a hadron gas or 
as nucleons in atomic nuclei. For low temperatures and \ib larger than nuclear values 
(conditions expected to be present in neutron stars) a variety of ordered quark matter 
phases is predicted. As described above, at high temperatures strongly interacting 
matter is expected to exist as a quark-gluon plasma. The nature of the transition 
from the QGP to the hadronic phase is not fully understood. Lattice QCD calcula- 
tions at fiB = tend to favor a crossover rather than a sharp phase transition, while 
other models for larger //g = favor a first-order phase transition. These results 
suggest the existence of a critical point in the phase diagram. This critical point may 
be in the region of the phase diagram accessible to heavy-ion collision experiments. 
An effort is underway at the Relativistic Heavy Ion Collider (RHIC) to find the crit- 



ical point. 10 This dissertation is concerned with the strongly interacting matter in 
heavy-ion collisions, which exists at high temperature and low baryo-chemical poten- 
tial. Such conditions are also predicted to have existed in the early universe, < 10/is 



after the Big Bang 11 and before the formation of hadrons. 
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At RHIC (see Section 2.1) gold nuclei (and copper nuclei) are collided at center- 
of-mass energies of 200 GeV per nucleon-nucleon pair (denoted as a /s/vtv=200 GeV). 
These collisions introduce a large amount of energy into a region of space the size of 
an atomic nucleus (and much larger than the volume of a single hadron). At RHIC, 



the highest energy heavy-ion collisions are estimated 12 to reach energy densities of 
~ 4.9 GeV/fm 3 . This is seven times the critical energy density of ~ 700 MeV/fm 3 
(predicted by Lattice QCD calculations (8)) necessary for the formation of a QGP. 



Figure 1.4 (upper) is a schematic spacetime diagram of the expected stages of a heavy- 
ion collision for two scenarios, with the right-hand side including the formation of a 



QGP. [13] Figure 1.4 (lower) is a visual representation 14 of the various stages of a 



heavy-ion collision. The QGP formed in a heavy-ion collision is generally assumed to 
reach thermal equilibrium (at least locally), although there is debate on this point. 
The QGP expands and cools, eventually reaching the critical temperature, where 
the system makes the transition to a hadron gas. This transition may be a sharp 
phase transition, or a crossover during which QGP and hadronic gas coexist (called 
the "Mixed Phase"). As described above, recent Lattice QCD calculations tend to 
favor a crossover for low baryo-chemical potentials (values of fis expected for RHIC 
collisions). After the transition from a QGP to a hadron gas, the hadrons continue to 
interact with each other inelastically until the system reaches the "chemical freezeout" 



temperature (T c h in Figure 1.4). When the system reaches chemical freezeout, the 
ratios of the yields of the various particle species are fixed. The hadrons continue 
interacting elastically until the system reaches the "thermal freezeout" temperature 



(Tf in Figure 1.4), after which there are no hadronic interactions. The shower of 
particles produced (for Au + Au collisions at ^3^^=200 GeV, approximately 1000 
particles per unit rapidity at mid-rapidity) continues outward from the collision vertex 
and may be detected by experiment. 
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Figure 1.4: Upper plot: spacetime diagram 
ion collision. The right-hand side includes the formation of a quark-gluon plasma. 
Lower plot: Illustration (not to scale) of some of the stages of a heavy-ion collision at 



various times. 14 The incoming nuclei are heavily Lorentz contracted (the Lorentz 
factor 7 « 100 for the highest RHIC energies). The acronym "sQGP" indicates 
"strongly interacting QGP." 
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1.3 Heavy-Ion Collisions: Basic Concepts 



Some of the experimental signatures of the formation of a quark-gluon plasma will 



be discussed in Section 1.4 However, before that discussion the standard coordinate 



system and kinematic variables will be defined (Section 1.3.1 ) and the geometry of a 



heavy-ion collision will be described (Section 1.3.2). 



1.3.1 Coordinate System and Kinematic Variables 

The beam line is defined to be the z-axis (along which the incoming nuclei travel), 
with the x- and y-ax.es of a Cartesian coordinate system oriented horizontally and 
vertically, respectively. In a spherical coordinate system, the polar angle is referred 
to as 8, while the azimuthal angle is referred to as (p. The transverse momentum (j>t) 
of a particle is the component of its momentum in the direction perpendicular to the 
z- axis: 



Pt = ypl + p\ = psm.9 = p z tan#, (1.6) 

where p is the magnitude of the particle's momentum, and p x , p y , and p z are the x-, 
y-, and z-components of its momentum, respectively. The transverse momentum is 
invariant under Lorentz boosts along the beam direction. The rapidity of a particle 
is generally defined to be 

rapidity = arctanh- = arctanh— = Mn — - , (1.7) 
F y c E 2 E-pc v ' 

where v is the particle's velocity, E is its energy, and c is the speed of light in 

vacuum. However, this dissertation uses a modified definition of the rapidity: the 

rapidity relative to the beam axis. This quantity, also denoted as y, is 

y = Mn — . 1.8 

y 2 E-p z c 1 ; 

For particles moving at ultrarelativistic velocities, the rapidity y is well approximated 
by the pseudorapidity r\: 



y 7] = | In - — -^ = — In tan-. (1-9) 



P + Pz , , 9 
= — In tan-. 

P-Pz 2 

Unlike the rapidity, the pseudorapidity is independent of the particle's mass. This 
thesis is primarily concerned with e ± with momenta > 2GeV/c; the pseudorapidity 
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will usually be used in place of the rapidity, and instead of the polar angle 6. For a 
particle with mass m, the "transverse mass" is defined to be 



rriT = W p^c 2 + m 2 . (1-10) 

For a heavy-ion collision, the energy quoted is typically the center-of-mass col- 
lision energy per nucleon-nucleon pair y/sNN- Collision energies referred to in this 
dissertation are values of ^snn unless otherwise noted, even if the symbol "^snn" 
has been omitted for the sake of brevity^] Furthermore, this dissertation frequently 
mentions RHIC data, many of which are for collisions at v /s^v r /v=200 GeV; if no en- 
ergy is mentioned in reference to a collision, it should be assumed that the energy is 
200 GeV. 



1.3.2 Collision Geometry 

When two particles traveling along the z-axis collide, the plane including the z-axis 
and the centers of the two particles is called the "reaction plane." The distance in the 
transverse plane between the centers of the colliding particles is called the "impact 
parameter." For heavy-ion collisions, many quantities depend upon the impact pa- 
rameter. Head-on collisions, those with small impact parameter and a large overlap 
between the nuclei in the transverse plane, are said to be "central" collisions, while 
collisions with large impact parameter and a small overlap are said to be "peripheral." 
In central collisions, the number of nucleons that participate in collisions (called N part , 
the "number of participants") is large, while in peripheral collisions, that number is 
small. Therefore, for a given pair of colliding nuclei, the amount of QGP produced 
is directly related to the centrality of the collision; many of the signatures of a QGP 
should be stronger in central collisions than peripheral collisions. Centrality is typi- 
cally quoted as a percent of the geometric cross-section, e.g. the "0-20% most central 
collisions." Related to N part is the number of binary nucleon-nucleon collisions, which 
is called N^n in this dissertation, but is also called N co u in some references. 



N par t and Nu n are estimated from simulation using a Glauber model 15 of 
nucleus-nucleus interactions. Nuclei are described 16 as a collection of nucleons 
distributed in space according to a density profile, typically (for spherical nuclei) 



2 For example, "200-GeV Cu + Cu collision" means "Cu + Cu collision at ^5^=200 GeV.' 
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P{r) = Po 



1 + exp 



i.ir 



a 

where R is the nuclear radius and a is the "skin depth." It is assumed that in high- 
energy collisions, the nucleons in the colliding nuclei carry enough momentum that 
they are undeflected when the nuclei pass through each other. The nucleons there- 
fore travel in independent linear trajectories along the z-axis. The nucleon-nucleon 
scattering cross-section is assumed to be the same as the inelastic p + p scattering 



cross-section, <7j ne /. = 42 ± 2 mb at y/s = 200 GeV. 17, 18 Given the nuclear density 
profile and the nucleon-nucleon cross-section, the number of collisions experienced by 
a single nucleon can be calculated. The STAR collaboration performs Glauber-model 
calculations using Monte Carlo simulations. Nucleus-nucleus impact parameters are 
randomly generated, and nucleons are randomly distributed in space according to 
the nuclear density profile. A nucleon-nucleon collision occurs if two nucleons come 
within a distance of \/<Ji ne i,/ir of each other. The mean values of N part and Nu n 
((Npart) and (Nu n )) for a given range of impact parameters (centrality class) can be 
determined after simulations of many collisions. The "mean nuclear thickness func- 
tion" is (Taa) — (Nbin) I Vinei- ft is a l so possible to calculate the area and eccentricity 
of the nuclear overlap region. 

In experiments, the impact parameters, and therefore centralities, of heavy-ion 
collisions cannot be controlled. The multiplicity of charged particles produced in a 
collision increases with decreasing nuclear impact parameter (i.e., as the collisions 
become more central). Collisions may be grouped into centrality classes based on 
the charged-particle multiplicity. The STAR collaboration defines centrality classes 
based on the "reference multiplicity," the number of charged tracks (not corrected 
for tracking efficiency) measured in the pseudorapidity range \r]\ < 0.5. Reference 
multiplicity values are ~ 700 for central Au + Au collisions and ~ 400 for central Cu 
+ Cu collisions at 200 GeV. 



1.4 Experimental Signatures of a QGP 

This section will describe a few of the experimental results that indicate that a quark- 
gluon plasma is formed in heavy-ion collisions at RHIC. 
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1.4.1 Jet Quenching 
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Figure 1.5: High-p^ dihadron azimuthal correlations: the difference in azimuth be- 
tween a charged, high-p^ trigger particle and other charged particles in the same 
event. In the left-hand plot are the azimuthal correlation data measured in 200-GeV 



p + p (black), central d+Au (red), and central Au + Au (blue) collisions. [19-21 The 



right-hand plot shows azimuthal correlation data for 200-GeV p + p (black) collisions 
and Au + Au collisions (colored). [22] The Au + Au data are divided into two classes: 
in the "in-plane" class (blue) the azimuth of the trigger particle is within 45° of the 
reaction plane; in the "out-of-plane" class (red) the azimuth of the trigger particle 
is > 45° from the reaction plane. The open symbols are reflections of the measured 
data (closed symbols) about A<p = and A0 = ir. 



Figure 1.5 shows measurements 20 21 by the STAR collaboration which, in the 
first few years of RHIC operation, were among the strongest indications that a QGP 
is formed in heavy-ion collisions. The distributions are A0, the difference in azimuth 
between a high-pr "trigger" particle (with px > 4 GeV/ c) and "associated" particles 
from the same event which have pt > 2 GeV/ c. When the trigger and associated 
particles are in the same jet, a peak near A0 = is produced, called the "near-side" 
peak. When the trigger and associated particles come from back-to-back jets, a peak 
near A0 = 7r is produced, called the "away-side" peak. Figure |1.5| shows that for 
200-GeV p + p and d+Au collisions (in which no QGP would be produced), both the 
near- and away-side peaks are present. However, as shown in the left-hand panel, for 
central Au + Au collisions the away-side peak is not present. This is interpreted as 
the absorption (or "quenching") of one member of a back-to-back jet pair in a QGP. 
Consider a qq or gg pair that becomes back-to-back jets and is produced near the 
surface of the medium. One parton is directed out of the medium into the vacuum and 
fragments into the near-side jet, largely unmodified by the presence of the medium. 
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The other parton passes through a large amount of QGP, which absorbs its energy, 
leading to a reduction in associated particles with px > 2 GeV/ c and the observed 
quenching of the away-side jet. Since the away-side peak is present in d+Au collisions, 
the quenching cannot be due merely to the presence of (cold) nuclear matter in the 
collision region. 



Jet quenching is further illustrated in the right-hand plot of Figure 1.5, which 



shows dihadron azimuthal correlations for two subsets of Au + Au events in the 20- 



60% centrality class. 22 In the "in-plane" class, the azimuth of the trigger particle 
is within 45° of the reaction plane, while in the "out-of-plane" class the azimuth 
of the trigger particle is > 45° from the reaction plane. For a non-central heavy- 
ion collision, the overlap region between the nuclei is azimuthally anisotropic, as 



is the QGP produced (see Figure 1.6). Projected into the transverse plane, the 
overlap region is roughly elliptical or "almond" shaped, with the long axis oriented 
perpendicular to the reaction plane. A parton produced perpendicular to the reaction 
plane will therefore have a greater path length through the medium, leading to greater 
energy loss. The jet quenching effect should be larger when the trigger particle is 
roughly perpendicular to the reaction plane than when the trigger particle is parallel 



to the reaction plane; as observed in the right-hand plot in Figure 1.5 



1.4.2 Elliptic Flow 

Results from RHIC also indicate that the medium produced in a heavy-ion collision 
can be described by ultrarelativistic hydrodynamical models with low or viscos- 
ity (i.e., a "perfect fluid"). In a non-central nucleus-nucleus collision, azimuthal 
anisotropy in the shape of the overlap region results in an azimuthal anisotropy in 
the momentum distribution of the emitted particles. The solid curves in Figure [l~6l are 



hydrodynamical calculations 23 of curves of constant energy density within a QGP 
that can be described as a perfect fluid. The pressure gradients within the QGP 
will impart additional transverse momentum to the particles emitted from it (a phe- 
nomenon called "radial flow"). Since the pressure gradient is greater in the reaction 
plane than perpendicular to it, particles with azimuthal angle near the reaction plane 
receive a greater momentum boost due to flow. This results in a momentum-space 
azimuthal anisotropy in the spectrum of emitted particles, referred to as "elliptic 
flow." The azimuthal distribution of particles within a given transverse-momentum 



bin may be written in the form of a Fourier expansion 29 : 
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Figure 1.6: A hydrodynamical model [23] describing the overlap region of two Au 
nuclei (dashed circles) moving along the z-axis (into and out of the page) and colliding 
with an impact parameter of 7 fm. The overlap region is "almond" shaped, with the 
long axis oriented perpendicular to the reaction (xz) plane. The solid red curves show 
hydrodynamical calculations of curves of constant energy density within the QGP. 



d 3 N 



1 d 2 N 



= s i^fy { 1 + £ 2v - cos [»(*-*»■)] ) • < L 12 > 

Here, the v n are the Fourier expansion coefficient^] and ^rp is the azimuthal angle 
of the reaction plane in the lab frame. The coefficient for the second harmonic t>2 
is used as a measurement of elliptic flow. In heavy-ion collisions at RHIC, tends 



to increase as collisions become more peripheral 24 28 . This is consistent with the 
assumption that the momentum-space azimuthal anisotropy is due to the differing 
pressure gradients in and out of the reaction plane in a QGP: the eccentricity of 
the overlap region increases as collisions become more peripheral, resulting in greater 



difference in the pressure gradients. Figure 1.7 shows t>2 as a function of pr for 



identified and (anti)protons 24 2627 along with hydrodynamical calculations 25 



3 For a system that is symmetric with respect to the reaction pla ne th e sine terms in a Fourier 
expansion will be 0, hence the omission of such terms from Equation 
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Figure 1.7: Measurements of v% by the STAR collaboration [24] for charged pions 
and (anti)protons in 200-GeV Au + Au collisions. The curves are from calculations 
assuming a hadron-gas equation of state (dotted) and calculations assuming a perfect- 



liquid QGP and a phase transition to a hadron gas. 25 



of V2- For the solid curves, it is assumed that there is a phase transition between a 
QGP with zero viscosity and a hadron gas, while the dashed curves assume that the 
system is only described by a hadron-gas equation of state. The data are better 
described by the hydrodynamical calculations that assume the presence of a QGP. 



Figure 1.8 shows Vi for a variety of particle species 28 ; in the left-hand panel at low pr 
particles with lower mass tend to have a higher v^- In the middle panel, v-i is plotted as 
a function of ttit — m (the kinetic energy, neglecting longitudinal momentum) rather 
than pt'- all baryons appear to follow one trend, while all mesons follow another, and 
at low pt all of the particle species shown appear to follow the same trend. This 
baryon-meson separation can be eliminated by scaling V2 and rriT — m by a factor 
of l/n q , where n q is the number of valence quarks in the hadron (right-hand panel). 
This is taken to be an indication that if a fluid is created, it is the individual quarks 
that are flowing. 
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Figure 1.8: Measurements of v i by the STAR 
for various particle species. In the left-hand panel, v-i is plotted as a function of 
transverse momentum. In the middle panel, v-i is plotted as a function of — m 
(the "transverse kinetic energy"). In the right-hand panel, both the abscissa and the 
ordinate have been scaled by the number of valence quarks in each hadron. 
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GLV parton energy loss (dNrVdy = 1 100) 
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Figure 1.9: The nuclear modification factor (Raa) for neutral pions, 30 31 rj mesons, 
32 and direct photons 33 in central 200-GeV Au + Au collisions. Also shown is the 



suppression of light-flavor hadrons predicted by the GLV model 34 , 35 with gluon 
density dW/dy = 1100. 



As a parton passes through a quark-gluon plasma, it is expected to lose energy 
to the medium through gluon bremsstrahlung and collisions with other partons in 



the medium. 36 , 37 This is expected to result in a suppression of high-pr particles 



in heavy-ion collisions (relative to p + p collisions, in which no medium is present). 
The "nuclear modification factor" is used to quantify the effect of the presence of 
nuclear matter on particle yields. This factor is often denoted as Raa, although the 
subscripts can be used to indicate the specific collision system being considered (e.g., 
RdAu for d + Au collisions, .RauAu for Au + Au collisions, and -RcuCu for Cu + Cu 
collisions). The nuclear modification factor is the ratio of the yield of a particle species 
in nucleus-nucleus collisions to the yield of the same particle species in proton-proton 
collisions, with the p + p data scaled to account for the fact that there are many 
nucleon-nucleon collisions in an A + A collision. For a given transverse-momentum 
bin, the nuclear modification factor for particle species X is defined as 
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/ Ed 3 N x 



Raa = )„7 y / , (1-13) 



\ dp 



AA 



&inel. 



Ed 3 a x 



i'i> 



(N bin ) 



dp 3 

where Ed 3 N x A /dp 3 is the invariant yield of particle species X in nucleus-nucleus 
collisions, (Nun) is the mean number of binary nucleon-nucleon collisions for those 
nucleus-nucleus collisions, Ed 3 a x p j 'dp 3 is the invariant cross-section for particle species 
X in p + p collisions, and collisions being studied, and (Ti ne i. =42 mb is the inelastic 



p + p scattering cross-section 17 . The value of a/s/vtv is assumed to be the same 
for the A + A and p + p collisions being compared. If the presence of nuclear matter 
in the collision region has no effect on the spectrum of particle species X, an A + A 
collision can be viewed as a superposition of (AT 6in ) nucleon-nucleon collisions and 
the nuclear modification factor will be 1. Values of Raa greater than 1 indicate that 
particle species X is enhanced in A + A collisions relative to p + p collisions, while 
values of R X A < 1 indicate a suppression of particle species X. 



Figure 1.9 shows measurements by the PHENIX collaboration of the nuclear mod- 
ification factor for neutral pions, rj mesons, and direct photons (photons produced in 
the collision itself and not through the subsequent decays of hadrons) in central Au 



+ Au collisions. The direct photons 33 are not suppressed, which is expected since 
photons do not interact via the strong interaction and should not be affected by the 



presence of a quark-gluon plasma. The hadrons 30 32 are observed to be suppressed 
by about a factor of five at high p?. The curve is a theoretical calculation of light- 



flavor-hadron suppression in one model, the GLV model. 34 35 The quantity dN g /dy 
is the gluon density per unit rapidity, a parameter in the GLV model that is related 
to the opacity of the medium. The observed values of Raa for light-flavor hadrons 
were used to constrain the model; the GLV model with dN 9 /dy ~ 1100 best describes 
the data. 



1.4.4 Thermal Models 

Measurements of the ratios of particle yields in heavy-ion collisions indicate that the 
system has reached thermal equilibrium at the time it reaches chemical freezeout 
(which occurs after the system has passed through the phase transition from a QGP 
to a hadron gas). In statistical thermal models, the system is described by a grand 
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200 GeV 197 Au + 197 Au central collision 

Figure 1.10: The ratios of the yields of various particle species measured [7] by the 
STAR collaboration in central 200-GeV Au + Au collisions (red points). Also shown 
is a thermal-model fit to those data, which is used to find the temperature and baryo- 
chemical potential of the system at chemical freezeout. The inset shows how the 
strangeness-suppression factor 7 S changes with collision centrality. 



canonical ensemble, with the Hamiltonian usually taken to be that of a hadron reso- 
nance gas. From a simple expression for the grand canonical partition function, the 
mean number of particles of species X (which has mass mx, baryon number Bx, 
strangeness Sx, and electric charge Qx) is calculated [38] to be 



(Nx) 



2tt 2 



2 \J^r) 



cxp 



. BxfiB + Sxl^s + QxHq 



;i.i4) 



V and T are the volume and temperature of the system. The quantities /is, 
and /j,q are the chemical potentials associated with baryon number, strangeness, and 
electric charge. K 2 is the modified Bessel function. The factor gx is the spin-isospin 
degeneracy factor (39 1 for particle species X. The upper (lower) sign is to be used for 



fermions (bosons) and Equation 1.14 is written in units where c = h = ks = 1- The 
chemical potentials are not independent: the initial charges of the incoming nuclei 
and the requirement of zero net strangeness restrict the possible configurations of the 
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system, leaving V, T, and /!# as the only independent parameters. The volume of the 
system is unknown, but is removed from consideration if the ratios of particle yields 
are considered. The measurement of two particle ratios is sufficient to determine 
the two unknown parameters; in practice, many particle ratios are measured and the 
set of those measurements is fit to determine the most likely values of T and the 
baryo-chemical potential. The "strangeness-suppression factor" 7 S is also introduced 
to account for the fact that strange quarks may not be in full thermal equilibrium. 



Figure 1.10 shows the results of one such fit for central 200-GeV Au + Au collisions 
measured by the STAR collaboration. [7| The thermal fit with T = 163 ± 4MeX 
fiB = 24±4 MeVand 7 S = 0.99±0.07 describes the ratios well. The inset shows how the 
value of 7 S changes with collision centrality, with strangeness approaching equilibrium 
values as collisions become more central. The chemical freezeout temperature of 163 
MeV places a lower limit on the critical temperature T c at = 24 MeV. 
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1.5 Heavy Flavor 



Heavy quarks, defined for the purposes of this dissertation to be the c and b quarks, are 
useful probes of the medium produced in heavy-ion collisions at RHIC. Because their 
masses (m c = (1.27+°;^) GeV/c 2 and m b = (4.19+°:oe) GeV/c 2 ) jl7| are much larger 
than the RHIC temperature scale, heavy quarks are not expected to be produced 
thermally in significant amounts in the QGP. At RHIC, heavy quarks are expected to 
be produced primarily in the hard (large momentum transfer) scattering of partons in 
the initial stages of the collision. Heavy quarks are therefore expected to be present 
in the QGP and the subsequent hadron gas throughout their evolution. Furthermore, 
because of their large masses, heavy quarks have been predicted to interact differently 
with the medium than light quarks. A satisfactory model of the QGP should be 
able to describe gluon, light-quark, and heavy-quark interactions with the medium 



using the same set of parameters (see Section 1.5.2). This section describes some 
of the theoretical issues regarding heavy quarks and presents some of the major 
experimental heavy-flavor measurements from RHIC. Because of their short lifetimes, 
hadrons containing heavy quarks do not survive long enough to leave the RHIC 
beam line (much less reach the tracking components of the RHIC detectors). Heavy- 
flavor hadrons must therefore be studied indirectly through their decay products. 
The yields of D mesons have been measured by reconstructing their hadronic decays 
(see 40-42 1 ) . The yields of the J/^p and T have been measured by reconstructing 



their decays to e e + pairs (see |43}j48] ) . Some measurements of heavy-flavor hadron 



yields have also been performed by measuring their decays to muons (see 49,50]). 

Several studies (of which this dissertation is one) have measured the suppression 
(relative to p + p collisions) of heavy quarks at high p? by measuring the nuclear 
modification factor of non-photonic e in various collision systems. A non-photonic 
e (sometimes abbreviated as "NPE" and sometimes referred to as a "single e ) is 
defined to be an produced with an (anti) neutrino in a charged-current 
weak decay. In contrast, "photonic e are produced in pairs from sources including 
photon conversions, Dalitz decays (e.g., tt° — > e~e + 7), J/ip — > e~e + , the Drell-Yan 
process, and others. While non-photonic e may also come from sources other than 
heavy-flavor decays (charged pions, muons, kaons, etc.) the spectrum of non-photonic 
e is dominated by the products of heavy-flavor decays (see section 8.2). In this 
dissertation, the yield and nuclear modification factor of non-photonic e in Cu + 
Cu collisions at ^/snn = 200 GeV will be measured. 
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1.5.1 Production 
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Figure 1.11: The measured total charm cross-section (<r cE ) scaled by (N bin ) as a 
function of (N bm ) for 200-GeV p + p [51] , d + Au |40|, Cu + Cu [42], and Au + 



Au |50|[52] collisions in various centrality classes. Also shown are NLO [53] and 
FONLL |54| calculations of a c5 . 



Since heavy quarks are primarily produced through parton hard scattering, the 
production cross-section may be calculated using perturbative QCD. These per- 
turbative calculations are typically done at the Next-to-Leading-Order (NLO) or 
Fixed-Order-Next-to-Leading-Log (FONLL) levels. NLO calculations include con- 
tributions of order o? s and a^, while FONLL calculations include contributions of 
those orders along with logarithmic contributions of the form cu" log n_1 (pr/Mg) and 



a" log" (pt/Mq) resummed to all orders n (where Mq is the heavy-quark mass). |53p4 



Since the heavy-quark pt appears in the logarithmic terms, FONLL calculations must 
be performed in transverse-momentum bins and then integrated to obtain the total 
cross-section. The heavy quark is treated as an "active" flavor for px 3> Mq and 
included in the calculation of the coupling constant (see Equation 1.4). Calculations 
using this framework give af E 07VLL = 0.256^;^ mb and af 5 07VLL = 1.87^ fib for 
the total charm and bottom cross-sections in nucleon-nucleon collisions, 54 with the 
uncertainties in the cross-sections due to uncertainties in the heavy quark masses 
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and the choice of the renormalization scale. NLO calculations may be performed as 
described above for FONLL, but they may also be performed in a different frame- 
work in which the total cross-section is evaluated for all transverse momenta (i.e., 

division into pt bins is not necessary) and the heavy quark is treated as massive (and 

nlo _ n Qm +i-ooo 

cc — u - ou± -0.210 

for the total charm and bottom cross-sections in 



not counted in the evaluation of a s ). One such calculation gives cr^ LO = 0.301. ' nn " 



53 



mb and a^ LO = 2.06±^;|f /d) 
nucleon-nucleon collisions. 

The STAR and PHENIX collaborations have measured the total charm cross- 
section in 200-GeV p + p and A + A collisions. The STAR collaboration measures 
the yield of \(D° + D°) by reconstructing the D° — > n + K~(D — > n~K + ) decays and 
finding the invariant masses of pion-kaon pairs. A simultaneous fit of the D° mea- 
surements (mostly at ^ 2 GeV/ c) and non-photonic e measurements (at higher 
Pt) is performed to obtain an estimate of the total charm cross-section. The PHENIX 
collaboration estimates a C c from non- photonic e ± (which PHENIX can measure for 



Pt > 200 MeV/c). Figure 1.11 shows the experimental measurements of a c5 / (Nun) f° r 



p + p[51],<i + Au[40j, Cu + Cu [42], and Au + Au (50j[52j collisions at v /s^7=200 
GeV, along with the two calculated values of a c5 described above. For each experi- 
ment, the measurements indicate that the total charm cross-section scales with the 
number of binary collisions, consistent with the assumption that charm quarks are 
produced predominantly in the initial hard scattering of nucleons. This indicates that 
it is reasonable to use perturbative QCD (FONLL) calculations of the c- and 6-quark 
spectra as initial spectra for theoretical studies of heavy-quark interactions with the 



QGP (as is done in the studies described in Section 1.5.2). However, there appears 
to be a systematic factor of ~ 2 difference between the values of a c5 measured by the 
STAR and PHENIX collaborations; the STAR measurements are also inconsistent 
with the FONLL prediction by a factor of ~ 4. 

1.5.2 Interactions With the Medium 

A parton passing through a strongly interacting medium is expected to lose en- 
ergy to the medium through gluon bremsstrahlung. For heavy quarks, the gluon- 



bremsstrahlung distribution is estimated 69 to be reduced from the light-quark case 



by a factor of [1 + (Mq/Eq) ■ #~ 2 ]~ 2 , where M Q and E Q are the mass and energy of 
the heavy quark and 9 is the angle of gluon emission (with respect to the momentum 
of the heavy quark). The suppression factor is always less than 1, but is closest to 
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Figure 1.12: The nuclear modification factor for non-photonic e R-aa E * n the 0" 
10% most central Au + Au collisions measured by the PHENIX collaboration [55] 
and predicted by the DGLV 56 59 and BDMPS 60,61] models. The green dashed 



line (yellow band) shows the predictions of the BDMPS (DGLV) model, including 
only radiative energy loss. The green band shows the prediction of the DGLV model 
including both radiative and collisional energy loss. The thin dashed curves show the 
prediction of the DGLV model for e from D-meson decays only. 



zero for small angles (8 < Mq/Eq); this suppression of small-angle gluon radiation 
is called the "dead-cone effect." 

Because of the dead cone effect, it was expected that heavy quarks would lose 
less energy to the medium than light quarks. This would lead to less suppression 
of heavy quarks (and therefore of non-photonic e ) at high transverse momenta. 
However, measurements of the non-photonic nuclear modification factor |51[|52|55l 



70 at RHIC indicate that the suppression of non-photonic e at high px is similar in 



magnitude to the suppression of light-flavor hadrons 30 -32 . Several models of parton 
energy loss in a QGP have been developed; the models describe heavy quark energy 
loss through different mechanisms, including (but not limited to) medium-induced 
gluon radiation and elastic collisions with the partons that make up the medium. 
Some of these models will be described below and their predictions will be compared 
to the PHENIX collaboration's measurements of Raa and V2 for non-photonic e ± in 



central 200-GeV Au + Au collisions. 55 
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Figure 1.13: In the upper panel is the nuclear modification factor for non-photonic e 
Raa E m the 0-10% most central Au + Au collisions measured by the PHENIX col- 



laboration 55 and predicted by the Moore/Teany [62] and van Hees 63,64 models. 
The predictions of the Moore/Teany model are shown for two different choices of the 
heavy-quark diffusion coefficient. Shown in the lower panel is the elliptic flow parame- 
ter v 2 for non-photonic e measured by the PHENIX collaboration for minimum-bias 
(0-92% centrality) collisions and predicted by the two models. 



In its original form, the DGLV model 56 58 describes parton energy loss through 
medium-induced gluon radiation, with the parton interacting with multiple scattering 
centers as it passes through the medium. Perturbative QCD is used to calculate the 
matrix elements for gluon emission in these interactions. The opacity of the medium 



is parametrized by the gluon rapidity density dN g /dy. The BDMPS model 60 61 



describes parton energy loss through medium-induced gluon radiation. As it passes 
through the medium, the parton undergoes multiple independent interactions with 
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Figure 1.14: The nuclear modification factor for non- photonic e ± Raa in the °- 10 % 
most central Au + Au collisions measured by the PHENIX collaboration 55 and 
predicted by the Gossiaux/Aichelin 65,66 model. 



a series of scattering centers, which have screened Coulomb potentials (the weak- 
coupling regime of QCD). The magnitude of the energy loss in the medium is deter- 
mined by the color charge and mass of the parton, the path length in the medium, 
and the time-averaged squared momentum transfer from the parton to the medium, 
which is represented by the "transfer coefficient" q. The DLGV and BDMPS radiative 
energy-loss models with parameters dN g /dy « 1000 and 4GeV 2 /fm < q < 14GeV 2 /fm 
give values of Raa that reproduce the suppression of light-flavor hadrons well (see Fig- 
ure 1.9 and 30-32]). However, as shown in Figure 1.12 those models tend to under- 
predict the suppression of non-photonic e ± {i.e., over-predict the values of Raa) in 



heavy-ion collisions at RHIC. An extension of the DGLV model 59 to include parton 
energy loss through elastic scattering in the medium is in better agreement with the 
data, although the model may still under-predict the suppression at high p?. 



The model of Moore and Teany 62 describes heavy-quark energy loss and flow 
through elastic scattering in a thermalized, expanding QGP, with the scattering am- 
plitudes calculated perturbatively. The heavy-quark diffusion coefficient, which is 
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Figure 1.15: The nuclear modification factor for non-photonic Raa E i n the 0-10% 
most central Au + Au collisions measured by the PHENIX collaboration |55j and 
predicted by the Collisional Dissociation 67 68 model. 



related to the temperature of the medium, determines the magnitude of the energy 
loss and elliptic flow. To avoid infrared divergences in some of the perturbative cal- 
culations, an infrared regulator proportional to the gluon Debye massn is introduced. 



The model of van Hees et al. 63 64 is similar to the Moore/Teany model, but heavy 
quarks are also allowed to scatter through D- and 5-meson-like resonance states in 
the QGP. For the perturbative calculations of the elastic scattering amplitudes, a 
non-running coupling constant in the range 0.3 < a s < 0.5 is used and the gluon 
Debye mass is used as an infrared regulator. Calculations of the non-photonic e nu- 
clear modification factor and flow parameter v-i using these two models are shown in 



Figure 1.13 The Moore/Teany model is shown for two different values of the diffusion 
coefficient; this model does not appear to describe both Raa an d v 2 for non-photonic 
e simultaneously. The van Hees model describes Raa and i>2 reasonably well, al- 
though there are some indications that it too may under-predict the suppression of 
non-photonic e at high pj>. The uncertainty in the calculation is due in part to 
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The Debye screening length Ad is the length scale over which charges (color or electric) are 



screened in a plasma 71 ; the Debye mass mo is the inverse of this length. With c = h = 1, 



m D = A^, 1 = SlTy \ n c + \nf in a QCD plasma. 
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uncertainty in the choice of the strong coupling constant and the masses and widths 
of the D and B meson resonance states. 

Like the Moore/Teany and van Hees models, the model of Gossiaux and Aiche- 



lin 65 66] describes heavy-quark energy loss through elastic scattering in a QGP. 
However, this model uses a physical running coupling constant and a different infrared 
regulatoi0 in the perturbative QCD calculations of the scattering matrix elements. 
This model appears to describe the observed suppression of non-photonic e well for 



all but the most peripheral centrality classes (shown in Figure 1.14 for the 0-10% 
centrality class; see (55) for the other centrality classes). 

In the Collisional Dissociation model [67|[68] , heavy quarks lose energy through 
fragmentation to D and B mesons and the subsequent dissociation of those mesons 
through collisions with partons in the medium. The formation times for those mesons 
are expected to be less than the time required to travel through the medium and 
repeated cycles of fragmentation and dissociation lead to an effective energy loss. A 
calculation of the non-photonic e ± nuclear modification factor using this model is 



shown in Figure 1.15 . The Collisional Dissociation model appears to describe the 
high-p T suppression of non-photonic e ± . 



5 The infrared regulator in the Gossiaux/ Aichelin model is « 0.2m£>, a choice guided by a QED 
calculation j72j of muon energy loss in a plasma. 
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1.5.3 Relative Abundances of Charm and Bottom 
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Figure 1.16: Perturbative QCD (FONLL) calculations [54] of the contributions of D- 
and .B-meson decays to the non-photonic spectrum in 200-GeV p + p collisions. 



FONLL calculations [54] indicate that the bottom-decay contribution to the non- 
photonic e spectrum in nucleon-nucleon collisions dominates over the charm con- 



tribution at high transverse momenta (see Figure 1.16). The calculations indicate 
that the crossover occurs near pt ~ 5 GeV/ c, though this estimate has large un- 
certainties. The relative contribution of c- and 6-quark decays to the spectrum of 
non-photonic e has been measured using the azimuthal correlations of non-photonic 



e ± with hadrons. The method [73, 74 is illustrated in Figure 1.17 Similar to di- 



hadron correlation measurements (see Section 1.4.1), the difference in azimuth 
between a "trigger" non-photonic e and "associated" particles in the same event is 
calculated. Due to decay kinematics, the width of the near-side peak is expected to 
be wider when the e ± originates from a S-meson decay than when it originates from 
a D-meson decay. If a D and a B meson are traveling with the same momentum, 
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Figure 1.17: An illustration 
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of the fitting method used to find re, the ratio 
of non-photonic from B decays to the total non-photonic yield (from both 
D and B decays) using measurements of the azimuthal correlation between e ± and 
hadrons. In red (blue) is the D (B) contribution to the A(f> distribution, shown as 
black circles. The near-side peak in the A0 distribution (black circles) is fit with 
a linear combination (green) of the .D-decay (red) and -B-decay (blue) contributions 
calculated from simulations. For this pr bin (2.5GeV/c < px < 3.5GeV/c), tb = 
0.28 ±0.06. 



the D will be traveling with a larger Lorentz boost and its decay products will tend 
to be separated by smaller angles in the reference frame of the laboratory. The A<p 
distributions for trigger e from D and B decays are simulated using PYTHIA. [75] 
The near-side peak in the measured distribution is fit with a linear combination of the 
simulated D- and -B-decay distributions to extract r#, the fraction of non-photonic 
that come from B decays in 200-GeV p + p collisions. The results of this mea- 



surement are shown in Figure |1.18| the measured ratio is consistent with the FONLL 
prediction. 
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Measurements of the non-photonic e nuclear modification factor R AA anc ^ r B 
will allow constraints to be placed on the values of R e AA and R e AA , the nuclear modifica- 
tion factors for non-photonic from D and B decays, respectively. The non-photonic 
e ± nuclear modification factor can be written as 



= (1 - r B )RT A + r B R e f A . (1.15) 



Figure 1.19 shows the allowed values |76|| of R\° A and R e AA based on the STAR 
collaboration's measurements of R AA E for px > 5GeV/c in Au + Au collisions (the 
0-10% centrality class) and r# in p + p collisions. This measurement would exclude 
the predictions of the DGLV model (with only radiative energy loss) at the 90% 
confidence level. Measurements of R AA and vb with smaller uncertainties will allow 
more of the R AA vs. R e A B A parameter space to be excluded. However, it should be 
noted that the value of R AA E used in the calculation shown in Figure 1.19 is suspect, 
as the spectra of non-photonic e ± in p+p and Au + Au [77] collisions used to calculate 
it are subject to an erratumj^] The figure is included only as an illustration of the 
method. 



6 This erratum 78 indicates that both the p + p and Au + Au spectra decreased from the 
published values by a factor of ss 2. Therefore, the corrected values of Raa are the same as the 
original values withi n uncertainties, although the uncertainties changed. The calculations used to 
produce Figure 1.19 have not yet been redone. 
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Figure 1.18: The fraction (tb) of non-photonic e ± that come from B- meson decays 



as a function of px predicted [54] by FONLL calculations and measured 74 by the 



STAR collaboration through e -hadron correlation studies. The open circle indicates 
the measurement of tb from an e -D° correlation measurement 79 . 
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Figure 1.19: An illustration [76] of the method which can be used to constrain the 
possible values of the nuclear modification factor for non-photonic e ± from D and 
B decays {R e f A and R e AA ) based on the measurement of the R^a E an d r B- The 
predictions of some of the models discussed in Section 1.5.2 are given ("Djodjevic et 
al." refers to DGLV, "Adil et al." refers to Collisional Dissociation). The inset shows 
the total non-photonic e ± Raa as a function of pt- 77,78 
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1.6 Remarks on Non-Photonic e 
in 200-GeV Cu+Cu Collisions 

In this dissertation, the yields of non-photonic e in Cu + Cu collisions at ^s^n = 
200 GeV for the 0-20% and 20-60% centrality classes (as well as the overlap 0-60% 
centrality class) are measured. The non-photonic e ± nuclear modification factor in 
Cu + Cu collisions Rq^Cu are calculated by comparing those yields to measurements 
of the non-photonic e ± spectrum in 200-GeV p + p collisions by the STAR [80 1 and 
PHENIX 51 collaborations]^] The Cu + Cu centrality classes studied in this analysis 



are expected to have values of {N part ) similar to the 40-60% centrality class of Au + 
Au collisions. A measurement of R^a E * n Cu + Cu and Au + Au (centrality 40-60%) 
collisions allows the suppression to be studied in two collision systems with similar 
QGP volumes but different collision geometries: the overlap region in a mid-central 
Au + Au collision should have a larger eccentricity than the overlap region in a central 
Cu + Cu collision. If a difference were to be observed between Rq^gL an d -RauAu' ^ 
could be an indication that the physical shape of the QGP plays a role in determining 
the suppression of non-photonic e ± . 



7 The STAR collaboration's measurement is a new measurement not subject to the erratum 
mentioned in the previous section. 
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Chapter 2 
Experiment 



This chapter describes the experimental facilities used to collide heavy ions and to 
record data on the products of those collisions. The Relativistic Heavy Ion Collider is 
described in Section |2~T , followed by an overview of the STAR detector in Section 2.2 



Sections 2^3 and 2A describe the STAR Time Projection Chamber and the Barrel 
Electromagnetic Calorimeter, the two components of the STAR detector that are 
used most directly in the analysis presented in the following chapters. 



2.1 RHIC 



2.1.1 Overview of RHIC 



The Relativistic Heavy Ion Collider (RHIC) (8182 is a collider /accelerator facility at 
Brookhaven National Laboratory in Upton, NY. It was conceived in 1983 to advance 
the study of nuclear physics at high energies. RHIC became operational in 2000 and 
until recently was the highest energy heavy-ion collider in the world. The facility 



(see Figure 2.1) consists of two rings that carry beams of ions ranging in mass from 
protons to gold nuclei. The maximum possible energy to which ions can be accelerated 
decreases as their charge-to-mass ratio increases: for example, the maximum energy 
is 250 GeV for protons, 114.9 GeV per nucleon for Cu ions, and 100 GeV per nucleon 



for Au ions. 83 RHIC is capable of accelerating polarized beams of protons, allowing 
the spin structure of the proton to be studied; RHIC is the highest energy polarized- 
proton accelerator in the world. 

The main part of the collider consists of two rings 3.8 km in circumference, which 
intersect in six regions equally spaced around the ring. Each ring consists of six 
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Figure 2.1: The Relativistic Heavy Ion Collider. 81 



curved sections of length ~ 356 m and six straight sections (through the interaction 
regions) of length ~ 277 m. In each of the curved sections, each ion beam passes 
through a series of 32 superconducting dipole magnets with field strengths of 3.458 
T at the maximum beam energy, which bend the beams in the horizontal plane. In 
the curved sections, the beams are horizontally separated by 90 cm, and each beam 
passes through a separate set of magnets. It is therefore possible for each ring to carry 
a different species of ion; this capability has been used to collide beams of deuterons 
with Au ions at a center-of-mass energy per nucleon-nucleon pair of y/s^N = 200 GeV. 
Each curved section also contains quadrupole and sextupole magnets (and magnets 
with octupole and decupole components) [84] to focus the beams and correct for 
imperfections in the fields of the other magnets. In the areas where the straight and 
curved sections meet, a set of dipole magnets directs the two beams into a single 
beam line to achieve head-on collisions. 

The journey of heavy ions (Au ions in this example) begins at a pulsed sputter 
ion source, which produces negative ions with charge Q — — 1. Those negative ions 
are accelerated through a Tandem Van de Graaff accelerator, reaching a potential of 
14 MV at the center of the Tandem. There, the ions pass through a foil which strips 



36 



electrons from them, leaving Au ions with a charge of Q = +12. Those positive ions 
are accelerated through a potential of 14 MV, leaving the Tandem Van de Graaff 
accelerator with an energy of 1 MeV/u. A second stripping foil at the exit of the 
Tandem leaves Au ions with a charge state Q = +32 upon entering the transfer line 
for injection into the Booster synchrotron. There are two independent Tandem Van 
de Graaff accelerators, allowing for the production of beams of different species of 
ions. A separate linear accelerator is used to produce beams of protons, which are 
injected directly into the Booster synchrotron. In the Booster, RF cavities are used 
to group the ions into bunches and to accelerate them to energies of 95 MeV/u. Upon 
exiting the Booster, the Au ions pass through another stripping foil, which removes 
all but the two most tightly bound (i^-shell) electrons, leaving the ions with a charge 
state of Q = +77 upon injection into the Alternating Gradient Synchrotron (AGS). 
There, the ions are accelerated to energies of 8.86 GeV/u, then stripped of their final 
two electrons (leaving them in a charge state of Q = +79). Bunches of ions are then 
injected into one of the two RHIC rings; the ion beams in these rings circulate in 
opposite directions. Upon entering the RHIC rings there are ~ 10 9 ions per bunch 
(~ 10 10 protons per bunch when proton beams are being collided). In the RHIC rings, 
the acceleration RF system uses two RF cavities (per ring) operating at a frequency 
of 28.15 MHz and per cavity voltages of 20-300 kV to accelerate ions to the desired 
energy. A second "storage" RF system operating at a frequency of 197 MHz and a 
voltage of 6 MV is used to prevent the bunches of ions from spreading out in the 
longitudinal direction. The RF systems maintain 360 "buckets" around each ring, 
of which 120 buckets are typically filled with ion bunches. When RHIC is filled, it 
typically contains ~ 1.2 x 10 11 ions in each ring. 



2.1.2 The RHIC Experiments 

Ions may collide in the six interaction regions around the RHIC ring. Four of those 
interaction regions have hosted experimental detectors. The Solenoidal Tracker At 
RHIC (STAR) detector, which was used to collect the data for this dissertation, will 



be described in the following section. The PHOBOS detector 85 was designed to 
measure charged-particle multiplicities over a large solid angle using a set of silicon de- 
tectors around the beam line, covering the full azimuth and pseudorapidity \r]\ < 5.4. 
The PHOBOS detector also included two "arms" opposite each other in azimuth 
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and spanning the pseudorapidity range < rj < 2, with a width of 0.2 rad in az- 
imuth. These arms used several layers of silicon detectors and time-of-flight detectors 
to identify and measure the trajectories of particles with low momenta (transverse 
momenta down to 30 MeV/c could be accurately measured). The BRAHMS detec- 



tor 86 consisted of two moveable magnetic spectrometer arms covering small solid 
angles. The "Forward Spectrometer," which covered 0.8 msr and could be rotated 
through the polar angle range 2.3° < 9 < 30°, used a combination of dipole magnets, 
time projection chambers, drift chambers, time-of-flight hodoscopes, and Cherenkov 
detectors to track and identify particles produced with small angles relative to the 
beam line. The "Mid-Rapidity Spectrometer," which covered 6.5 msr and could be 
rotated through the polar angle range 30° < 9 < 95°, used a dipole magnet, time pro- 
jection chambers and a time-of-flight hodoscope to track and identify low-momentum 
particles produced at mid-rapidity. BRAHMS also included a barrel of silicon-strip 
and scintillator-tile detectors to measure charged particle multiplicities and collision 
vertex position around the nominal interaction point, covering the full azimuth and 
the pseudorapidity range \rj\ < 2.2. Both PHOBOS and BRAHMS have been decom- 
missioned. 

The detector for the Pioneering High Energy Nuclear Interaction experiment 
(PHENIX) [87] is a large, multi-component detector with four arms, two for par- 
ticle tracking and identification at mid-rapidity, and two used to track muons at large 



rapidities. The central arms 88-90 , which span the pseudorapidity range |^| < 0.35 



and two back-to-back sectors of width 90° in azimuth, consist of a set of multi-wire 
chambers, time-of-flight detectors, Cherenkov detectors, and calorimeters that are 



used to track and identify particles in an axial magnetic field. 91 The muon arms [92] , 
which cover the full azimuth and the pseudorapidity ranges — 1.15 > r\ > —2.25 and 
1.15 < rj < 2.44, consist of three multi-plane drift chambers that track charged 
particles in a radial magnetic field [9l] and a multi-layer muon identifier (alternating 
layers of steel absorbers and streamer tubes) to distinguish muons from other charged 



particles. PHENIX also includes a two- layer silicon strip detector 93 close to the 
beam, covering the full azimuth and the pseudorapidity range \r]\ < 2.6 to measure 
charged-particle multiplicities, collision vertex position, and the orientation of the 
event reaction plane. 



The Zero-Degree Calorimeters (ZDCs) [94 95 are a set of eight identical hadronic 
calorimeters, with one placed near the beam on each side of each of the four RHIC 
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experiments. The ZDCs act as minimum-bias triggers and luminosity monitors for 
heavy-ion collisions, and measure event characteristics in ways that are the same for 
each of the four experiments. The ZDCs are described in Section 2.2.1 



2.2 STAR 
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Figure 2.2: A cross-sectional view (in the yz-plane) of the STAR detector. 96 Note 
that the BEMC is labeled "EMC" here, the EEMC is not shown, and the SSD is not 
labeled. 



The STAR detector |96] is a large, multi-component detector primarily designed 
to measure hadron production over a large solid angle while operating in an envi- 
ronment with a high multiplicity of charged particles (~ 1000 charged particles per 
unit pseudorapidity for central Au + Au collisions). The major components of the 



STAR detector (see Figure 2.2) will be briefly described in this section. The analysis 
described in this dissertation is most directly concerned with the Time Projection 
Chamber (TPC) and the Barrel Electromagnetic Calorimeter (BEMC); those detec- 
tors will be described in greater detail in subsequent sections. The primary tracking 
component of the STAR detector is the Time Projection Chamber (TPC), a large 
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gas detector that sits inside a room-temperature solenoidal magnet that produces 
a uniform magnetic field of 0.5 T. Silicon detectors provide particle tracking near 
the nominal interaction point, while electromagnetic calorimeters sit outside of the 
TPC (but within the magnet) to provide energy measurements and assist in particle 
identification. 



2.2.1 Overview of STAR 



The TPC 97 is a large-volume, cylindrical gas detector with multi-wire proportional 
chambers as endcaps. It extends from cylindrical radius 50 cm to radius 200 cm in 
the range \z\ < 210 cm, covering the full azimuth. At its outer (inner) radius it covers 
the pseudorapidity range \t)\ < 0.9 (\rj\ < 1.8). The TPC will be described in greater 
detail in Section l2~3l 

Outside the TPC is the Barrel Electromagnetic Calorimeter (BEMC) [98], a large 
lead-scintillator sampling calorimeter that extends from cylindrical radius 223 cm to 
263 cm, covering the pseudorapidity range \rj\ < 1. The BEMC is ~ 20 radiation 
lengths deep, allowing for nearly full containment of electromagnetic showers. A 
sub-detector of the BEMC, the Barrel Shower Maximum Detector, sits 5-6 radiation 
lengths from the inner face of the BEMC and provides measurements of the shapes 
and sizes of electromagnetic showers. The BEMC will be described in greater detail 
in Section EH 

The TPC, BEMC, and most other components of the detector sit inside a magnetic 
field |99| produced by a solenoidal magnet of inner diameter 5.3 m and length 6.85 
m. This magnet produces a uniform magnetic field parallel to the beam line. The 
field strength may be adjusted within the range 0.25 T < B < 0.5 T; the polarity 
may be changed so that the field points along either the +z or — z directions. The 
magnet is typically operated with the maximum field strength, as was the case when 
the data analyzed in this dissertation were recorded. The field is produced by 10 main 
coils, which are connected in series and draw a current of 4500 A when producing 
the maximum field. The magnet also has four additional coils, called "trim coils," 
two at each end of the solenoid, that are used to improve field uniformity. The 
STAR magnetic field has been mapped using an array of Hall probes; the azimuthal 
component of the magnetic field is observed to have a magnitude of < 3 gauss, while 
the radial component is observed to have a magnitude < 50 gauss at the outer radius 
of the TPC (with even smaller deviations from closer to the center of the detector). 
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Between the beam pipe and the inner boundary of the TPC were two silicon 
detectors (the SVT and the SSD), which were present in STAR at the time the data 
analyzed in this dissertation were recorded, but were removed a few years ago. Closest 
to the beam pipe was the Silicon Vertex Tracker (SVT) |100 , which consisted of 216 
silicon drift detectors arranged in three concentric cylindrical barrels with radii of 6.9 
cm, 10.8 cm, and 14.5 cm. The SVT was segmented into 13 million pixels, providing 
20-/zm spatial resolution and energy-loss resolution of ~ 7%. The SVT contributed 
about 0.06 radiation lengths to the total material budget inside the TPC. The SVT 
was intended to provide precise tracking in the inner part of the STAR detector, 
allowing the position of the primary collision vertex to be reconstructed accurately 
and allowing the decay vertices of short-lived particles (such as strange and multi- 
strange baryons) to be reconstructed and distinguished from the collision vertex. 

Between the SVT and the inner boundary of the TPC, at a radius of 23 cm, was 
a fourth layer of silicon, the Silicon Strip Detector (SSD). [101] The SSD covered 
the pseudorapidity range \r]\ < 1.2 and consisted of ~ 500,000 silicon strip readout 
channels, providing spatial resolution of 740 /zm in the ^-direction and 20 jum in 
the transverse direction. By providing a fourth layer of position and energy-loss 
measurements, the SSD was intended to improve tracking in the inner portion of the 
STAR detector and improve the extrapolation of TPC tracks to hits in the SVT. 

The Forward Time Projection Chambers (FTPCs) |102 are two cylindrical gas 
detectors which surround the beam pipe and span the pseudorapidity ranges 2.5 < 
1 97 1 < 4, ranges not covered by the TPC. While the electric field in the TPC points in 
the indirection, the FTPCs have a radial electric field. Signals are read out in five 
ring-shaped multi-wire chambers on the outer surface of each cylinder. 



A single Endcap Electromagnetic Calorimeter (EEMC) |103 sits beyond the west 
end of the TPC. This annular lead-scintillator sampling calorimeter covers the full 
azimuth and the pseudorapidity range 1 < rj < 2, a range not covered by the BEMC. 
Similar to the BEMC, the EEMC is approximately 20 radiation lengths deep and has 
a shower maximum detector at a depth of about 5 radiation lengths. The Endcap 
Shower Maximum Detector (ESMD) is composed of scintillator strips (cf. the Barrel 



Shower Maximum Detector, which is a set of wire chambers, see section 2.4.2). 



The Zero-Degree Calorimeters (ZDCs) [94 95 are hadronic calorimeters which 
sit along the z-axis ±18 m from the center of the STAR detector and subtend an 
angle of 2.5 mrad. The ZDCs are located beyond the RHIC dipole magnets which 
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bend the beams to merge them into a single beam line through the interaction region. 
The ion beams and charged collision fragments are bent away from the ZDCs by the 
dipole magnets, while neutral collision fragments (mainly neutrons) strike the ZDCs, 
producing hadronic showers. Each ZDC consists of alternating layers of tungsten 
absorbers and wavelength-shifting fibers, longitudinally segmented into three mod- 
ules. The fibers collect the Cherenkov radiation produced as the shower particles pass 
through the calorimeter and direct that light to photomultiplier tubes for conversion 
into an amplified electronic signal. The STAR collaboration uses a coincidence in the 
ZDCs as a minimum-bias trigger for heavy-ion collisions, with the signal in each ZDC 
required to be greater than 40% of a single neutron signal. A precise measurement 
of the time difference between ZDC signals after a collision allows for a measurement 
of the z-position of the collision vertexj^J Identical ZDCs sit in similar positions near 
each of the other RHIC experiments. 



2.3 The TPC 




2.3.1 Overview of the TPC 

The STAR Time Projection Chamber (TPC) |97] is the primary tracking component 
of the STAR detector. It allows the trajectory of a charged particle moving through 
the STAR magnetic field to be recorded in three dimensions, permitting the momen- 
tum of the track to be reconstructed. The TPC is also used to measure the ionization 
energy loss (typically denoted dE/dx) of the particle, information which is useful in 
determining the identity of the particle. 

The TPC, shown in Figure 2.3, is a large volume of P10 gas (a mixture of 10% 



methane and 90% argon) contained between two concentric cylindrical "field cages" 
(50 cm and 200 cm from the beamline) and two endcaps (at z = ±210 cm). Within 
this volume, a uniform 140-kV/cm electric field is maintained. A thin, conductive 
central membrane (cathode) sits at z = at a potential of 28 kV, while planes of 
anode wires at the endcaps are held at ground. Field uniformity between these planes 
is maintained by the inner and outer field cagesQ the structures of which include a 
series of closely spaced equipotential rings. The rings are held at potentials that 
decrease linearly from the central-membrane potential to ground. 

The endcaps of the TPC are multi-wire proportional chambers with pad read- 
out. The endcaps consist of three planes of wires (the anode grid, the shielding grid, 
and the gating grid) and a plane of readout pads. When a charged particle passes 
through the TPC, it ionizes the gas. The ionization electrons drift parallel to the 
electric field towards the endcaps with a drift velocity of ~ 5cm/^s. The transverse 
(longitudial) diffusion of the drifting cloud of electrons is 230/zm/ y/cm (360/im/ y/cm). 
The fact that the electric and magnetic fields in the TPC are parallel helps to reduce 
the transverse diffusion. The electron cloud drifts to the anode-wire plane, where 
the large electric field near each wire accelerates the electrons enough to induce sec- 
ondary ionization, resulting in an avalanche of charge and amplification by a factor 
of 1000-3000. The positive ions produced in the avalanche induce an image charge 



on the nearby readout pads; the signal on each pad is then measured 104 by a 
preamplifier/shaper/waveform digitizer system. The ionization track produced by 
the passage of a charged particle is thus broken into a set of discrete "clusters." The 
x- and ^/-positions of each cluster are determined by the position(s) of the affected 
readout pad(s). The z-position is reconstructed by combining measurements of the 
drift velocity and the time at which the cluster signal was recorded. 



2 Thc field cages serve two purposes: containing the P10 drift gas and maintaining field uniformity. 
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Each endcap is divided azimuthally into twelve sectors; each sector is divided 
into inner and outer sub-sectors. There is a 3-mm gap between sectors. Within each 
sector, the wire planes and rows of readout pads run roughly in the azimuthal direction 
(i.e., perpendicular to the radial direction at the center of each sector). There are 
45 padrows in each sector. The pad width along the wire direction is 2.85 mm (6.20 
mm) for the inner (outer) sub-sectors. With these pad widths, the image charge 
from a single avalanche is spread out over several adjacent pads; measurements of the 
relative signal amplitude in each pad allow the centroid of a cluster to be determined 
with a spatial resolution less than the width of a single pad. 

The function of the shielding grid, which sits 2 mm into the gas volume from the 
anode grid, is to shape the electric field as it makes the transition from the avalanche 
region to the uniform field of the drift volume. The gating grid, which sits 6 mm 
into the gas volume from the shielding grid, acts as a "shutter" to control the entry 
of electrons from the drift volume and to prevent positive ions from the avalanche 
region from entering the drift region (where they could distort the electric field). 
In its "open" configuration, the wires of the gating grid are held at 110 V. In the 
"closed" configuration, alternate wires are held at potentials of 35 V and 185 V; the 
150- V potential difference between adjacent wires creates an electric field roughly 



perpendicular to the z (drift) direction. The gating grid is opened when a trigger 95 



signal is received; the 2.1-/is combined response time of the trigger and gating grid 
reduce the effective length of the TPC's tracking volume by 12 cm. 



The STAR tracking algorithm 105 reconstructs charged-particle tracks starting 
with clusters in the outermost few padrows of the TPC (where the track density is 
lowest). The algorithm works inward, finding nearby clusters and then sets of three 
clusters that lie along a straight line. As the track length increases, the algorithm 
begins fitting the track with a helix and adding clusters that fall near that helix. The 
full track is reconstructed once the algorithm reaches the innermost padrows. If it 
is assumed that the particle that produced the track had charge ±e, the transverse 
momentum can be determined [3] from the magnetic field strength B and the radius 
of curvature r of the helix: = \eB\r. The sign of the charge can be determined 
from the directions of the track's curvature and the magnetic field. The momentum 
of the track is determined from the measured transverse momentum and polar angle. 

The location of the collision (the "primary vertex" ) is found by extrapolating all 
tracks towards the center of the detector and finding the point V for which \ 2 — 
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Figure 2.4: The resolution of the position in the xy-pl&ne of the primary vertex; the 



resolution improves with increasing charged-particle multiplicity. 97 
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Figure 2.5: Relative pr resolution in the TPC for a 0.25-T magnetic field. 97 



^2j{dj/(Tj) 2 is at a minimum, where dj is the distance between the point V and the 
jth track helix and <jj is the uncertainty of that distance. 106 This procedure is 
repeated several times with large, outlying values of dj removed in each iteration. 
Figure 2.4 shows the resolution of the position of the primary vertex in the xy-pl&ne. 
The primary vertex resolution is less than 1 mm in the xy-pl&ne and in the ^-direction. 
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The position resolution for an individual cluster in the xy-plane and in the z- 



direction is < 4 mm. The transverse momentum resolution, shown in Figure 2.5 



is typically < 10%. The degradation in momentum resolution at low px is due to 
multiple Coulomb scattering^ while the degradation at high p T is due to the difficulty 
in determining the radius of curvature. 



2.3.2 Energy Loss 

Moderately relativistic charged particles (with 0.1 < f3j < 100) passing through mat- 
ter typically lose energy through ionization. The mean energy-loss rate is given fl7] 
by the Bethe-Bloch equation. For a particle with mass M and charge z traveling with 
speed (3c and Lorentz factor 7, 



/dE\_ e 4 P N A Z 2m e c 2 ^ 2 T max 2 5(^)\ 

\ dx I 4ire 2 m e c 2 A ' (3 2 \2 I 2 P 2 J' 1 ] 

where m e is the electron mass and Na is Avogadro's number. Z, A, and p are the 
atomic number, atomic mass, and density of the matter being traversed. T max is the 
maximum kinetic energy that can be imparted to a free electron in a single collision: 

2m e c 2 /3V r<? . 

max l + 2 7 m e /M+(m e /M) 2 - 1 ] 

The mean ionization potential per electron is / ~ Z-(10 eV). The quantity 5 accounts 



for the screening of the projectile's electric field in dense materials. 107 The only 



dependence on the projectile particle's mass enters Equation 2.1 through the factor 
of T max . The values of energy loss are distributed about this mean in a non-normal 
distribution, with a "tail" at high values of dE/dx. 

In the STAR TPC, energy loss is measured through the amount of charge collected 
for each TPC cluster. Ideally, the mean energy loss could be measured by calculating 
the average energy loss for all of the TPC clusters on a track. Since the energy loss is 
measured on a cluster-by-cluster basis, a track is effectively divided into segments of 
length ~ 2 cm (the width of a readout pad in the radial direction) for the purposes of 
energy-loss measurements. That segment length is too short to average out the large 
fluctuations in ionization energy loss (due to the "tail" in the dE/dx distribution). 



3 



The difference between pions and antiprotons is due to the fact that the RMS scattering angle 



is proportional to /3 1 17 and for a given momentum, /3 is larger for pions than for protons. 
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It was found that the energy loss could not be measured accurately by simply finding 
the mean of the cluster-by-cluster measurements. Instead, the STAR collaboration 
finds the most probable energy los^] by excluding the fraction (typically 30%) of 
clusters having the largest signals. This "truncated mean" value is found |108| through 
simulations to be a good estimate of the most probable value of the energy loss. The 
energy-loss resolution of the TPC is < 8% for tracks that cross 40 padrows. 

Instead of using the truncated mean energy loss directly, the STAR collaboration 



typically uses the normalized energy loss. |109| The measured energy-loss value (the 
truncated mean, which will be denoted as (dE/dx)) is compared to the expected 
value for a given particle species. The difference between the measured and expected 
values is divided by the energy-loss resolution. 

na x = in ±, (2.3) 

where Bx is the most probable value of (dE/dx) for particle species X and ax is the 
expected width of that distribution. For pions, the distribution of na w should be a 
Gaussian with mean and width 1 for all momenta (and the same should hold true 
for other particle species). A single value of (dE/dx) corresponds to multiple values 
of na, one for each particle species. For example, a value (dE/dx) = 2.5 keV/cm 
might correspond to na e = —5 and na n = 0. The values of Bx and ax are found 



through simulations. 108 For particle species X, the value of Bx as a function of 
momentum or velocity is called the "Bichsel function." The Bichsel functions have 
the same general shape as the Bethe-Bloch function and depend upon the number of 
track segments used in the calculation of (dE/dx) and the lengths of those segments. 



2.4 The BEMC 
2.4.1 BEMC Towers 

The STAR Barrel Electromagnetic Calorimeter (BEMC) [98] is a sampling calorimeter 
that sits outside of the TPC, covering the full azimuth and the pseudorapidity range 
\r)\ < 1. The BEMC is segmented into 4800 towers, 120 towers in <fi by 40 towers in 77; 
for each tower, A77 x A0 « 0.05 x 0.05. A cross-section of a two-tower BEMC module 

4 The most probable value is not the same as the mean due to the skewed shape of the energy- loss 
distribution. 
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Figure 2.6: A cross-sectional view (in the rry-plane) of a BEMC module (two towers 
wide). (98 



is shown in Figure |2.6| Each tower consists of 21 plastic scintillator tiles, most of 
which are 5 mm thickj^] separated by 5-mm-thick layers of lead absorber. Photons or 
e striking the inner face of the BEMC will initiate electromagnetic showers in the 
lead absorber layers. Hadrons striking the calorimeter will initiate hadronic showers, 
which develop over longer distance scales than electromagnetic showers; hadronic 
showers may have an electromagnetic component due to the production of neutral 

5 The first two scintillator layers are 6 mm thick to accomodate additional wavelength-shifting 
fibers for the Preshower Detector. 
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pions, which decay predominantly to photons. The charged particles in the shower 
interact with the scintillator, producing ultraviolet light. This light is collected by 
wavelength-shifting fibers which run in a a-shaped groove around the outside of each 
scintillator tile. These fibers guide the scintillation light out of the calorimeter and 
beyond the STAR magnet to a set of photomultiplier tubes (PMTs). The light from 
all 21 scintillator layers in a single tower is directed to a single PMT, which converts 
the light pulse to an electronic signal and amplifies it by a factor of ~ 10 5 . The 
wavelength-shifting fibers increase the wavelength of the scintillation light into a 
range to which the PMTs are sensitive. The signal from each PMT is integrated by a 
charge integrator and converted to a digital value by an Analog-to-Digital Converter 
(ADC). The ADC value is then converted to an energy using a linear conversion 
function. The relative energy resolution of the BEMC is 12% for an e ± with energy 
1.5 GeV; the energy resolution scales with l/\/E and is degraded by a few percent 
for events with high charged-particle multiplicity. 

The Preshower Detector is a sub-detector of the BEMC which is designed to 
measure energy deposition near the beginning of a shower and thus improve e ± -hadron 
discrimination. The first two scintillator layers each contain two wavelength-shifting 
fibers. The "extra" fibers direct scintillation light to illuminate one pixel of a multi- 
anode PMT, providing a separate measurement of the energy deposition in the first 
two scintillation layers. 



Electromagnetic showers typically 17 reach maximum energy deposition at a 
depth of 

t max ^[\n(E/E c )±0.5}X , (2.4) 

where E is the energy of the particle that initiated the shower, E c is the critical 
energy (7.43 MeV for lead), and Xq is the radiation length (637g/cm 2 = 0.5612 
cm for lead). The upper (lower) sign is used for showers initiated by photons (e ). 
The average depth needed for containment of 98% of a shower's energy is Lq ~ 
{tmax + 13.6)X . 110 For electrons with E = 5GeV/c, t max = 6X and L c = 



19.6Ao. The BEMC is approximately 20 radiation lengths (11 cm of lead) deep 
at 1] w 0, approximately the depth needed for near-full containment of an elec- 
tromagnetic shower in the energy range of interest for many STAR measurements 



(lGeV< E < 10 GeV). In contrast, one parametrization 110 of the longitudi- 



nal development of hadronic showers estimates the shower maximum to occur at a 
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depth of X m ax ~ [0-2 In (£7/1 GeV) + 0.7] A/, where Xj is the nuclear interaction length 
(199.6g/cm 2 = 17.59 cm for lead). The depth for 95% containment of a hadronic 
shower is Ac ~ Xmax + 2.5Xi(E/l GeV) 0,13 . Since the depth of the calorimeter is 
less than one interaction length, hadronic showers are not well contained within the 
BEMC. 

2.4.2 Shower Maximum Detector 



&\ Electromagnetic 




Front plane 



Figure 2.7: A schematic diagram of a portion of the Barrel Shower Maximum Detector 



(BSMD). 98 



The Barrel Shower Maximum Detector (BSMD) is a sub-detector of the BEMC 
which sits 5-6 radiation lengths from the inner face of the calorimeter, approximately 
the location of the maximum of an electromagnetic shower. The BSMD allows the 
shape and size of a shower to be measured with greater resolution than the BEMC 



towers. Figure 2.7 shows a schematic diagram of the BSMD. The BSMD is two layers 
of gas wire proportional counters with strip readout. Two layers of wires separated 
by an aluminum extrusion run in the z-direction, with 15 wires running over the 
width of a tower. In an electromagnetic shower, will cause avalanches (with a gas 
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amplification factor of ~ 3000) in the high electric field around the anode wires, this 
induces an image charge on the cathode (readout) strips. There are 36,000 cathode 
strips in two planes. The (inner) BSMD-r? plane is divided into 300 strips in the z- 
direction (and 60 strips azimuthally), allowing a measurement of the pseudorapidity 
(or, equivalently, z or 9). The (outer) BSMD-0 plane is divided into 900 strips 
azimuthally (and 20 strips in the ^-direction) allowing for a measurement of the 
azimuth. The strip widths and pitches are all between 1 and 2 cm (see [98] for the 
exact values), a small fraction of the dimensions of a tower (which is approximately 
12 cm x 12 cm). In each plane, adjacent strips that record a signal are grouped 
into "clusters." The centroid of a cluster can be determined with a resolution smaller 
than the width of a strip. The position resolution of the BSMD, measured using test 
beams in the AGS, is 2.4 mm + 5.6 mm /^E/lGeV (3.2 mm + 5.8 mm/^/1 GeV) 
for the rj (0) plane. 
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Chapter 3 

Analysis of Non-Photonic e^ 1 



This chapter describes the method used to extract the yield of non-photonic from 
Cu + Cu collisions at y/s^N = 200 GeV. Section 3.1 describes the criteria used to 



select the events from which the e yields will be extracted. Section 3.2 describes the 
selection cuts used to identify e ± in those events and section 3.3 describes the method 



used to extract the yield of photonic e . Section 3^ summarizes the correction factors 
that are applied to the measured data to find the fully corrected non-photonic e ± 
spectra; those correction factors will be discussed in detail in subsequent chapters. 



Table |3.1| summarizes the selection criteria and cuts used in this analysis. The 
symbols will be defined in the following sections. Note that the high-tower trigger 
condition is not always applied. 
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Table 3.1: Summary of cuts used in this analysis. See the following sections for more 
details. 



System/Measurement 


Condition 


Run Selection Criteria 


run number 


6031103 < run number < 6065058 


STAR detector 


TPC and BEMC included in run 


BEMC geometrical acceptance 


Abemc > 0.45 


Event Selection Criteria 


ZDC 


coincidence in both STAR ZDCs 


trigger ID 


trigger id = 66007 (MB) and/or 66203 (HT) 


high-tower trigger 


E T o Wer >3.75GeV 


z position of primary vertex 


v z < 20 cm 


reference multiplicity 


refmult > 19 


Track Quality Cuts 


radius of first TPC point 


Ritpc <102 cm 


number of TPC fit points 


N T pcfit > 20 


TPC fit points ratio 


R FP > 0.52 


Non-Photonic e ± Selection Cuts 


DCA to primary vertex 


GDC A < 1.5 cm 


track pseudorapidity 


-0.1 < rj < 0.7 


e ± Selection Cuts 


transverse momentum 


p T > 2 GeV/c 


TPC energy loss 


— 1.5 < na e < 3.5 


energy in EMC 


p/E Tower < 2/c 


BSMD cluster displacement 


A 5 md < 0.02 


BSMD cluster sizes 


Nsmdij > 2 and N S MD<t, > 2 


BEMC status 


"good" for tower, BSMD-?], and BSMD-0 


Partner Track Selection Cuts 


TPC fit points ratio 


RFp(partner) > 0.52 


partner track px 


Pt (partner) > 0.3GeV/c 


TPC energy loss 


(dE/dx) (partner) > 2.8 keV/cm 


Pair Selection Cuts 


DCA 


DCA(e-e + ) < 1.5 cm 


invariant mass 


M inv {e-e+) < 150MeV/c 2 
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3.1 Event Selection 



This section describes the criteria used to select the events to be analyzed. The 
data analyzed in this dissertation were recorded from 31 January to 6 March 2005, 
when RHIC was colliding beams of Cu ions at a/saT/v = 200 GeV The STAR TPC 
and BEMC are both used to identify e ± in this analysis, so only runs during which 
both of those detectors were operational were analyzed. A run is a period of data 
recording usually lasting from a few minutes to a half hour, during which thousands of 
collisions are recorded. 200-GeV Cu + Cu data were recorded before 31 January, but 
the STAR BEMC was not properly calibrated during those runs. A small number of 
runs exhibited corruption in the BEMC data and were excluded from the analysis. For 
some runs, an abnormally low BEMC geometrical acceptance (Abemc) was observed; 



only runs with Abemc > 0.45 are included in this analysis. See Chapter |5.1| for a 
more detailed discussion of the BEMC geometrical acceptance. 

For the Cu + Cu collisions analyzed in this dissertation, event readout is triggered 
by a coincidence in the STAR Zero-Degree Calorimeters during a bunch crossing. This 
is called the "minimum-bias" trigger (abbreviated as "MB"). During some runs, an 
additional trigger condition was applied: the Barrel Electromagnetic Calorimeter is 
required to have at least one tower with energy above a threshold value, Exower > 3.75 
GeV. This is called the "high-tower" trigger (abbreviated as "HT"). The actual 
trigger threshold is not in units of energy, but in ADC counts: the amplified signal 
produced by the photomultiplier tube attached to a tower is integrated over time 
by a charge integrator and then converted to a digital value by an Analog-to-Digital 
Converter (ADC). [98] Since not all towers have the exact same calibration, a single 
ADC threshold corresponds to a range of energy values. Typically, only one tower per 
event satisfies the HT trigger condition. The HT trigger results in an enhancement 
of particle yields at high p?. In this dissertation, the spectra of e ± from the MB and 
HT data sets will be combined, with MB data used for px < 4 GeV/ c and HT data 
used for p T > 4 GeV/ c. 

For each event, it is required that a primary vertex was found and that the .ex- 
position of that vertex is within the range \v z \ < 20 cm. This cut is applied in order 
to reduce the photonic e ± background. The support structures for the Silicon Vertex 
Tracker begin at z ~ ±30 cm and extend outward in the direction of the beam line. 
Photons produced in Cu + Cu collisions with \v z \ > 30 cm therefore have a higher 
probability of converting to photonic e~e + pairs than photons produced closer to the 
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200-GeV Cu+Cu, real data, vertex. 
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Figure 3.1: The ^-position (v z ) of the primary collision vertex for minimum-bias 
(black) and high-tower-triggered (red) events. 
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Figure 3.2: The position of the primary collision vertex in the xy-plane for MB events 
that pass the event selection cuts. 



center of the STAR detector. Figure 3.1 shows the distribution of v z for the two 



trigger types. Figure 3.2 shows the position of the primary vertex in the xy-plane. 
No cut is applied to v x or v y \ this figure is included only for reference. 
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200-GeV Cu+Cu, real data, reference multiplicity 



MB, triggered 
MB, accepted 
HT, triggered 
HT, accepted 
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Figure 3.3: The distribution of reference multiplicity for MB and HT triggered events. 
The distributions for all triggered events (gray and magenta) and for the subset of 
those events that were accepted by the event selection cuts (black and red) are shown. 



In addition, the reference multiplicity of each event is used to group events into 
centrality classes. Table |3.2 gives the reference-multiplicity ranges that correspond 
to each of the three centrality classes used in this analysis |111 . Also given for 
each centrality class are the values of (Nf, in ) and (N part ) (the mean number of bi- 
nary nucleon-nucleon collisions and the mean number of nucleons participating in 
collisions, respectively) for each of these centrality classes. These values were found 
from simulations based on the Glauber model of nucleus-nucleus collisions, see Sec- 
tion 1.3.2 Figure 3.3 shows the reference-multiplicity distributions for MB (black) 



and HT (red) triggered events. Note that the high-tower-triggered data set is biased 
towards more central (higher reference multiplicity) events. 



Table |3.2| also gives the total number of events analyzed for each trigger type and 
centrality class. N^ t B indicates the number of minimum-bias events that passed the 
event selection cuts. For high-tower triggered events, three types of event totals are 
given. N^J (accepted) indicates the number of events that passed the selection cuts 
and were labeled as having satisfied the HT trigger condition. However, in about 
14% of those events, the tower that recorded the signal that fired the HT trigger is 
recorded as having energy. This is due to the fact that whether or not the HT 
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Table 3.2: Reference-multiplicity ranges, collision characteristics, and event totals for 
the centrality classes of 200-GeV Cu + Cu collisions used in this analysis. 



Centrality 


0-20% 


20-60% 


0-60% 


Ref. Mult. 


98 < ref mult 


19 < ref mult < 97 


19 < ref mult 


(N bin ) 


156.351-;- 


42.68™ 


80.41™ 


{-^part) 


86.841^ 


33.611^8 






3,708,138 


8,960,915 


12,669,053 


N^J (accepted) 


1,323,950 


1,240,857 


2,564,807 


N?J(true) 


1,139,315 


1,076,621 


2,215,936 


N^ V J (normalized) 


331,754,960 


868,834,688 


1,200,589,648 



trigger fires is determined by the raw data (ADC counts) read out from each tower. 
In contrast, the recorded energy values are determined using each tower's calibration 
equation. In the event of a bad calibration, or if a tower is later determined to 
have bad status]]] a tower which recorded an ADC count above the threshold may 
subsequently be assigned energy. N^(true) is the number of high-tower triggered 
events for which the tower energy was nonzero. 

During runs in which the HT trigger was active, the MB trigger was prescaled: 
only a small fraction (usually < 1%) of events that satisfy the MB trigger conditions 
were actually recorded. This was done to conserve data storage space. For each run, 
one out of every F presca i e events was recorded. The value of F presca i e , the prescale 
factor, was not constant, but rather decreased as the RHIC luminosity (collision 
rate) decreased. For this analysis, the values of the prescale factor are within the 
range 23 < F presca i e < 595, though for most runs 100 < F presca i e < 400. When the 
HT non-photonic e spectra are normalized, the correct number of events by which 
the particle yield must be divided is not N^J (accepted) or N^J (true). Rather, 
the correct number of events for normalization is the number of MB events that were 
observed (whether recorded or not) during the time that the HT events were recorded. 
This quantity is the sum over all runs (during which the HT trigger was active) of 
F pr escaie times the number of recorded MB events that passed the appropriate event- 
and centrality-selection cuts. 

N*J (normalized) = ^[F prescale (run) x N^ t B (run)} (3.1) 



The values of N^ V J (normalized) are also given in Table 3.2 



1 see Section 5.1 
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3.2 e Identification 

This section describes the cuts used to identify e ± . 

Several cuts are used to select good-quality TPC tracks. Longer tracks with 
more TPC hits will tend to have their momenta reconstructed and their energy loss 
measured more accurately than shorter tracks. It is required that the radial position 
(in STAR's cylindrical global coordinate system) of a track's first recorded TPC hit 
(Ritpc) be less than 102 cm (i.e., the first TPC hit is in the innermost 9 pad rows of 



the TPC). Figure 3.4 shows the distribution of Ritpc f° r tracks that have passed all 
other e ± identification cuts. For most of the tracks, the first TPC hit is in the first 
few pad rows. 

Cuts are also applied to the number of TPC fit points (Nxpcfit > 20) and the fit 
points ratio (Rfp > 0.52). The fit points ratio (Rfp) is the ratio N T p C fit/N T p C max, 
where N PP cmax is the maximum possible number of TPC fit points. The value of 
Nrpcmax can vary from track to track, but is usually equal to 44, the numbei^] of 
functioning pad rows in the TPC. Lower values of Nxpcmax are possible, such as 



when a track lies at the edges of a TPC sector. Figures 3.5, 3.6, and 3.7 show 
distributions of N TPC fit, N TPCma x, and R P p for 200-GeV Cu + Cu collisions; in all 
cases all other e identification cuts have been applied. 

Cuts on track pseudorapidity and the distance of closest approach of a track with 
the primary collision vertex (this quantity is referred to as the "global DC A" or 
GDC A) are intended to increase the fraction of non-photonic e in the e sample. 
The radiation length of the STAR detector increases with the absolute value of the 
pseudorapidity for photons originating at the center of the detector. This is because 
photons with larger \rj\ will pass through the layers of the Silicon Vertex Tracker 
(SVT) and other inner structures at larger angles. In addition, support structures 
for the SVT begin at z ±30 cm and extend outward in the direction of the beam 
line. Photons with higher \t]\ will have a greater chance of striking these structures 
and converting to produce background e~e + pairs. The combination of the cut on the 
z-position of the collision vertex (\v z \ < 20 cm, discussed above) and a cut on track 
pseudorapidity (\r}\ < 0.7) reduces the number of e ± produced by photon conversions 
that are included in the sample. Due to the fact that only one half of the Barrel 
Electromagnetic Calorimeter (BEMC) was functioning when the data analyzed in 



2 The TPC has 45 p ad rows 97 , but pad row 13, the innermost pad row of the outer sub-sectors, 
was masked out. [l!2| 
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MB Trigger, inclusive er, R 1Tpc Distributions, cent. 0-60% 
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Figure 3.4: The radial position (in cylindrical coordinates) of the first TPC hit for 
tracks that pass the other e identification cuts. Most of these tracks have the first hit 
in the first few TPC padrows and satisfy the R\tpc < 102 cm cut. The distributions 
shown are for MB triggered events in several pt bins; the distributions are similar for 
HT triggered events. The gap near Rytpc = 120 cm is due to a masked-out padrow 
in the TPC. 



this dissertation were recorded, an asymmetric pseudorapidity cut (—0.1 < i] < 0.7) 
is used. The functioning half of the BEMC covers the half of the STAR detector with 
z > in STAR's global coordinate system. A track with a slightly negative value 
of pseudorapidity originating from a primary vertex near v z = +20 cm can strike 
the BEMC at a positive z coordinate; this is the reason for the use of -0.1 as the 



lower bound on 77. Figure 3J3 shows the pseudorapidity distributions of e ± in real MB 
events and simulated data. An increase in the photonic e yield for larger values of 
[-77 1 is observed. 

A cut on the distance of closest approach between a track and the primary colli- 
sion vertex (GDC A < 1.5 cm) is applied in order to preferentially select non-photonic 
e from the decays of heavy-flavor hadrons. A D ± meson with momentum 2 GeV/c 



'0 



(20 GeV/c) produced at the collision vertex has a probability of only 3 x 10 (1.1 
of decaying farther than 1.5 cm from the vertex (based on mass and lifetime measure- 
ments in [17]); the other D and B mesons have even lower probabilities of decaying 
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far from the vertex. In contrast, a K with momentum 2 GeV/c (20 GeV/c) produced 
at the collision vertex has a probability of 99.9% (1 — 10 -4 ) of decaying farther than 
1.5 cm from the vertex. Photon-conversion e will not be produced in any significant 
amount within the radius of the beampipe (4 cm). However, photonic e~e + pairs 
from 7T° Dalitz decays and the decays of other neutral mesons will all originate from 
near the collision vertex and will not be excluded by the GDC A < 1.5 cm cut. It 
should be noted that even if an e originates far from the vertex, the helix of its 
track may still have a DC A of less than 1.5 cm with the vertex, a situation which 
becomes more common as transverse momentum increases and the track radius of 



curvature becomes larger. Figure 4.20 (page 100) shows a GDC A distribution for 
simulated non-photonic e and for photonic e from both real and simulation data. 
The distribution for non-photonic e ± is narrower than the distributions for photonic 
e , although all the distributions lie predominantly below GDC A = 1.5 cm. The 
simulations used to create these distributions will be discussed in Chapters [4] and [7j 
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MB Trigger, inclusive er, N Distributions, cent. 0-60% 
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Figure 3.5: The number of TPC fit points for tracks that pass the other e identi- 
fication cuts. Most of these tracks satisfy the N TPC fit > 20 cut. The distributions 
shown are for MB triggered events in several px bins; the distributions are similar for 
HT triggered events. 
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Figure 3.6: The maximum possible number of TPC fit points for tracks that pass the 
other e ± identification cuts. Most tracks have N T p Cma x near 44. 
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MB Trigger, inclusive er, R pp Distributions, cent. 0-60% 
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Figure 3.7: The ratio Rfp = Nxpcfit/NTPCmax for tracks that pass the other 
identification cuts. The tracks shown here all satisfy the Rfp > 0.52 cut. The 
distributions shown are for MB triggered events in several bins; the distributions 
are similar for HT triggered events. 



MB Trigger,^ Distributions, cent. 0-60%, 2.0<p t <3.0 



in 

§ 0.08 
o 

"S 0.07 

N 

E 0.06 



0.05 
0.04 
0.03 
0.02 
0.01 



* real inclusive e ± 

* real pho. e* 

* sim. single e ± 
O sim. pho. e ± (from 7i°) 




Figure 3.8: The pseudorapidity distributions for inclusive (in the range 2 GeV/c < 
Pt < 3GeV/c) in real events (black), photonic e ± in real events (red), simulated 
non-photonic e ± (green), and simulated photonic e ± from photon conversions (blue). 
The yield of photonic e ± is slightly larger for larger values of pseudorapidity. The 
distributions for real data are from MB triggered events. The method used to extract 
the photonic e ± distribution from the real data is described in Section 3.3 
simulations will be discussed in Chapters [4] and [7J 
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Figure 3.9: Normalized TPC energy loss (na e ) vs. particle momentum for tracks in 
MB events that have passed the track quality cuts; the e identification cuts using 
data from the BEMC have not been applied. The curves indicate the estimated most 
probable value of na e for each of five particle species. The horizontal lines show the 
range of the na e cut used in this analysis. 



Several quantities measured using the TPC and BEMC are used to distinguish 
e from hadrons. The single most important cut used to identify e is a cut on 



particle energy loss in the TPC. Figure 3J) shows na e , the normalized TPC energy 
loss, versus momentum for tracks that pass the track quality cuts; the e identification 



cuts using data from the BEMC have not been applied. Figures |3.10| and |3.11| show 
na e versus momentum for tracks that pass all other e identification cuts for MB and 



HT triggered events. See Section |2.3.2| for a more detailed discussion of na e . The 
curves in Figures 3.9, 3.10, and 3.11 show the estimated^ most probable values of na e 
for five different particle species. For e ± , the most probable values of na e are around 
0. For this analysis, a cut on the normalized TPC energy loss of —1.5 < na e < 3.5 
is applied to exclude particle species other than e . This cut on na e approximately 
corresponds to a cut of 3.35 keV/cm < (dE/dx) < 5 keV/cm, where (dE/dx) is the 



70% truncated mean value of the TPC energy loss (see Appendix [A]) . Figures 3.12 
3 See Appendix A for a description of the calculation of those curves. 
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and 3.13 show na e distributions in single transverse- momentum bins for MB and HT 
triggered events with various combinations of e selection cuts applied (the track- 
quality cuts have been applied in all cases). The BSMD cuts are observed to be more 
effective at removing hadron contamination than the p/Ej< ower cut. When the HT 
trigger condition is applied, the p/E Tower is largely redundant. 

For momenta less than 1.5 GeV/c, it becomes difficult for the TPC energy-loss 
cut to distinguish between e ± and charged hadrons. For this reason, and due to 
problems in the simulation data for photonic e (data which are used in calculating 
the background rejection efficiency Eb), this analysis is restricted to particles with 
p T > 2 GeV/c. 
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200-GeV Cu+Cu, MB trigger, na e for charged particles, p/E and SMD cuts, cent. 0-60% 
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Figure 3.10: Normalized TPC energy loss (na e ) vs. particle momentum for tracks 
in MB triggered events that have passed all other e identification cuts. The curves 
indicate the estimated most probable value of na e for each of five particle species. 
The horizontal lines show the range of the na e cut used in this analysis. 



200-GeV Cu+Cu, HT trigger, no,, for charged particles, p/E and SMD cuts, cent. 0-60% 




10 12 
p [GeV/c] 

Figure 3.11: Normalized TPC energy loss (na e ) vs. particle momentum for tracks in 
HT triggered events that have passed all other e ± identification cuts. 
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Figure 3.12: Normalized TPC energy loss (na e ) for particles with momenta in the 
range 2GeV/c < p < 2.5GeV/c in MB triggered events. The effects of various sets 
of cuts on the na e distribution are shown (the track-quality cuts have been applied 
in all cases). The black histogram shows the na e distribution for all particles (no 
e selection cuts), while the red histogram shows the na e distribution for hadrons 
(identified by requiring that p/E Tower > 2.5 and that N S MDr) = 1 or Nsmd4> = !)• 
Also shown are the na e distributions for particles that pass the p/ETower < 2 cut 
(orange), particles identified as e ± using the BSMD (cyan), and particles identified 
as e ± using all cuts (green). These BEMC-related cuts are described in the text. The 
BSMD cuts are observed to be more effective at removing hadron contamination than 
the p/E Tower cut. 
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200-GeV Cu+Cu, HT trigger, na e for charged particles, cent. 0-60%, 4.0<p<5.0 



10 3 



10 2 



10 




o all 

^ hadrons 

TOW 

SMD 
■ TOW && SMD 

no. cut 



-15 



-10 



10 15 



Figure 3.13: Normalized TPC energy loss (ner e ) for particles with momenta in the 



range 4 GeV/c < p < 5 GeV/c in HT triggered events. See Figure 3.12 (caption) for a 
description of each histogram. All tracks are required to point to a tower that satisfies 
the HT trigger condition, which renders the p/ETower cut redundant. 
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200-GeV Cu+Cu, MB trigger, p/E for e*, cent. 0-60% 

Tower 




Figure 3.14: The ratio p/E Tower vs. px for tracks in MB triggered events that pass 
all other e ± identification cuts. Most e ± have p/ETower ~ 1/c. 
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Figure 3.15: The ratio p/ETower vs. pt for tracks in MB triggered events that are 
identified as hadrons (requiring that na e < —1.5 and without BSMD cuts). Hadrons 
have a much broader p/E Tower distribution than do e ± . 

The STAR Barrel Electromagnetic Calorimeter is also used to distinguish from 
hadrons. The BEMC is approximately 20 radiation lengths deep, the depth required 



200-GeV Cu+Cu, HT trigger, p/E Tower for e ± , cent. 0-60% 




Figure 3.16: The ratio p/ETower vs. pt for tracks in HT triggered events that pass all 
other e identification cuts. Since all tracks were required to point to a tower that 
satisfied the HT trigger condition, the p/ETower distribution is heavily biased and the 
p/ETower < 2/c cut is largely redundant. 



to contain approximately 95% of the energy of an electromagnetic shower with energy 



< 20 GeV (see Section 2.4). In contrast, the BEMC is ~ 0.6 nuclear interaction 
lengths deep, much less than the depth (several interaction lengths) needed to contain 
a hadronic shower with energy < 20 GeV. Therefore, e ± will deposit, on average, more 
energy in the calorimeter than hadrons with the same momentum. 

In this analysis, a cut of pjE^ ower < 2/c is used, where p is the momentum of a 
track in the TPC and E^ ower is the energy measured in the BEMC tower to which 
that track projects. An additional requirement that Exower > is implied whenever 
the p/ETower cut is mentioned. It is also required that the tower to which the track 
projects has "good" status: i.e., it has been determined that the tower was operating 



properly (see Section 5.1). 

Figure 3.14 shows the ratio p/ETower versus transverse momentum for e candi- 
dates; these particles are mostly e ± , with some (a few percent) hadron contamination. 
A hadron that appears as contamination in this distribution must have had both large 
TPC energy loss and also large electromagnetic component early in the development 
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Figure 3.17: Distributions of p/ETower for particles with transverse momenta in the 
range 2 GeV/c < pt < 2.5 GeV/c in MB triggered events. The effects of various sets of 
cuts on the p/E Tower distribution are shown (the track-quality cuts have been applied 
in all cases). The black histogram shows the p/ETower distribution for all particles 
(no e ± selection cuts). The solid red histogram shows the p/ETower distribution for 
hadrons identified by requiring that n<j e < —1.5 (and without BSMD cuts), while 
the open red histogram shows the p/ETower distribution for hadrons that have passed 
the e identification cuts in the BSMD. Also shown are the p/ETower distributions 
for particles that pass the TPC energy- loss cut (magenta), particles identified as e 
using the BSMD (cyan), and particles identified as e ± using all cuts (blue). The 
BSMD cuts are described in the text. The TPC energy-loss cut is observed to be 
more effective at removing hadron contamination than the BSMD cuts. 



of its shower, leading to an increased likelihood that it would pass the BSMD cuts 



(see below). Figure 3.15 shows the ratio p/E ToW er versus pt for particles identified as 
hadrons by requiring that na e < —1.5. 

An e ± may also fail the p/ETower cut if it struck the calorimeter near the edge 
or corner of a tower and its energy was deposited in multiple towers. This would 
be expected to skew the p/E To wer distribution, producing a "tail" at large values of 
p/E To wer- Transverse to an electromagnetic shower's direction of propagation, 95% 
of the shower's energy is expected to be contained within a cylinder of diameter 



4pM, 107 where pu is the Moliere radius (1.602 cm for lead) |17|. The towers in 
the BEMC are approximately 12 cm x 12 cm (at rj = 0) while the expected shower 
diameter is about half of those dimensions, 6.4 cm. Most electromagnetic showers 
will therefore be contained within a single tower. Any losses in the efficiency of the 
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Erower/p cut due to the energy of a shower being shared between multiple towers are 
accounted for in the e ± reconstruction efficiency (sr, see Chapter [4]). 



Figure 3.15 indicates that many hadrons pass the p/Exower < 2/c cut. Those 
hadrons likely produced neutral pions in collisions near the beginning of their showers, 
leading to large electromagnetic components and more of their energy being deposited 



in the calorimeter. 107 In the high-tower triggered data set (but not the minimum- 
bias data set), the EMC tower to which a track projects is required to have satisfied 
the high-tower trigger condition. Figure 3.16 shows the distribution of p/E Tower 
versus transverse momentum for e in HT triggered events. The distribution is biased 



towards low values of p/Exower and the cut on that ratio has little effect. Figure 3.17 
shows p/Exower distributions in a single transverse- momentum bin for MB triggered 
events with various combinations of e ± selection cuts applied (the track-quality cuts 
have been applied in all cases). 
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MB Trigger, inclusive er, A SMD Distributions, cent. 0-60% 
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Figure 3.18: The BSMD cluster displacement A$md for tracks that pass the other 
e identification cuts. The peaks of the distributions are below A$md — 0.02. The 
distributions shown are for MB triggered events in several bins; the distributions 
are similar for HT triggered events. 



The STAR Barrel Shower Maximum Detector (BSMD, see Section 2.4.2 ) is used to 
measure the shape and size of electromagnetic showers with finer resolution than the 
sizes of the towers. The BSMD sits 5-6 radiation lengths inside the calorimeter, near 
the maximum density for electromagnetic shower with energies greater than about 



1-2 GeV p38l; for a 2-GeV shower one parametrization 17 predicts a maximum at 
[ln(E/7A3 MeV)— 0.5]X = 5.1X or 2.9 cm of lead. Hadronic showers are expected to 
develop more slowly on average, reaching a maximum around one nuclear interaction 



length inside a calorimeter; for a 2-GeV shower, one parametrization 110 predicts a 
maximum at [0.2 \n(E/l GeV) + 0.7]A/ = 0.83A/ = 26X in lead. At the depth of the 
BSMD, many hadronic showers will not be fully developed. 

The STAR collaboration uses an algorithm to find groups of adjacent strips (called 
clusters) for which signals were recorded. Each track in this analysis is projected to 
the location of the BSMD and the nearest clusters in the r\ and <ft planes of the BSMD 
are found. A cut is applied on the distance between the track projection and the two 
clusters. The BSMD cluster displacement A S md is defined as 
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MB trigger, SMD els. sizes for e , cent. 0-60%, 2.0<p r <2.5 




MB trigger, SMD els. sizes for If, cent. 0-60%, 2.0<p. r <2.5 




Figure 3.19: Cluster sizes in the rj and <p planes of the BSMD for tracks with transverse 
momenta in the range 2GeV/c < px < 2.5GeV/c in MB events. The center of the 
each black cross indicates the mean BSMD cluster sizes in each plane, while the bars 
indicate the standard deviation in the BSMD cluster sizes. The left plot shows Nsmd</> 
vs. Nsmdti for tracks that satisfy all other e identification cuts. The right plot shows 
NsMDcf) vs. Nsmdt) f° r tracks identified as hadrons by requiring that na e < —1.5. 



&smd = V {Ar] SMD f + (A(f) SMD ) 2 , (3.2) 

where ArjsMD is the difference in pseudorapidity between the track projection and 
the centroid of the nearest cluster in the BSMD-r? plane, and A<Psmd is the difference 
(in radians) of the azimuth between the track projection and the centroid of the 
nearest cluster in the BSMD-0 plane. It is required that A$md < 0.02, equivalent 



to a distance of m 4.6 cm. Figure [3718] shows distributions of A$md for particles in 
minimum-bias events that have passed all other e identification cuts. Most of these 
distributions lie well below A S md = 0.02. It is also required that the centroid of each 
cluster lies in a strip that had "good" status (was determined to have been working 
properly) . 

A cut on the sizes of the nearest cluster in each plane of the BSMD is applied. The 
size of a cluster is measured by the number of strips in which a signal was recorded. 
It is required that N S md v > 2 and N SMD( p > 2, where N SMDr} (N S md<i>) is the size 
of the nearest cluster in the rj (0) plane of the BSMD. Figure 3.19 shows Nsmd</> vs. 
Nsmdt] for particles in the transverse-momentum range 2GeV/c < pr < 2.5GeV/c 
for e (particles that have passed the other e identification cuts) and for hadrons 
(identified by requiring that na e < —1.5). Figure 3.20 shows the mean values of 



73 



MB trigger, BSMD mean cluster sizes for charged particles, cent. 0-60% 
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Figure 3.20: Mean cluster sizes in the r\ and planes of the BSMD versus transverse 
momentum for tracks in MB events. The dark (light) green histogram indicates 
(Nsmdti) ((NsMDtj>)) for tracks that pass all other e ± identification cuts. The red 
(magenta) histogram indicates {Nsmdti) ((Nsmd4>)) f° r tracks identified as hadrons 
by requiring that na e < —1.5. The vertical bars indicate the standard deviations 
of Nsmdt/ and Nsmd</>- It is observed that have higher mean cluster sizes than 
hadrons and that (NsMD-q) is consistently larger than (Nsmd4>) for a given particle 
type. 



Nsmdti and NsMD<f> for e and hadrons as functions of px- The e are observed to 
have larger values of (Nsmdti) and (Nsmd^)- A difference in the performance of the 
two planes of the BSMD is also observed: (Nsmdti) is observed to be larger (0.3-0.4 
strips) than (Nsmd</>)- This is due to the loss of the low-energy component of showers 
in the inner (77) plane of the BSMD, leading to reduced spatial resolution and lower 



energy deposition in the outer (0) plane. 113 



3.3 Photonic e Background Subtraction 

Once e have been identified, it is necessary to identify the large photonic e ± back- 
ground, which will be subtracted from the inclusive e ± yield to give the non-photonic 
e ± signal. Photonic are identified by reconstructing low- invariant-mass e~e + 
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pairs. Photon conversions tend to produce e e + pairs with invariant mass near 
0, and the pairs produced in neutral-pion Dalitz decays e + 7) have in- 

variant masses less than the tt° mass. In this analysis, a cut on the pair invariant 
mass of M inv Xe-e + ) < 150MeV/c 2 , slightly above the 7T° mass, is applied. Other 
sources (including rj and rf meson Dalitz decays) may also produce e~e + pairs with 
M inv .(e~e + ) < 150MeV/c 2 . Photonic e~e + pairs with invariant mass greater than the 
150MeV/c 2 cut are not identified using this method; the yield of photonic e from 
those sources ("residual photonic e ±n ) will be estimated in Chapter [8] and subtracted 
from the efficiency-corrected non-photonic e yield. 



200-GeV Cu+Cu, MB trigger, M (pair) distributions for pairs, cent. 0-60%, 2.0<p (e*)<3.0 



400 



300 



200 



100 




photonic 

photonic (from sum bkg. 



_±_ 



_±_ 



_l_ 



_L 



_J_ 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

NUpair) [GeV/c 2 ] 

Figure 3.21: Distributions of invariant mass for pairs of e ± candidates and partners 
in MB triggered events. The transverse momenta of the candidates are in the 
range 2 GeV/ c < pt < 3 GeV/ c; all other e identification, pair, and partner selection 
cuts have been applied. In black is the distribution for unlike-charge pairs (n un uk e )- 

The combinatorial background calculated as 2^/n n ++ (see Equation 3.3) is shown 

in blue, while an alternate calculation of the combinatorial background (n. + n ++ ) 

is shown in green. The red and gray histograms are the invariant-mass distributions 
for photonic e~e + pairs: the combinatorial background subtracted from the number 
of unlike-charge pairs. The invariant mass distribution for photonic e~e + pairs is 
peaked near 0. The alternate method of calculating the combinatorial background 



gives essentially the same result as Equation 3.3 
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200-GeV Cu+Cu, MB trigger, DCA(pair) distributions for pairs, cent. 0-60%, 2.0<p (e t )<3.0 
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Figure 3.22: Distributions of the distance of closest approach for pairs of candidates 
and partners in MB triggered events. The transverse momenta of the candidates 
are in the range 2 GeV/ c < px < 3 GeV/ c; all other identification, pair, and partner 



selection cuts have been applied. See the caption of Figure 3.21 for a definition of 
each histogram. Photonic e~e + pairs have pair DCA distributions peaked near 0, 
although most random pairs also have DCA near 0. 



The number of photonic e~e + pairs is found as follows. Each particle that has 
passed the e ± identification cuts described in the previous section (called an Li e ± 
candidate") is paired with "partner" tracks from the same event. The number of pairs 
of oppositely charged tracks for which Mi nv .(pair) < 150MeV/c 2 and DCA(pair) < 
1.5 cm is recorded. However, it is possible that a random combination of an e 
candidate with a partner track could satisfy these conditions. The combinatorial 
background is constructed from pairs in which both tracks have the same charge. The 
combinatorial background is then subtracted from the number of oppositely charged 
pairs to give the number of photonic e"e + pairs. Mathematically, this is expressed as 

P = n unHke - 2 v /n__n ++ , (3.3) 

where P is the number of true (non-random) photonic e~e + pairs and n un uk e is the 
number of oppositely charged pairs that pass the invariant-mass and DCA cuts. The 
quantity n ("<++) is the number of pairs of electrons (positrons) with negatively 
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(positively) charged partners that pass the invariant-mass and DC A cuts. It should 
be noted that when both members of a pair pass the e identification cuts (i.e., both 
tracks qualify as candidates) that pair is counted twice. Given this, and the fact 
that photonic e can have only one true conversion partner, the quantity P in Equa- 
tion 3.3 is the same as the number of photonic e ± (modulo efficiency corrections). 



Figure |3.21| illustrates the combinatorial background subtraction procedure using the 
invariant mass distribution. It is possible to calculate n un uk e , the combinatorial back- 



ground, and the photonic e yield as functions of other quantities. Figure 3.22 illus- 



trates the calculation of the DCA distribution for photonic e~e + pairs. Figure 6.3 
(page 111) shows those three quantities as functions of the e ± candidate's value^Jof 
na e . 

A few selection cuts are applied to the partner tracks. Simulations by the STAR 



collaboration 97 indicate that for transverse momenta less than a few hundred MeV, 
the pt resolution of the STAR detector is degraded due to multiple Coulomb scattering 
in the inner structures of the detector. The fractional difference between the simulated 
Pt and the transverse momentum reconstructed by the STAR detector is observed 
to increase sharply with decreasing transverse momentum for p? < 300MeV/c for 
low- mass particles (charged pions in the simulations). Therefore, in this analysis 
partner tracks are required to have pT{partner) > 300MeV/c. A track-quality cut 
is also applied to partner tracks: it is required that Rpp(partner) > 0.52 (the same 
fit-points-ratio cut used for e candidates). Finally, a TPC energy-loss cut is used to 
exclude some charged hadrons from the set of partner tracks and reduce the number 
of random pairs in the sample. It is required that (dE/dx) < 2.8 keV/cm, which 
generally corresponds to a normalized TPC energy loss of na e < 4, far enough from 
the e ± peak that the number of e ± excluded by the cut is expected to be negligible. 



3.4 Correction of Spectra 

The uncorrected non-photonic e yield is simply the difference between the number 
of e identified (the inclusive e yield) and the number of photonic e identified. 
A number of factors have been computed to correct for inefficiencies in the analysis 

4 When both members of a pair qualify as e candidates, the pair is counted twice: once at the 
value of na e of the first candidate and once at the value of na e of the second candidate. 
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method described above. Subsequent chapters describe the calculation of each correc- 



tion factor in detail; see Table |3.3| for a short summary of the correction factors and 
the chapters in which their calculations are discussed. Simulations of non-photonic e 
and their interactions with the STAR detector are used to find the e ± reconstruction 
efficiency {er), the efficiency with which tracks are reconstructed in the STAR 



TPC and then identified as e ± using the selection cuts described in Section 3J2 The 
mean BEMC geometrical acceptance ((Abemc)) is used to correct for temporary 
losses in BEMC geometrical acceptance. The e spectra in the high-tower triggered 
data set are corrected by the HT trigger efficiency et- The simulations used to cal- 
culate Er do not properly reproduce TPC energy loss; therefore, the efficiency of 
the TPC energy- loss cut (£dE/dx) is calculated separately using the distributions of 
na e in real data. To correct for the hadron contamination that remains after the 
application of the identification cuts, the purity of the inclusive sample (K inc ) 
is calculated, also using the the distributions of na e in real data. The background 
rejection efficiency (eb) is the efficiency with which photonic e are identified using 



the method described in Section |3.3| this efficiency is calculated using simulations of 



photonic e and their interactions with the STAR detector. 
The efficiency-corrected non-photonic e yield is 

Nec = 7i (K inc I - — P 1 , (3.4) 



\A-BEMC) e R " £ T ■ £dE/dx \ £ B 

where / is the uncorrected inclusive e yield and P is the uncorrected photonic e 



yield given in Equation 3.3 



The residual background is the sum of the estimated yields of e ± from non-open 
heavy-flavor sources, which are subtracted from the efficiency-corrected non-photonic 

yields. The components of the residual background are photonic e ± from light- 
flavor sources with Mi nv .(e~e + ) > 150MeV/c 2 (denoted Blfp), non-photonic e from 
light-flavor sources (mainly the semileptonic decays of kaons, denoted B>lfn), pho- 
tonic e ± from J ftp and T decays (Bj/^ and B r , respectively), and photonic e ± from 
the Drell-Yan process (Bd-y)- Most of the components of the residual background 
are estimated using simulations of particles decaying into e ; the Drell-Yan contribu- 



tion is calculated 114 using perturbative QCD. The yield of non-photonic e ± from 



open-heavy-flavor sources (mostly D and B meson decays) is 

Yd,b = Nec — B LF p — B LF n — Bj/^ — B r — B D _ Y - (3.5) 
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Table 3.3: A summary of the correction factors used in this analysis and the chapters 
or sections in which their calculations are described. 



Name 


Symbol 


Corrects for 


Chapter/ 
Section 


pi T-ppnnfitriiction eff 

v. X V_y V_y X X kj> U X I., 1. V. U XV/ XX V^XX • 


F E> 

c H 


track reconstruction 
e ± ID cuts 


a 


mean BEMC 
geometrical acceptance 


(Abemc) 


losses in BEMC 
geometrical acceptance 




5.1 




HT trigger eff. 


f t 


HT trigger "turn-on" 




5.2 




energy-loss cut eff. 


e dE/dx 


TPC energy- loss cut 




6 




inclusive e ± purity 


K 


hadron contamination 




6 




background rejection eff. 


£b 


photonic e ± ID 




7 




Residual Background Components 


residual photonic e ± 
from light-flavor sources 


Blfp 


photonic e ± from 
77, p°, w, V, 0, and Kg 




8.1 




residual non-photonic e ± 


Blfn 


non-photonic e ± from 
non-heavy-flavor sources 




8.2 




e ± from J/ij} 


Bj/^j 


photonic e ± from J/^ decays 




5.3. 


L 


e from T 


By 


photonic e from T decays 




5.3.1 


> 


Drell-Yan e 


Bd-y 


photonic e from Drell-Yan process 




i.3.! 


3 



This correction procedure will be discussed further in Chapter [9j 
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Chapter 4 

e^ 1 Reconstruction Efficiency 



4.1 Introduction 

This section describes the calculation of the e ± reconstruction efficiency (sr), the 
efficiency with which an e ± track is identified as an e ± using all identification cuts 
except the TPC energy-loss cut. The efficiency is calculated with the use of simulated 
non-photonic e ± embedded into real 200-GeV Cu + Cu events. 

The efficiency Er is the pr-dependent ratio of N rec . (the number of reconstructed 
e ) to N sim . (the number of simulated e ). For a given transverse-momentum bin 
with p™ n ' < pj< < p™ ax \ N rec , is the number of reconstructed tracks that satisfy the 
following conditions: 

• The reconstructed track shares at least one half of its TPC points with a sim- 
ulated e ± track. 



Table 4.1: The combined efficiency of these e identification cuts is found using the 



method described in this chapter. These cuts are discussed in Section [372} For reasons 
discussed in the text, the TPC energy-loss cut is not included here. 



radius of first TPC point 


Ritpc <102 cm 


number of TPC fit points 


Nrpcfit > 20 


TPC fit points ratio 


R FP > 0.52 


DCA to primary vertex 


GDC A < 1.5 cm 


track pseudorapidity 


-0.1 < 7] < 0.7 


energy in BEMC 


p/E Tower < 2/c 


BSMD cluster displacement 


A 5 md < 0.02 


BSMD cluster sizes 


Nsmdv > 2 and N SMD(j} > 2 


BEMC status 


"good" for tower, BSMD-??, and BSMD-0 
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• The track passes all e identification cuts (see Table 4.1 ) except the TPC energy- 
loss cut. 

• The reconstructed transverse momentum and pseudorapidity of the track are 
within the ranges p™ m ' < pT(rec.) < p™ ax ' and —0.1 < r)(rec.) < 0.7. 

For a given pt bin, iV s j m , is the number of simulated e ± tracks for which the 
simulated transverse momentum and pseudorapidity are within the ranges p™ m ' < 
p T {MC) < p^ ax - and -0.1 < r](MC) < 0.7. Note that the tracks in N rec , are not 
strictly a subset of the tracks in N s i m _. Imperfect momentum resolution can cause 
a reconstructed track to be in a different pt bin than the corresponding simulated 
track. The factor er therefore allows the non-photonic e spectrum to be corrected 
for the STAR detector's momentum (and pseudorapidity) resolution. 

The calculation of Er presented in this chapter uses embedding data. In embed- 
ding, simulated particle tracks are generated and their interactions with the STAR 



detector are simulated using GEANT 115 , producing hits in the TPC, BEMC, and 



other detector components. These simulated hits are then mixed with real data and 
the combination of the two is run through STAR's event reconstruction (tracking, 
vertex-finding, etc.) algorithms. This allows the efficiency of the STAR detector to 
be studied in a realistic (high-particle-multiplicity) background. 

The BEMC simulator uses the calibration and status of the BEMC for the real 
event into which the simulated data are embedded. Therefore, any losses in BEMC 
acceptance in the real data will also appear in the embedding data. For this reason, 
the reconstruction efficiency is corrected by the mean BEMC acceptance (the 



calculation of this factor is described in Chapter 5.1) 



£r[Pt) = j~ A 1-. — rr • J , r , (4.1) 

{A BEMC {sim.)) N sim \p T ) 

where {ABEMc(sim.)) is the BEMC acceptance for the set of events used in embed- 
ding, which is a subset of all real Cu + Cu events. Since the BEMC acceptance 
changes with time, the value of (A BEMC ) may be different for the embedding subset. 
Three embedding data sets are used in this chapter. 

• Data Set Ni. This is the primary embedding data set used to calculate Er. 
Simulated single e ± are embedded into « 670,000 real 200-GeV Cu + Cu events. 
The e were generated uniformly in the transverse-momentum range 1 GeV/ c < 
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p T (MC) < 15GeV/c and the pseudorapidity range —0.3 < rj(MC) < 1.3, 
covering the full azimuth. The e were embedded at a rate of ~ 0.05 times the 
reference multiplicity of the underlying event. For this data set, (Abemc(Ni)) = 
0.606. This data set was created using an incorrect function in the vertex-finding 
algorithm; as described below, this error did not have an effect on the calculated 
efficiency 

Data Set N 2 : This data set is used to study the effect of a flaw in data set Ni. 
This data set is very similar to data set N±, but was generated using the correct 
function in the vertex-finding algorithm. This smaller data set included only »s 
150,000 events. For this data set, (A BEMC (N 2 )} = 0.613. 

Data Set Pi. This data set is not directly used in the calculation of Er, but some 



distributions of measured quantities from this data set are shown in Section 4.4 
Simulated photonic e (from 7r° — > 77 decays and the subsequent conversion 
of those photons to e~e + pairs in the inner structures of the STAR detector) 
are embedded into ~ 800,000 real Cu + Cu events. This data set is used to 
calculate the background rejection efficiency (see Chapter [7j which contains a 
more detailed description of this data set). 

Note that the embedding data sets used for this calculation did not include a 
satisfactory simulation of TPC energy loss; the efficiency of the TPC energy-loss cut 
is calculated separately in Chapter |6j Unless otherwise noted, the statistical uncer- 
tainties in Er shown in Figures in this chapter are estimated uncertainties calculated 



using the method described in Section B.l 
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4.2 Weighting 



In embedding data sets N\ and N2, the simulated e are distributed uniformly in 
transverse momentum over the range 1 GeV/c < pt < 15 GeV/c. If the correction fac- 
tor Er is to properly account for the STAR detector's momentum resolution, the e ± 
spectrum must be weighted to have a realistic dependence on p?- For example, con- 
sider two e tracks, one with simulated transverse momentum pt(MC) = 3GeV/c, 
the other with pt{MC) = 5GeV/c, with both tracks having the same reconstructed 
transverse momentum p^(rec) = 4GeV/c. In a realistic spectrum, the parti- 
cle yield decreases with increasing px- Thus, there should be fewer electrons with 
Pt(MC) = 5 GeV/c than with px(MC) = 3 GeV/c. If px weights are not applied, the 
two tracks will make equal contributions to the calculation of Er and the value of the 
efficiency at px = 4GeV/c will be too large (given that Er increases with p T ). 



PHENIX p+pNPE(x10) 
fit F[e N ,1]x10 

STAR p+p NPE 




p r [GeV/c] 



Figure 4.1: Non-photonic e ± cross-sections (vs. pt) for 200-GeV p + p collisions 



measured by the PHENIX 51 and STAR 1 80 1 collaborations. Statistical uncertainties 



are indicated by error bars; systematic uncertainties are indicated by the gray shaded 
boxes. Also shown are power-law fits to the data (for the fits, the statistical and 
systematic uncertainties of each data point were added in quadrature). The fit to 
the PHENIX data, F[ej^, 1] is used as the initial weighting function in the calculation 
of Er. The dashed cyan curve is F[eN, 1] without the scaling factor of 10. The 



parameters for the fit to the STAR data (F[ejy, 2]) are given in Section D.6 
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The weighting function should have the same shape as the spectrum of non- 
photonic e ± in 200-GeV Cu + Cu collisions]]] However, this spectrum is not known 
a priori, so the following iterative weighting procedure is used. An initial weighting 
function is chosen as the expected shape of the non-photonic e ± spectrum. The e re- 
construction efficiency is calculated with each e ± weighted according to its simulated 
transverse momentum. The values of Er calculated with this initial weighting func- 
tion are used to find the fully corrected non-photonic e spectrum for 200-GeV Cu 
+ Cu collisions. A fit to that spectrum is then used as a weighting function in a new 
calculation of Er. This procedure is repeated until the values of the e ± reconstruction 
efficiency stabilize. 

Figure |4.1| shows the initial weighting function used in the iterative momentum 
weighting procedure: the non-photonic e cross-section found by the PHENIX col- 
laboration [5l] in 200-GeV p + p collisions. That spectrum is fit with a power-law 
function F[cnA}(Pt) = A(l +Pt/po) u - The parameters of this fit function, as well 



as others discussed in this dissertation are given in Appendix D The e spectra used 
to calculate Er are weighted by factors of px(MC) x F[e^, 1]. 



Figure 4.2 shows the corrected non-photonic e spectrum for the — 60% most 
central 200-GeV Cu + Cu collisions after the first two iterations of the momentum 
weighting procedure. See Chapter [9] for a detailed description of the correction pro- 
cedure. The spectrum from iteration is used as the input weighting function to find 



Er for iteration 1. Figure 4.3 shows er calculated for various weighting functions, 
including three iterations of the weighting procedure. The efficiencies for iteration 1 
and iteration 2 are virtually identical, so the iteration procedure was stopped and the 
non-photonic e spectra from iteration 1 were used in subsequent calculations of the 
nuclear modification factor Raa- 



Figure 4.4 shows the effect of the STAR detector's momentum resolution on er. 



The black histogram shows Er calculated according to its usual definition given in 



Section 4.1| The green histogram shows Er calculated assuming perfect momentum 



resolution (z.e., the spectrum in the numerator of the calculation is a function of the 
simulated transverse momentum pt{MC), instead of pT{rec.)); the weighting function 
is F[e N , 1]. 



1 In data set Pi, each photonic was weighted according to the transvers e mo mentum of the 



pion that produced it. The weighting function, F[7r°,2], is described in Section 7.2 
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Non-Photonic (e'+e + )/2, 200-GeV Cu+Cu, cent. 0-60% 




p r [GeV/c] 



Figure 4.2: Two iterations of the corrected non-photonic e ± spectrum (vs. px) in 
200-GeV Cu + Cu collisions. Statistical uncertainties are indicated by error bars; 
systematic uncertainties are indicated by the boxes around each data point. Also 
shown are two of the weight functions used in the calculation of Er. The initial 
weighting function is .F[ejv, 1], the fit to the PHENIX collaboration's non-photonic 
e spectrum in p + p collisions. It was used to generate iteration of Er, which was 
used to generate iteration of the non-photonic e spectrum for Cu + Cu collisions 
(black). That spectrum was fit with a power-law function F[eN,Io], which was used 
as a weighting function in the calculation of iteration 1 of Er. This was used to 
generate iteration 1 of the non-photonic e spectrum for Cu + Cu collisions (red). 
Note that F[ew, Iq] was fit to the data using the function's integral over each px bin, 
not its value at the bin center. 
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e± Rec. Eff. (e R ) for Various p t Weighting Functions 
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Figure 4.3: The e reconstruction efficiency as a function of pt, calculated using vari- 
ous pt weighting functions. The differences between iterations 1 and 2 are negligible. 
The decreasing efficiencies at high px are due to momentum resolution and the fact 
that the simulated transverse momentum does not exceed 15 GeV/c. 
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Figure 4.4: The changes in er as a function of pt when the STAR detector's mo- 
mentum resolution is taken into account. The black histogram (which includes mo- 
mentum resolution effects) is calculated as defined in Section 4.1 To generate the 



green histogram, the numerator of the calculation was taken to be a function of the 
simulated transverse momentum (pt{MC)) rather than the reconstructed transverse 
momentum (pT(rec.)). 
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4.3 Calculated Efficiencies from Embedding 

In this section, various calculations of the e reconstruction efficiency are shown. 
Unless otherwise stated, the efficiencies are extracted from embedding data set Ni, 
with weighting function F[e^, 1]. 



Figure |4.5| shows the e reconstruction efficiency for various collision centrality 



classes. No correlation between er and collision centrality is observed. 



Figure [476] shows the values of Er extracted from embedding data sets iVi and N 2 . 
As was mentioned earlier, data set N 2 was created due to the existence of a potential 
flaw in data set Ni. In both of those embedding data sets, the simulated non-photonic 
e tracks originate from the primary vertex of the underlying event. In data set N\ 
the STAR vertex-finding algorithm was applied to all reconstructed tracks (including 
tracks in the underlying event and tracks associated with simulated e ) and a new 
primary vertex was found. The global DCA for each track was calculated relative to 
this new vertex, which would not be at the same position as the initial vertex of the 
embedded particles. If the two vertices were separated by more than the global DCA 
cut (1.5 cm), a simulated non-photonic e would not be counted in the numerator of 
Er even if its track were reconstructed perfectly. In data set N 2 , the global DCA of 
each reconstructed track was calculated relative to the original vertex of the collision 
(the vertex from which the simulated tracks also originated)]^] No difference in Er 



based on the choice of vertex finding method is observed in Figure 4.6 Since data 
set Ni contains a greater number of events, that data set is used for subsequent 
calculations of Er. 



Calculations of Er, with the values of various cuts changed (see Table 4.1), are 



shown in Figures 4.7 through 4.9, with some additional Figures presented in Sec- 



tion C.l.l In these Figures, the black circles indicate the calculation of Er with all 



cuts having their standard values. 



Figure |4.7| shows Er with various sets of cuts applied. Note that with no BSMD 
cuts applied, the efficiency has little dependence on transverse momentum. 



Figures 4.8 and 4.9 show er calculated with various values for the cuts on p(rec.) / 'Erowe 
and the minimum shower sizes in the BSMD-r? and BSMD-0 planes, respectively. 
Figure 4.10 shows the e reconstruction efficiency fit with a function of the form 



2 For a reader familiar with the STAR collaboration's nomenclature: data set Ni was generated 
using the VFMinuit vertex-finding option, while data set N2 was generated using the VFMCE option. 
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Table 4.2: Parameters of fits of the e reconstruction efficiency. 



Parameter 


Value 


a 


0.97539 ± 0.02464 


ai 


i.32f8±0.0i42 


bo 


0.37873 ± 0.00811 




(2.6216 ± 0.2498) -10" 2 


b 2 


(-9.5544 ± 1.7578) -10" 4 


X'/NDF 


57.2/16 



|{tanh [a Q [p T - ai)] + 1} (b + hp T + b 2 p 2 T ) 



(4.2) 



The fit parameters are given in Table 4.2 The dashed curves in Figure 4.10 show the 



±8% relative systematic uncertainties assumed for er (see Section 4.4). 
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Figure 4.5: The e reconstruction efficiency as a function of pt for various collision 
centrality classes; Er is not observed to change with the collision centrality of the 
underlying event. 
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Figure 4.6: The e ± reconstruction efficiency as a function of pr, calculated using two 
different vertex finders (see the text for an explanation). Using the incorrect vertex 
finder (data set Ni) does not affect the value of Er. 
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e± Rec. Eff. (e R ), MB trigger, cent. 0-60% 
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Figure 4.7: The e reconstruction efficiency as a function of p?, calculated when 
various sets of identification cuts are applied. The GDC A cut is always applied. 
For the "associated" (gray) histogram, no other cut is applied. "TPC" indicates 
the combination of the cuts on Ritpc, Ntpc/u, an d Rfp (but not TPC energy 
loss); "TOW" indicates the cuts relating to BEMC towers; "SMD" indicates the cuts 
relating to the BSMD. Note that the correction factor of (Abemc{N\)) has not been 
applied to these calculations of Sr. 
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e± Rec. Eff. (e R ), MB trigger, cent. 0-60% 
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Figure 4.8: The reconstruction efficiency as a function of pt with various upper 
and lower cuts on p(rec.) / E-rower- 
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Figure 4.9: The e reconstruction efficiency as a function 
on the minimum shower sizes N S mdt) and N S MD<j>- In the le. 
{NsMDr, > X AND N SMD(f) >Y). 



of pt for various cuts 
gend (X, Y) stands for 



91 



Fit of e 1 Rec. Eff. (e R ), 200-GeV Cu+Cu, cent. 0-60% 
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Figure 4.10: Fit of Er as a function of px with the form of equation 4^2, a polynomial 
multiplied by a sigmoid function. The dashed curves indicate ±8% relative system- 
atic uncertainties, see Section 4.4 See Table 4.2 (page 88) for the values of the fit 
parameters. 



4.4 Comparison of Embedding Data to Real Data 
and Systematic Uncertainties 

In this section, the e in the embedding data are compared to the e found in the 
real data; discrepancies between the two data sets are quantified and incorporated 
into the systematic uncertainties of Er. The effect of each e identification cut on 
Er (when all other cuts are applied) is defined to be the the partial efficiency of that 
cut. For example, the partial efficiency of the NTPCfu cut is a ratio: the numerator 
is the number of reconstructed e ± tracks that pass all e identification cuts; the 
denominator is the number of reconstructed e ± tracks that pass all identification 
cuts except the Ntpcju cut. This is equivalent to the ratio of Er (with the Ntpcju 
cut applied) to the value of Er calculated without the Nxpcfit cut. 



Figure 4.11 shows four distributions of p/Exower- for reconstructed non-photonic 



e from embedding data set N% (green), reconstructed photonic e from embedding 
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data set Pi (blue), inclusive e in real data (black), and photonic e 1 * 1 in real data (red, 
found by subtracting the combinatorial background from the number of unlike-charge 
pairs). Note that the distributions have been normalized so that their integrals are 
all unity. To produce these distributions, all other e ± identification cuts were applied 
(the TPC energy-loss cut was applied to find e in the real data only). Note that the 
distribution for inclusive in the real data does include some hadron contamination 
(although the purity is high for the pr bin shown); the distribution for photonic e 
in the real data should not be contaminated by hadrons. 



As shown in Figure 4.11, the p/Exower distributions for simulated e are slightly 
narrower and have maxima at higher values oip/ Exower than the distributions for real 
e^^\ For each transverse-momentum bin, distributions like those shown in Figure 4.11 



can be used to calculate the partial efficiency of the p/E Tower < 2/c cut; the results 



of these calculations are shown in Figure 4.12 If the embedding data sets precisely 
reproduce the interactions of e with the STAR detector, there should be no difference 
between the partial efficiencies calculated for real data and embedding. For the 
purposes of this section, it is most relevant to compare the partial efficiency for 
photonic e in the real data (red), which are very pure, to the partial efficiency for 
e in data set N% (green). The differences between the real and simulated p/Exower 



distributions (Figure 4.11) for minimum-bias events result in noticeable (though not 
large) differences in the partial efficiencies]^] These discrepancies will be incorporated 
into the systematic uncertainty of Er (explained below). 



3 These observations apply only for minimum-bias triggered e . 

4 For high-tower triggered events, the partial efficiencies are essentially 1 until pt becomes greater 
than approximately twice the trigger threshold (the value of the threshold is approximately 3.75 
GeV). 
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The partial efficiencies have been calculated for each of the e identification cuts. 
Some of these calculations are presented in this section, while additional material is 



presented in Section |C.1.2[ Figures |4.11| through |4.20| show example distributions 
for each quantity (in one pt bin) and the calculated partial efficiencies for each cut. 
Partial efficiencies for minimum-bias e ± are shown for < 4 GeV/ c, while partial 
efficiencies for h igh-tower triggered e are shown for pt > 4GeV/c. 
shows the distributions of A$md = 



Figure 



4.13 



(At] S md) + (A0 5 a/d) ; the 
distributions for simulated e have maxima at lower values of A$md than the dis- 
tributions for real e . Figure 4.14 shows the partial efficiencies of the A$md < 0.02 
cut; the noticeable differences between the real and simulated Asmd distributions do 
not result in large differences in the partial efficiencies. 



Figure 4.15 shows the distributions of Nsmdt)] the distributions for simulated e 
have maxima at lower values of NsMD-q than the distributions for real e ± . These 



differences are present for all transverse-momentum bins. Figure 4.16 shows the 
partial efficiencies of the Nsmd v > 2 cut; the obvious differences in the Nsmdti 
distributions do not lead to large differences in the partial efficiencies (this would not 
be the case for a more restrictive cut, e.g. N SMDrj > 3). 



Figure 4.17 shows the distributions of Nsmd</>', the differences among the Nsmd<(> 
distributions are less pronounced than the differences among the Nsmdtj distributions. 
Figure 4.18 shows the partial efficiencies of the N S MD<f> > 2 cut. 



Figure 4.19 shows the distributions of GDC A; there are visible differences between 
the GDC A distributions for real and simulated photonic e ± . Figure 4.20 shows the 
partial efficiencies of the GDC A < 1.5 cm cut. At low transverse momenta, the 
differences in the GDC A distribution between non-photonic and photonic e produce 
a small difference in the partial efficiency of the GDC A cut. Note that photonic e 
would be expected to have a wider GDC A distribution than non-photonic e . 

For each cut, the discrepancy between the real and embedding data can be quan- 
tified by finding the ratio of partial efficiencies rx(pr)/ex(pr), where rx is the partial 
efficiency for cut X from the real data (photonic e ) and ex is the partial efficiency 
for cut X from embedding data set N\. If the e identification cuts were all indepen- 
dent of each other (i.e., applying or removing cut X has no effect on the efficiency of 
cut Y), the true value of the e ± reconstruction efficiency could be found by correcting 
the value of Br found from embedding: 
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£R,true(PT) = £r{Pt) J [ 



X 



rxjpi 
exiP' 



•T) 



e R (p T )D(p T ), 



(4.3) 



where D denotes the product for all cuts of the ratios of the partial efficiencies. 
In practice, the cuts are probably not all independent of each other and the partial 



efficiencies for the real data exhibit large statistical fluctuations, so Equation 4.3 would 
not be expected to give a reliable value for the true e reconstruction efficiency. Also, 
in some cases (such as the GDC A distribution) differences between the photonic and 
non-photonic e ± would be expected, and the use of photonic e ± to represent all real 
would be misleading. Nevertheless, the product D can be used as an estimated upper 
limit on the deviation of the true e ± reconstruction efficiency from the calculated Sr. 



Figure 4.21 shows the product D, which is generally within a few percent of unity. 
A relative systematic uncertainty of ±8% is assumed for the calculated values of Sr 



(see also Figure 4.10). 
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MB Trigger, p/E Distributions, cent. 0-60%, 2.0<p <3.0 
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Figure 4.11: Distributions of p/E Tower for real and simulated e ± . 
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Figure 4.12: Partial efficiencies of the cut p/E Tower < 2/c. 
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MB Trigger,A SMD Distributions, cent. 0-60%, 2.0<p T <3.0 



en 
+^ 

g 0.25 

o 

~o 

0) 
N 

=5 0.2 
E 



0.15 



0.1 



0.05 




* real inclusive e 1 

* real pho. e 1 

* sim. single e 

O sim. pho. e* (from jt° 



o.c 



0^4 



,rffo .L. 



0708 



Figure 4.13: Distributions of A$md for real and simulated e 
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Figure 4.14: Partial efficiencies of the cut Asmd < 0.02. 
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MB Trigger, N Distributions, cent. 0-60%, 2.0<p <3.0 
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Figure 4.15: Distributions of Nsmdt) for real and simulated e d 
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Figure 4.16: Partial efficiencies of the cut Nsmd^ > 2. 
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MB Trigger, N Distributions, cent. 0-60%, 2.0<p <3.0 
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Figure 4.17: Distributions of Nsmd^ for real and simulated e 
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MB and HT triggers, Effect of N smd ^ Cut on e R , cent. 0-60% 
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Figure 4.18: Partial efficiencies of the cut Nsmd4> > 2. 
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MB Trigger, GDCA Distributions, cent. 0-60%, 2.0<p T <3.0 
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Figure 4.19: Distributions of the global DCA for real and simulated e 
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MB and HT triggers, Effect of GDCA Cut on e R , cent. 0-60% 
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Figure 4.20: Partial efficiencies of the cut GDCA < 1.5 cm. 
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Product (D) of Ratios of Partial Efficiencies, MB and HT triggers, cent. 0-60% 
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Figure 4.21: The product (D) of the ratios of partial efficiencies for all cuts (D is 



defined in equation 4.3). For reasons discussed in the text, \D — 1 is treated only as 



an approximate upper limit of deviations of the calculated value of Er from its true 
value. With the exception of one pt bin, D lies within the range 1 ± 0.08, therefore 
±8% relative systematic uncertainties are assumed for sr. 
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Chapter 5 



BEMC Geometrical Acceptance 
and High- Tower Trigger Efficiency 

5.1 BEMC Geometrical Acceptance 

This section describes the calculation of correction factors to account for time-dependent 
changes in the geometrical acceptance of the STAR Barrel Electromagnetic Calorime- 
ter (BEMC). The data presented in this dissertation were recorded over the course of 
33 days. During that time the geometrical acceptance frequently changed as problems 
were encountered and the various towers and BSMD strips were turned on and off. 
It is necessary to correct the e yields to account for these changes in acceptance. 

The BEMC geometrical acceptance A BEMC is calculated for each run (a data- 
recording period usually lasting from a few minutes to one half hour, during which 
thousands of collisions are recorded). For each run, the status of every tower and 
every BSMD strip is recorded. The following method is used to determine the status 
of BSMD strips (a similar method was used to determine the status of BEMC towers). 
For every event, the signal in each strip is read out by an Analog-to-Digital Converter 
(ADC). The rate at which each strip records a signal greater than five standard 
deviations above the pedestaQ value is recorded and the mean rate rate for all strips 
is calculated. Strips that deviate from this mean rate by less than a factor of two are 
assigned "good" status, while strips that deviate above or below the mean rate by 
more than a factor of two are assigned "bad" status. 

lr The pedestal value is the mean ADC signal when no shower is present in the detector. 
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To calculate A BEMC , 10 6 points are generated uniformly distributed in pseudo- 
rapidity and azimuth in the ranges < rj < 0.7 and < (f) < 2n. The BEMC 
geometrical acceptance is the fraction of those points for which the corresponding 
tower and both BSMD strips all have good status. An example of this procedure is 



shown in figure 5.1 The mean BEMC geometrical acceptance {A BEM c) is calculated 
by finding the weighted average of Abemc for all runs, with the weight for each run 
being the number of events in that run that passed the event selection cuts used 
in this analysis. For minimum-bias (high-tower) triggered runs, the mean BEMC 
geometrical acceptance is 0.634 (0.642). 

To avoid possible calibration problems, towers/strips with bad status were ex- 
cluded from the analysis of the embedding data. To produce the embedding data 
sets, simulated events were embedded into a subset of sixteen runs of real events. 
The mean acceptance (A BEMC ) was calculated for two embedding data sets (see 
Chapter [4]) in the manner described above, with the weight for each run being the 
number of simulated events that passed the event selection cuts. For data sets Ni and 
N 2 the values of the mean BEMC geometrical acceptance are (A BEMC (Ni)) = 0.606 
and {A BEM c{N 2 )) = 0.613. These values will be used to correct the calculated values 
of the e reconstruction efficiency (er) for the BEMC geometrical acceptance in the 
embedding data. 
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BEMC status, run 6041118 
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Figure 5.1: The status of the BEMC for one run. "TOW" indicates that the tower at 
the given location has good status; "SMDE" ( "SMDP" ) indicates that the BSMD-77 
(BSMD-0) strip at the given location has good status. 
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5.2 High- Tower Trigger Efficiency 



As described in Section 3.1, the high-tower (HT) trigger is used to select events in 
which at least one BEMC tower has energy > 3.75 GeV The HT triggered data set 
therefore contains an enhancement of high-momentum e . The trigger enhancement 
is corrected by dividing the e yield by N^J (normalized) , the number of minimum- 
bias (MB) triggered events that were observed during the time the HT triggered 



events were recorded (see Section 3.1 ). However, it is still necessary to correct for the 
HT trigger efficiency (et), which is momentum-dependent. 

Figure 5.2 illustrates the calculation of Et- The efficiency-corrected inclusive e 



yields for the MB (black) and HT triggered (red) data sets are shown (for the 0- 
60% centrality class). For the HT triggered spectrum, all correction factors except 
Et have been applied. A combined spectrum, consisting of the MB spectrum for 
Pt < 5.5GeV/c and the HT spectrum otherwise, is fit with a function of the form 
A{1 + pt/B) u + exp[— (pt — C)/D]. The HT spectrum is divided by this fit (lower 
panel) and fit with a function of the forr^ 



Et = A — exp sinh ^ T — . (5.1) 

The HT trigger efficiency is assigned an absolute ±2.5% systematic uncertainty, which 
covers the statistical scatter of the scaled HT histogram. 

The method described above permits the shape of the trigger efficiency to be 
determined, but does not determine the overall magnitude of Et- The high-p^ value 
of Et is estimated using embedding (data set Ni, see Chapter [5]). The simulated Et 
is the (pr-dependent) ratio N rec ^ trig /N rec . Here, N rec . is the number of simulated e ± 
tracks that are reconstructed and identified as e ± . N rec . jt rig. is the subset of those e ± 



that fulfill the HT trigger condition. The simulated Et is shown (circles) in Figure 5.3 



The HT trigger efficiency calculated using real data (stars in Figure 5.3) is scaled so 
that it matches the simulated efficiency at high p?. 



2 This function is not a sigmoid and does not describe the scaled HT spectrum at low pt- It does, 
however, describe the scaled spectrum for p T > 4GeV/c. 
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Calculation of Trigger Efficiency (e T ) 
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Figure 5.2: Illustration of the calculation of Et, the high-tower trigger efficiency. The 
upper plot shows the inclusive e spectra for the minimum-bias (black circles) and 
high-tower triggered (red stars) data sets. All correction factors have been applied 
to the HT spectrum except Et- A combined spectrum (MB for pt < 5.5GeV/c and 
HT for p T > 5.5GeV/c) is fit with a fit function (blue curve). The lower plot shows 
the ratio of the HT spectrum to the spectrum fit. The trigger efficiency is fit with 



a function of the form of Equation f^l; the systematic uncertainties assigned to the 
trigger efficiency are shown as dashed lines. See the text for further explanation. 
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HT Trigger Efficiency (e T ) 




p r [GeV/c] 



Figure 5.3: The HT trigger efficiency (et) vs. transverse momentum calculated using 
real data (stars) and embedding data (circles). The efficiency from real data (and the 
fit thereof) is scaled to match the embedding data at high pt- 
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Chapter 6 

Energy-Loss Cut Efficiency and 
Purity 



6.1 Introduction 

The energy-loss cut efficiency (e dE/dx) is calculated by fitting the normalized energy- 
loss (na e ) distribution for photonic e ± with a Gaussian. The inclusive e ± purity 
(K inc ) is calculated using the normalized energy- loss distribution of particles that 
pass all other e identification cuts: that distribution is fit with multiple Gaussians 
to account for the contributions of e and various hadron species. These calculations 
will be described in greater detail below. 



As described in Section 2.3.2, na e} the normalized energy loss (relative to the 
expected value for e*) is defined as 

= 1. (6.1) 



Here, (dE/dx) is the measured 70% truncated mean value of the TPC energy loss, B e 
is the most probable value of the 70% truncated mean, and a e is the expected width of 
the In ((dE/dx) / B e ) distribution. The factors B e and a e are both derived from Bich- 
sel's simulations [108] of particle energy loss in the TPC gas. Figures 6.1 and 6.2 show 



na e versus momentum for minimum-bias and high-tower triggered triggered events, 
respectively. The curves indicate the approximate locations of the most probable 
value of na e as a function of momentum for five particle species (see Appendix [A] for 
more information on the calculation of these curves). Also shown (horizontal lines) 
are the values of the cut (—1.5 < na e < 3.5) used to identify e ± . 
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200-GeV Cu+Cu, MB trigger, na, for charged particles, p/E and SMD cuts, cent. 0-60% 
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Figure 6.1: Normalized TPC energy loss (na e ) vs. particle momentum for MB trig- 
gered events. The tracks shown here have passed all e identification cuts except the 
energy-loss cut. The curves show the estimated most probable value of na e for each 
of five particle species. The horizontal lines show the range of the na e cut used in 
this analysis. 



200-GeV Cu+Cu. HT trigger, no e for charged particles, p/E and SMD cuts, cent. 0-60% 
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Figure 6.2: Normalized TPC energy loss (na e ) vs. particle momentum for HT trig- 
gered events. The tracks shown here have passed all e ± identification cuts except the 
energy- loss cut. 
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6.2 Energy-Loss Cut Efficiency 



The efficiency of the energy-loss cut is calculated by fitting the na e distribution for 
photonic e with a Gaussian. Figure 6.3 illustrates the calculation of the na e distri- 
bution for photonic e . In black is the na e distribution for pairs of oppositely charged 
tracks; in blue is the combinatorial background calculated from the na e distributions 
for like-charge tracks. The red histogram shows the difference between the other two: 
the distribution of na e for photonic e^ 1 ^ Note that the contribution from hadrons 
(tracks with negative na e ) is largely cancelled out. This is the reason for using pho- 
tonic e ± to determine SdE/dx'- a fit of the e ± peak is less likely to be biased in the 
negative direction due to hadron contamination. 



Figure 6.4 shows the na e distribution for photonic e in minimum-bias collisions 
(0-60% centrality class). For each centrality class and trigger type, the e ± peak is 
fit with a Gaussiar0 with mean fi and width 5. These fit parameters, which are 
given in Table [fTTj deviate significantly from their expected values (/i = and 5 = 1, 



respectively), indicating differences in the TPC calibration between the real data 
and the simulations used to determine the normalization factors B P and a P . The fit 



parameters for different momentum bins are shown in Figures 6.5 (for minimum-bias 
triggered events) and 6.6 (for high-tower triggered events) for the 0-60% centrality 
class. No momentum dependence is observed]^] 

The energy-loss cut efficiency is the ratio of the area of the Gaussian fit within 
the range —1.5 < na e < 3.5 to the total area of the Gaussian. For each trigger 
type, the fit for a single large momentum bin (2 GeV/c < p < 6.5 GeV/c for minimum 
bias, 4GeV/c < p < 11.3GeV/c for high tower) is used to calculate SdE/dx- The 
uncertainties in the efficiency are calculated by changing the fit means and widths 
from their central values until the fit x 2 has increased by one and again calculating 
the efficiency. The uncertainties in the efficiency due to changing the fit mean and 
changing the fit width are added in quadrature to give the quoted uncertainty. 



3.3 



^or more information on the background-subtraction method, see Section 
2 of the form Aexp{— ^[(x — /i)/<5] 2 }, where A is the amplitude 

3 For minimum-bias triggered events, the na e distributions for 3.5 GeV/c < p < 4 GeV/c have low 
numbers of entries, leading to fluctuations in the fit parameters. 
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Table 6.1: Fit parameters for the photonic e peak and energy-loss cut efficiencies 
for various centrality classes and event triggers. 





Mean (//) 


Width (S) 


e dE/dx 


MB, 0-20% 


-0.41932 ± 0.03007 


0.87788 ± 0.02519 


0-890841™ 


MB, 20-60% 


-0.33308 ± 0.02638 


0.86441 ± 0.02268 


0.91148+^™ 


MB, 0-60% 


-0.37130 ± 0.02070 


0.86616 ±0.01769 


u.yuo/o_ 01444 


HT, 0-20% 


-0.46847 ± 0.01945 


0.92948 ±0.01587 


8664S +U - Uiiti4 
u.ooo^o_ .01160 


HT, 20-60% 


-0.41376 ± 0.02059 


0.91234 ±0.01689 


0.88309j#K 


HT, 0-60% 


-0.43650 ± 0.01435 


0.91620 ±0.01174 


n 8771 O+0.00882 
U.OI 1 l^_o.00879 



200-GeV Cu+Cu, MB trigger, na e distributions for pairs, 0-60%, 2.0<p<2.5 
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Figure 6.3: The calculation of the na e distribution for photonic e ± . The abscissa is 
the value of na e for the e candidate (i.e., the member of the pair that has passed all 
other identification cuts). Shown are the distribution of na e for unlike-charge pairs 
(black), the combinatorial background computed from the na e distributions for like- 



charge pairs (see Section 3.3), and the na e distribution for photonic e ± (the difference 
between n un iike an d the combinatorial background). 
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200-GeV Cu+Cu, MB trigger, na e for photonic e , cent. 0-20%, 2.0<p<6.5 
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Figure 6.4: The na e distribution for photonic e ± is fit with a Gaussian. The cuts on 
na e used for this analysis are shown by the vertical lines at na e = —1.5 and na e = 3.5; 
the area of the Gaussian fit within those limits is the energy- loss cut efficiency SdE/dx- 
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200-GeV Cu+Cu, MB trigger, photonic nc e fit mean and width, cent. 0-60% | 
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Figure 6.5: Means (/i) and widths (5) of fits to the na e distribution for photonic e ± 
as a function of momentum for minimum-bias triggered events. The vertical error 
bars indicate the uncertainties in the fit means. The boxes indicate the widths of the 
fits {n ±8). 
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Figure 6.6: Parameters of fits to the na e distribution for photonic e ± in HT triggered 
events. 
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6.3 Purity 



The inclusive e purity (Ki nc ) is calculated by fitting with three Gaussians the na e 
distribution for all tracks that pass the other e identification cuts, as illustrated in 



Figure 6.7 The Gaussian centered near na e = accounts for the contribution from 
and its mean and width are constrained to the values found from the (very pure) 



photonic e ± na e distributions described in Section 62 The other two Gaussians 
account for the contributions from hadrons, expected to be predominantly 7r ± , K ± , 
p, and p. The middle (blue) peak is intended to account primarily for charged pions 
(although it is not possible to guarantee that this peak accounts for all pions and only 
pions). Such three-Gaussian fits are performed for multiple momentum bins. When 
possible, the hadron peak means and widths are allowed to vary freely. However, it is 
occasionally necessary to constrain hadron fit parameters to avoid nonsensical results 
or deviations from the expected behavior of the fit (e.g., the mean of the pion peak 



suddenly decreasing as p increases, or a sudden change in a peak width). Figures 6.8 



and 6.9 show the parameters of these fits for all momentum bins for the minimum-bias 
and high-tower-triggered data sets, along with approximate expected peak locations 
for five particle species. The uncertainties of the hadron fit parameters are found by 
varying the value of each parameter until the fit \ 2 value increases by one. 



Figure 6.1 indicates that the (anti)deuteron band crosses through the range —1.5 < 
na e < 3.5 above p = 2GeV/c, resulting in contamination of the e that is not 



accounted for in the three-Gaussian fit. Figure |6.10| shows a four-Gaussian fit to the 
na e distribution for the momentum range 1.3GeV/c < p < 1.7GeV/c, where there 
is still a decent separation between the and (anti)deuteron peaks. The d(d) peak 
is allowed to have a large width to account for the fact that the mean value of na e 
for deuterons is not constant and changes significantly across that momentum bin. 
The ratio of the areas of those two peaks is 0.011. This ratio is used to estimate the 
(anti)deuteron yield in the momentum range 2 GeV/c < p < 2.5 GeV/c. 
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Figure 6.7: The na e distribution for tracks that have passed all other e identification 
cuts is fit with three Gaussians, to account for the contributions of e ± (green peak), 
tt^~ (blue), and other hadron species (red). 
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For a given momentum range, the na e distributions for the 0-20% and 20-60% 
centrality classes are not fit independently of each other. Instead, the fit for the 
full 0-60% centrality class is used to constrain the fits for the 0-20% and 20-60% 
centrality classes. As an example, consider a single momentum bin for the minimum- 
bias triggered data set. The na e distribution in the full 0-60% centrality class is fit 
with the function /, which consists of three Gaussians with amplitudes Aj (where 
j = 1,2,3), means fij, and widths Sj. The na e distribution in the 0-20% centrality 
class is fit with the function with parameters A'j, fij, and Sj. The values of some 
of the parameters of /' are constrained by the values of the parameters of /. Each 
amplitude Aj is allowed to vary freely. Each width Sj is constrained to lie within the 
range Sj — A(Sj) < Sj < Sj + A(Sj), where A(5j) indicates the uncertainty of Sj. Each 
mean fi'j is constrained to be within the range 

Hi - A(nj) - \s\ < fi'j < ^ + A(fij), (s < 0) (6.2) 

or 



fij - A(fLj) < fij < fij + A(fij) + s, (s> 0). (6.3) 

The variable s is the shift in the electron peak mean between the 0-60% centrality class 
and the 0-20% centrality class, which is determined by fitting the na e distribution for 



photonic e (see Section 6.2). The same procedure is also used to constrain the fit 
parameters for the 20-60% centrality class. 

The inclusive e ± purity, shown in Figure 6.11, is the ratio of the area of the e 
peak within the range —1.5 < na e < 3.5 to the total area of the three peaks within 
that range (the purity in the range 2GeV/c < p < 2.5GeV/c is reduced to account 
for the presence of (anti)deuterons). To find the lower uncertainty in K inc , each 
fit parameter is changed from its central value by its uncertainty in a manner that 
decreases the purity (the hadron peak means, widths, and amplitudes and the e 
peak width are increased by their own uncertainties, while the e peak amplitude 
and mean are decreased). The changes in K inc . caused by changing each parameter 
are then added in quadrature to obtain the total uncertainty. The upper uncertainty 
in Ki nc , is calculated in a similar fashion by changing the values of the fit parameters 
in the opposite directions. In the momentum range 2GeV/c < p < 2.5GeV/c, the 
lower uncertainty due to the (anti)deuteron contribution is found by calculating the 
purity with the d(d) yield doubled; the upper uncertainty is found by removing the 
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d(d) contribution. These d(d) uncertainties are added in quadrature with the other 
uncertainties. 

The calculations of the purity described above have been divided into bins in 
momentum, but the non-photonic e ± spectrum is presented in transverse momentum 
bins. It is useful (but not necessary) to convert the momentum-dependent purity 



calculation described above to a pr-dependent measurement (see Figure 6.12). This 



is accomplished by taking a weighted average. For each transverse momentum bin, 

Ys K inc.{p)I{PiPT) 

K '- M = ' EW) ' (6 - 4) 

V 

where I(p,Pt) is the measured (uncorrected) inclusive e yield in each bin in momen- 
tum and transverse momentum. 



Additional material regarding the calculation of K inc . can be found in Section C.2 



That section includes plots similar to Figure [677] showing the results of multi-Gaussian 
fits to the na e distributions for various momentum bins. That section also includes 
tables with the values of K inc {p T ) for all centrality classes and trigger types. 
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Figure 6.8: Means (/i) and widths (S) of the three Gaussians used to fit the na e 
distributions for particles that pass the other e identification cuts for minimum-bias 
triggered events. The vertical error bars indicate the uncertainties in the fit means. 
The boxes indicate the widths of the fits (// ± 5). 
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Figure 6.9: Means and widths of the three Gaussians used to fit the na e distributions 
for particles that pass the other e identification cuts for high-tower triggered events. 
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200-GeV Cu+Cu, MB trigger, Four-Gaussian Fit to na e Distribution, 0-60%, 1.3<p<1.7 
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Figure 6.10: In a procedure similar to that illustrated in Figure |6.7[ the na e distribu- 
tion in the momentum bin 1.3 GeV/c < p < 1.7 GeV/c is fit with four Gaussians, with 
the fourth (magenta) peak used to account for the contributions of (anti)deuterons. 
The ratio of the areas of the e ± and d(d) peaks is used to estimate the d(d) contami- 
nation in the e sample for the momentum bin 2 GeV/c < p < 2.5 GeV/c. 
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Purity (K. nc ) vs. p, cent. 0-60% 
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Figure 6.11: The inclusive e ± purity (i^ nc .) as a function of momentum for minimum- 
bias (black circles) and high-tower triggered (red stars) events. 
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Figure 6.12: The inclusive e purity (Ki nc ) as a function of px for minimum-bias 
(black circles) and high-tower triggered (red stars) events. 
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Chapter 7 

Background Rejection Efficiency 



7.1 Introduction 

The background rejection efficiency (eb) is the efficiency with which photonic e ± are 



identified, using the method described in Section 3.3 This efficiency is calculated from 
simulations of photon conversions and Dalitz decays of neutral pions (tt° — > 7e~e + ). 
The efficiency is found by considering the set of all reconstructed tracks that are 
identified as e and that are associated with simulated photonic e tracks. The 
fraction of those e that is identified as photonic e is Eb- 

More precisely, reconstructed are defined to be those that pass the following 
cuts. 



• The reconstructed track shares at least one half of its TPC points with a sim- 
ulated photonic e track. 

• GDCA(rec.) < 1.5 cm. 

• -0.1 < rj(rec.) < 0.7 

• Ritpc < 102 cm 

• N T p Cfl t > 20 

• R FP > 0.52 

• p(rec.) I E Tower < 2c' 1 

• Asmd < 0.02 
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• NsMDrj, NsMD(j> > 2 

• The BEMC has good status at the point to which the reconstructed track 
projects. 

Rejected e are defined to be the subset of reconstructed e ± for which the following 
conditions are also satisfied. 

• The true conversion (or Dalitz-decay) partner of the simulated e is also recon- 
structed (i.e., a reconstructed track shares at least one half of its TPC points 
with the simulated photonic e ± partner). 

• The reconstructed true partner track has RFp{partner) > 0.52. 

• The reconstructed true partner track has pxipartner) > 300 MeV/c. 

• The reconstructed tracks of the photonic e~e + pair have opposite charges. 

• The reconstructed tracks of the photonic e~e + pair have a distance of closest 
approach of DC ' A[eT e + ) < 1.5 cm. 

• The reconstructed tracks of the photonic e~e + pair have an invariant mass of 
M inv Xe~e + ) < 150 MeV/c 2 at their point of closest approach. 

The efficiency is the ratio of the number of rejected e to the number of recon- 
structed e . In the analysis of real data, a cut on the normalized TPC energy loss 
of —1.5 < na e < 3.5 is used to identify e ± , while a cut of (dE/dx) > 2.8 keV/cm is 
used to increase the fraction of e ± in the set of partner tracks. Due to difficulties in 
simulating TPC energy loss, these cuts are not applied in the analysis of simulation 
data. The energy- loss cut efficiency (e dE/dx) is used to correct for the first cut; the 
second cut is assumed to be 100% efficient (See Chapter [6]). 

Four separate simulation data sets are used in the calculation of the background 
rejection efficiency. 

• Data Set Pi: This data set is used to calculate Eb for e from photon conver- 
sions, with the photons originating from 7r° decays. Simulated 7r° — > 77 decays 
and the subsequent conversion of those photons to e~e + pairs in the inner struc- 
tures of the STAR detector (simulated using GEANT (TT5] ) were embedded into 



800,000 real Cu + Cu events. The simulated tt° were distributed uniformly 
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in transverse momentum from 0.5 GeV/c < Pri 71 ' ) <15 GeV/c and uniformly 
in pseudorapidity from —0.5 < 77(71"°) < 1.5. The pions were embedded at a rate 
of ~ 0.05 times the reference multiplicity of the underlying event. Neutral-pion 
Dalitz decays were also embedded into this data set at the natural branch- 
ing ratio (1.174%) [17], but these decays were simulated incorrectly and their 
products are not considered in this analysis. 

Data Set P<i- This data set is used to calculate Eb for e from 7T° Dalitz decays 
e + 7). Neutral-pion Dalitz decays (with the correct |116| invariant- 



mass distribution and form factor) were simulated in 10 6 p + p events gener- 



ated using PYTHIA. 75 The interactions of the decay products with the STAR 
detector were simulated using GEANT. The simulated pions were distributed 
uniformly in transverse momentum for Pt{^°) < 20 GeV/c and uniformly in 
pseudorapidity for \rj\ < 1.5. Calculations of Eb were done for the e~e + pairs 
produced directly in these Dalitz decays as well as for the e~e + produced by 
the conversion of the Dalitz-decay photons (the former contribution is larger). 
Two-photon decays of neutral pions were also simulated in this data set, but 
are not included in the final calculation of the background rejection efficiency 

• Data Set P3: This data set is used to calculate Eb for e from photon con- 
versions, with the photons from sources other than ir° decays. Photons were 
simulated in rs 10 6 p + p events generated using PYTHIA for the same pt and 
T) ranges used for the simulated ir° in data set P^. The interactions of the decay 
products with the STAR detector were simulated using GEANT. The conver- 
sions of these photons were used to calculate the effect on eb of photons that are 
not 7r° decay products (specifically, direct photons and photons from 77 — > 77 
decays) . 

• Data Set P4: This data set was used to determine the expected shape of the 
spectrum of photons from ?7-meson decays. Approximately 5 x 10 5 two-photon 
decays of rj mesons were simulated using PYTHIA. Unlike data sets P±, P2, and 
P3, the particles in this data set were not embedded into a real or simulated 
event. The interactions of the decay photons with the STAR detector were not 
simulated, so photonic e ± and the background rejection efficiency could not be 
studied directly for this data set. 
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7.2 Weighting 




Figure 7.1: STAR |m|[ll7][ll8l and PHENIX measurements |119}|i2l] of neutral and 
charged pion spectra: Plotted on the vertical axis is Ed 3 a/dp 3 forp + p spectra and 
Ed 3 N/dp 3 for the Cu + Cu spectra. Fits to the spectra are also shown; the dashed 
curves are the unsealed functions. 



In the three simulation data sets described above, the photonic e are produced 
from simulated particles (tt° or 7) that are uniformly distributed in transverse mo- 
mentum. This is, of course, not a realistic distribution for pions and photons; the 
photonic e spectra used to calculate eb must be weighted to correct for this. The 
weighting procedure is described below. 

In the calculation of the background rejection efficiency, each photonic e is as- 
signed a weight according to the pt of the ir° (for data sets P\ and P2) or photon (for 
data set P3) from which it originates. The choice of weighting function can have a 
large effect on the calculated value of the background rejection efficiency Consider 
two different photonic electrons, each with Pt{g~) = 2 GeV/c. One of these came 
from a photon with Pt{i) = 2.2 GeV/c and the other came from a photon with 
Pt(t) = 10 GeV/c. The conversion partner of the first electron must have low trans- 
verse momentum, pT{partner) m 0.2 GeV/c, and will likely fail to be reconstructed or 
fail the 0.3-GeV/c cut on p^partnev) . The first electron is unlikely to be identified 
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as photonic. The conversion partner of the second electron will have high transverse 
momentum, pT(partner) w 8 GeV/c; the second electron has a greater chance of 
being identified as photonic. A photon spectrum that decreases with pr{l) will give 
the first electron greater weight than the second, resulting in a lower efficiency than 
the un-weighted case. 

For data sets Pi and P2, each is assigned a weighting factor equal to the 



expected 7T° cross-section multiplied by pt(tt )- Figure 7.1 shows one STAR 117 and 
two PHENIX 119 120 measurements of the ir° cross-section in p + p collisions. Also 



shown are the PHENIX collaboration's measurement 121 of the n yield in Cu + Cu 



collisions and low- and high-pr STAR measurements 111 118 of the charged-pion 
((•7T - + 7r + )/2) yield in Cu + Cu collisions. Five fit functions are used to describe 
these cross-section measurements. The fit parameters of this function, as well as all 
other weighting functions discussed in this section, are given in Appendix [Dj 

• Function F[k , 1]: The PHENIX p+p — > tt° cross-section published in 2003 was 



fit [119 using a power-law function with the form A(l + pr/po 



Function F[vr°, 2]: The PHENIX p+p ->• tt° cross-section published in 2007 [120 
was fit using a function with the form A(l + p\j ' B) n + C exp(—pT/D). 



Function F[7r°,3]: The STAR p + p — > 7r° cross-section was fit |l!7| using a 
power-law function with the form A(l + Pt/pq)™- 



Function F[vr°, 4]: The PHENIX Cu + Cu ^ tt° yield (0-10% centrality class) [m 
was fit using a power-law function with the form A(l + pr/po) n - The data were 
scaled by 1/((T AA )R AA ) = 42mb/(182.7 x 0.515) before fitting. 

Function F[vr ± , 1]: The combined STAR Cu + Cu ^ tt ± yield (0-60% centrality 
was fit using a function with the form A{l+p\/ ' B) n + Ce~ PT ^ D '. 



Ill 
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class) 

The data were scaled by 1/((T AA )R AA ) = 42mb/(80.41 x 0.7) before fitting. 

The p + p neutral-pion weighting functions appear to match the shape of the Cu 
+ Cu pion spectra. Function F[7r°,2] is used as the primary weighting function for 
pions because the measurement from which it is derived has the smallest uncertainties 
and spans the largest transverse-momentum range of the measurements shown in 



Figure 7.1 The effect of the different pion weighting functions upon Eb is explored 
in Sections E3] and EU 
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RHIC p+pri-meson spectra 
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Figure 7.2: STAR 117 and PHENIX 32 measurements of ?7-meson spectra in p + p 



collisions. Fits to the spectra are also shown. 
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Figure 7.3: Generation of a simulated photon spectrum from fit F[rj, 2]: The simulated 
photon spectrum is scaled by the 39.31% branching ratio for the 77—7-77 decay. The 
simulated rj mesons are distributed uniformly in the transverse-momentum range 
lGeV/c < Pt(v) < 20GeV/c, which causes the photon spectrum to behave non- 
physically below Pt(j) = 1 GeV/c and near Pt(i) = 20 GeV/c. 



For data set P3, several different weighting functions are added together to describe 



the shape of the spectrum of //-meson-decay photons and direct photons. Figure 7.2 
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shows STAR |117| and PHENIX 32 measurements of the 77-meson cross-section in 
p + p collisions (PHENIX measured this cross section using two different 77-meson 
decay modes). Three different functions are used to describe the shape of the 77- 
meson spectrum. 

• Function F[r], 1]: The PHENIX cross-sections (from both measurements) are fit 
with a function of the form A(l +PT/po) n - 

• Function F[t],2]: The shape of the 77-meson spectrum is estimated from the 
PHENIX p + p — > 7T° spectrum by making the substitution 

Pt(tt ) — > ^/pt(v) 2 + rn(r]) 2 c 2 — m(7r°) 2 c 2 in function F[7T°,2]. The resulting 
function is multiplied by the measured [32] 77/71"° ratio of 0.48. 



• Function F[rj, 3]: The STAR cross section is fit 117 with a function of the form 
A(l+p T /p ) n . 

Function F[r), 2] is used as the primary weighting function for 77-mesons. That 
function is derived from the PHENIX tt° cross-section measurement in p+p collisions, 
which has smaller uncertainties and covers a larger transverse-momentum range than 
either of the //-meson cross-section measurements shown. The spectrum of photons 
expected from the decays of rj mesons is estimated by applying all three of these 
weighting functions to the simulated photons in data set P 4 , a set of simulated 77—7-77 
decays. The resulting photon spectra were fit with functions of the form A(l + 



Pt/po)" x Pt- This is illustrated by Figure 7.3, in which F[r], 2] is used as the rj- 
meson weighting function. The simulated photon spectrum is scaled by the 39.31% 
branching ratio [TtJ for the 77 — > 77 decay. 

• Function F[y v , 1] is the fit to the photon spectrum generated using F[rj, 1] as a 
weighting function. 

• Function F[y v , 2] is the fit to the photon spectrum generated using F[rj, 2] as a 
weighting function. This function is used as the primary weighting function for 
77-meson-decay photons because it is derived from weighting function F[7r°,2], 
the primary weighting function used for neutral pions. 

• Function F[y n , 3] is the fit to the photon spectrum generated using F[rj, 3] as a 
weighting function. 
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The direct-photon components of the photon weighting function were determined 
from the PHENIX measurements of the direct-photon spectra for p + p and Au + Au 



collisions. The low-pr 122 and high-py 123 cross-section measurements for p + p 



collisions are shown in Figure 7.4 Two fits to the p + p cross-section were studied. 



Function F[jd, 1] is a fit to the cross-section which was generated by the author 
of this dissertation. 

Function F[jd, 2] is a fit to the cross-section generated by the PHENIX collab- 
oration 122 . This function is used as the primary weighting function for direct 



photons for this analysis. 
Additional direct-photon weighting functions were generated from PHENIX direct 



photon cross-section measurements in Au + Au collisions. Figure 7.5 shows low 



Pt |122| and high-pr 33 measurements of the direct-photon yields in Au + Au 
collisions in three centrality classes: 0-92% (minimum bias), 0-20%, and 20-40%. 

These spectra were fit for pt > 4 GeV/c to determine three additional direct- 
photon weighting functions. 

• Function F[7£>,3] is a fit to the minimum-bias spectrum. 

• Function F[7^,4] is a fit to the centrality 0-20% spectrum. 

• Function F[jd,5] is a fit to the centrality 20-40% spectrum. 

In Au + Au collisions, an additional component of the direct-photon yield has been 



observed 122 . This component, which is most important at low pt, is believed to 



be thermal photons produced in a quark-gluon plasma. The three cross-section mea- 



surements shown in Figure 7J3 were also used to generate three weighting functions 
describing the thermal-photon component. To describe the shapes of the spectra, the 
fit of the p + p direct-photon spectrum (F[jd, 2]) is scaled by the nuclear thickness 
function (T^a)^] A decaying exponential is added to account for the thermal-photon 
component; only this component was allowed to vary during fitting. Three fits to the 
low-pr Au + Au direct-photon spectra were generated. 

1 (T AA ) = 6.14 mtT 1 for the 0-92% centrality class, 18.6 mb -1 for the 0-20% centrality class, and 
7.07 mtr 1 for the 20-40% centrality class [30l. 
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• Function F[jt, 1] has the form Cexp(— p^/T)). The minimum-bias spectrum 
was fit with a function of the form (Taa) x {F[yD, 2] + F[jt, with only the 
parameters of F[jt, 1] allowed to vary. This function is used as the primary 
weighting function for thermal photons for this analysis. 

• Function -Fpyr, 2] was generated in the same fashion by fitting the spectrum in 
the 0-20% centrality class. 

• Function F[jt, 3] was generated in the same fashion by fitting the spectrum in 
the 20-40% centrality class. 

The contributions of photons from //-meson decays, direct photons, and thermal 
direct photons are added together and used as a weighting function for data set P3. 



Figure 7J3 shows the eleven different weighting-function components used for data 
set P3: three to describe photons from ?7-meson decays, five to describe non-thermal 
direct photons, and three to describe thermal photons. 
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RHIC p+p direct-photon spectra 
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Figure 7.4: STAR 117 and PHENIX 122,123 measurements of direct-photon spec 



tra in p + p collisions: Fits to the combined PHENIX spectra are also shown. 
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Figure 7.5: PHENIX measurements 



33 



122 



of direct-photon spectra in Au + Au 



collisions for three centrality classes: The solid curves are fits to the spectra above 
Pt = 4 GeV/c. The dashed blue curves are the function F[jd, 2], a fit to the PHENIX 
p + p direct-photon spectrum, scaled by the value of (Taa) for each centrality class. 
The other dashed fit curves have the form of decaying exponentials, which describe 
the thermal-photon component, added to the scaled p + p fit. 
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y Weighting Function Components 
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Figure 7.6: Components of the photon weighting function used in this analysis: The 
components used to generate the primary weighting function are marked with aster- 
isks. 
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7.3 Background Rejection Efficiency for 
Data Set i\: 7r° — ► 77 Decays 

In this section, various calculations^] of the background rejection efficiency for data 
set Pi are presented. Unless otherwise stated, the weighting function was F[7r°,2]. 
Figure |7.7 shows the efficiency for five different centrality classes. No dependence of 



the background rejection efficiency upon collision centrality is observed. Figure 7.8 



shows Eb for two different types of underlying events: events with trigger id 66007 
(the minimum-bias trigger) and events with trigger id 66203 (high-tower triggered 
events). No strong dependence of £5 on the trigger id of the underlying event is 
observed. Taken together, these results support the assumption that the background 
rejection efficiency does not depend on the characteristics of the underlying event in 
which the photonic e is observed. 



Figure 7.9 shows Eb for that hit a tower which satisfied the high-tower trigger 
condition. Above 3 GeV/c transverse momentum, Eb does not appear to depend on 
whether or not the trigger condition was satisfied. The high-tower trigger condition 
removes almost all particles with transverse momentum less than 3 GeV/c, so there 
are not enough particles in the calculation to determine whether or not Eb is affected 
by applying the trigger condition. Low-p^ e from the real high-tower triggered 
data set are excluded from the calculation of the non-photonic e ± yield and nuclear 
modification factor. 



Figure 7.10 shows eb for various combinations of partner and pair selection cuts. 



Figure 7.11 shows the effect upon the background rejection efficiency of various 
values of the pT(partner) cut. 



Figure |7.12| shows the effect upon the background rejection efficiency of various 
values of the invariant-mass cut. 



Figure 7.13 shows the effect upon Eb of various values of the pair DCA cut. 



Figures 7. 14 7.15, and 7.16 provide a demonstration of the effect of various weight- 



ing functions on the background rejection efficiency. Figure [7. 14| shows a set of exam- 
ple weighting functions used in efficiency calculations. In addition to weighting func- 



tion F[tt , 2], several exponential weighting functions are shown. Figure 7.15 shows 



Eb for each of these functions, as well as the un- weighted efficiency. Figure 7.15 shows 



2 See the note on statistical uncertainties in Section 



B.l 
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the efficiency calculated for each of the five pion weighting functions defined in Sec- 
tion 7J2 The systematic uncertainties in the background rejection efficiency due to 
different weighting functions will be explored in Section 7.6 



MB trigger, Background Rejection Efficiency (e B ), weighted 
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Figure 7.7: Background rejection efficiency as a function of px for five centrality 
classes: The efficiency does not appear to depend strongly (if at all) upon the collision 
centrality of the underlying event. 
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200-GeV Cu+Cu, MB trigger, Back. Rej. Eff. (e B ), cent. 0-60% 
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Figure 7.8: Background rejection efficiency as a function of px for the two different 
STAR trigger IDs considered in this analysis: The efficiency does not appear to 
depend upon the trigger ID of the underlying event. 
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Figure 7.9: Background rejection efficiency as a function of px'- The efficiency is cal- 
culated for tracks that point to towers which satisfied the high-tower trigger condition 
(red). Also shown are efficiency calculations when that condition is not required. 
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Background Rejection Efficiency (e B ) for Various Cuts 
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Figure 7.10: Background rejection efficiency as a function of pr for various combina- 
tions of partner and pair selection cuts: For partner reconstructed, the true conversion 
partner shares at least one half of its TPC points with a reconstructed track, but no 
other cuts are applied. 
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Figure 7.11: Background rejection efficiency as a function of f° r various values of 
the pT(partner) cut: All other cuts have been applied. 
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Figure 7.12: Background rejection efficiency as a function of pr for various cuts on 
the invariant mass of the reconstructed photonic e~e + pair: All other cuts have been 
applied. 
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Figure 7.13: Background rejection efficiency as a function of for various cuts on 
the distance of closest approach between the two reconstructed photonic e ± tracks: 
All other cuts have been applied. 
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Figure 7.14: Different 7r° weighting functions used in the calculations of eb in Fig- 
[TTTol 
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Figure 7.15: Background rejection efficiency as a function of pt for different weighting 
functions. 
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Effect of Weighting Functions on e E 
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Figure 7.16: Background rejection efficiency as a function of pt, with five different 



fits to measured 7r° spectra (see Figure 7.1) used as weighting functions. 



138 



7.4 Background Rejection Efficiency for 
Data Set 7r° Dalitz Decays 

In this section, calculations of the background rejection efficiency for data set P2 are 



presented. Figure 7.17 shows the background rejection efficiency calculated for data 



set P2 for photonic e ± from three sources. 

• conversions of photons from the ir° — > 77 decay 

• conversions of the real photons produced in 7r° Dalitz decays 

• produced directly in 7r° Dalitz decays 

In order to increase the number of e ± in the calculations, the efficiencies shown in 
Figure 7.17| were calculated without BEMC cuts and with a different pseudorapidity 



cut (\rj\ < 0.7). For comparison, Figure 7.17 also shows the background rejection 



efficiency for data set P\ (calculated with the same set of cuts). The calculations 
of Eb for the two-photon decay mode appear to agree for the two simulation data 
sets Pi and Pi- At high p?, Eb for Dalitz-decay e appears to be higher than the 
efficiency for e from photon conversions. At low px, the opposite appears to be true. 



Figure 7.18 shows calculations of Eb (with the standard set of cuts) for electrons from 



7T° Dalitz decays along with the efficiency calculation for data set P\ for comparison. 
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Background Rejection Efficiency for tc — >yy and n Dalitz decays, No BEMC Cuts 
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Figure 7.17: Background rejection efficiency as a function of px for various sources 
of e : To increase the numbers of tracks in the calculation, no cuts using the BEMC 
have been applied. 
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Figure 7.18: Background rejection efficiency as a function of px for various sources of 
e ± \ the standard set of cuts has been applied. 
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7.5 Background Rejection Efficiency for 
Data Set P3: Other Sources 

Data set P3 is used to determine the effect on the efficiency of photonic e ± sources 
other than pions. The photon weighting function has three components: 

• (7^): photons from 77 — > 77 decays 

• (jd)'- direct photons (non-thermal) 

• (7t): thermal direct photons 

Unless otherwise noted, these components are described by functions F[j v ,2], 
F[j£>, 2], and F[y T , 1], respectively The relative strengths of these components should 
be different for Cu + Cu collisions than for p + p or Au + Au collisions. For this 
analysis, the combination Rc uCu (lv) + {Id) + 0-5(7t) was primarily used. The values 
of -RcuCu are assumed to be the same as the pion nuclear modification factor at high 
which varies with the collision centrality: -R^uOu 0-20% = 0-64, -RcuCu 20-60% = 



111 



Pt 

0.94, and -RcuCuo-60% = 0-75. These values of the nuclear modification factor will 



be allowed to vary by ±50%. The PHENIX collaboration has shown 33 that non- 
thermal direct photons are not suppressed in Au + Au collisions with respect to p + p 
collisions. This analysis assumes that the direct photon yield is also not suppressed 
in Cu + Cu collisions. A factor of 0.5 is used as a rough estimate of the suppression 
of thermal photons in Cu + Cu collisions relative to Au + Au collisions. However, 
the thermal-photon contribution is only important at low pt and does not affect the 
background rejection efficiency for the transverse-momentum range of interest for this 
dissertation. 



Figure |7.19| shows the background rejection efficiency calculated separately for 
each of the eleven alternate photon weighting function components (shown in Fig- 
ure 



7.6) 



Figure 7.20 shows the background rejection efficiency for mixtures of photon 
weighting-function components, with the relative strengths of the three components 
varied. 
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Figure 7.19: Background rejection efficiency for the photon weighting components 
shown in Figure 7.6 Also shown is the efficiency calculated using the primary mixture 



of these components. 
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Figure 7.20: Background rejection efficiency for different mixtures of photon weighting 
function components, accounting for photons from ?7-meson decays, direct photons, 
and thermal direct photons. 
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7.6 Combined Background Rejection Efficiency 



The background rejection efficiency for all sources is calculated by adding together 
the reconstructed and rejected e ± yields from 7r° — > 77 decays (data set Pi), tc° Dalitz 



decays (data set P 2 ), and other photon sources (data set P 3 ), as shown in Figure 7.21 
For each data set, the simulated e ± spectra are scaled by the number of simulated n° 
or photons with 1 GeV/c < pt{tt°, 7) < 15 GeV/c. Photonic e from tt° — > 77 decays 
(data set Pi) dominate the combined spectrum at low p?, while the contribution from 
direct photons (in data set P3) is larger at high p?. 



Figures 7.22 and 7.23 show the background rejection efficiency calculations (for 
the 0-60% centrality class) for various mixtures of photonic-e 111 sources. Calculations 
with a greater direct-photon contribution tend to have a higher efficiency, due to 
the fact that the (non-thermal) direct-photon spectrum has a less steeply falling 
slope than the neutral-pion and rj- meson spectra. The effects on eg of the 7r°-Dalitz- 



decay and thermal-photon components are negligible. Figure 7.24 shows fits to the 



efficiency calculations shown in Figure 7.23 for the transverse-momentum range 2 



GeV/c < px < 11 GeV/c. The fits have the form 

C + Aexp sinh sinh {b{j>T — po)) ■ (7.1) 



The fit parameters for the primary efficiency calculation are given in Table [7TT] for 
three centrality classes. The systematic uncertainty curves are obtained by adding or 
subtracting 2% from the primary mixture (black curve). These uncertainties include 
the fits for almost all of the alternate mixtures shown. The exceptions are the mixture 
with no direct photon contribution and the mixtures with meson contributions much 
smaller than expected. 



Figure [7.25| shows the systematic variation in eb for different mixtures of all of 
the alternate weighting function components. The efficiency is calculated for each of 
the following seven mixtures of components, with the relative strengths of the meson 
(7r°,7^), direct photon (715), and thermal direct photon (7^) contributions varying 
within reasonable limits. 

. P£ uCu (7T°,7r,) + (7D)+0.5(7T) 

. 0.5P5 uCu (7r ,7,) + (7D) + 0.5(7T) 



1.5P£ uCu (7r ,7,) + (7D) + 0.5(7r) 
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• £WA7,) + 0.5(7d) + 0.5( 7 t) 

• i?SuCu(A7,) + 2( 7 z>) + 0.5( 7 T) 
. i?£ uCu (7r )7r ,) + ( 7D ) + 0.25(7T) 

• ^CuCuK,7r,) + (7D) + (7T) 

For each of these mixtures, the efficiency is calculated with all permutations of 



the alternate weighting function components described in Section 7.2 



• five possible neutral-pion weighting functions 

• three possible weighting functions describing ?7-meson-decay photons 

• five possible weighting functions describing non-thermal direct photons 

• three possible weighting functions describing thermal direct photons 



Figure 7.25 shows the results of all 1575 of these efficiency calculations; fits to each 



result are shown in Figure 7.26 Figure 7.27 shows the fits of eb for three centrality 
classes; these fits are used to find the efficiency-corrected non-photonic e ± yield (see 
Chapter g). 



144 



Simulated Photonic e 1 Spectra, Mixture 0.75(jt°,y n )+1 .0(y D )+0.5(y T ), cent. 0-60% 
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Figure 7.21: Spectra of reconstructed and rejected photonic e from all three simu- 
lation data sets. 
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Figure 7.22: Background rejection efficiency for various mixtures of simulated pho- 
tonic e ± : Completely removing the Dalitz-decay or //-meson contributions has little 
effect on the efficiency. 
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Background Rejection Efficiency (e B ), Various Mixtures, cent. 0-60% 



1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0, 



~i 0.75(n"T„)+1.0(7 D ) + 0.5(y T ) 
0.00(n°,y„) + 1.0(y D ) + 0.5(y T > 
X 0.20(it,y„) + 1.0(y D ) + 0.5(y T ) 
0.38(ir°,y„) + 1.0(y D ) + 0.5(y T ) 
• 1.12(/, Tn )+1.0(Y D )+0.5(Y T ) 
0.75(it°,y,) + 0.0(y D ) + 0.5(Y T ) 
0.75(ji°,y„) + 0.5(y d ) + 0.5(y t ) 
0.75(n ,y„)+2.0(y d )+0.5(y t ) 
0.75(1 ,y„)+1.0(y d )+0.0(y t ) 
0.75(/,Y>1.0(Yn)*1-0(TT) 



_L 




10 



12 



14 

p r [GeV/c] 



Figure 7.23: Background rejection efficiency for various mixtures of simulated pho- 



tonic e 



± 



Fits of Background Rejection Efficiency (e B ), cent. 0-60% 
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Figure 7.24: Fits to the background rejection efficiency calculations shown in Fig- 



ure 



7.23 Also shown are the primary efficiency measurement (points) and the ±2% 



estimate for the systematic uncertainty of the fits. 
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Systematic Variations in e B , MB Trigger, cent. 0-60% 
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Figure 7.25: Background rejection efficiency (for the 0-60% centrality class) calcu- 
lated for various permutations of weighting functions and various mixtures of func- 
tions: The primary calculation is shown in black; the dashed curves are the estimated 
systematic uncertainties to the primary fit. 
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Figure 7.26: Fits to the measurements of eb (for the 0-60% centrality class) shown 



in Figure 7.25 The ±2% systematic uncertainty on the primary fit includes nearly 



all of the fits shown here. 
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Table 7.1: Fit parameters for the background rejection efficiency (primary calculation) 
for three centrality classes. The parameter b has units of (GeV/c)" 1 ; po has units of 
GeV/c. 



Parameter 


0-20% 


20-60% 


0-60% 


C 


0.68025 ± 0.00644 


0.68601 ±0.00679 


0.68215 ±0.00579 


A 


-0.25726 ± 0.09650 


-0.19389 ±0.04691 


-0.21004 ±0.04182 


b 


-0.34214 ±0.07763 


-0.35903 ± 0.04248 


-0.35923 ± 0.03442 


Po 


0.87181 ±0.57172 


1.1478 ±0.3312 


1.0928 ±0.2731 


X 2 /NDF 


21.9/9 


13.3/9 


15.9/9 



Centrality Dependence of e B Fits 

>,0.75 r 




p r [GeV/c] 



Figure 7.27: Fits to the measurements of sb for three centrality classes. The dashed 
curves indicate the ±2% systematic uncertainty assigned to these fits. 
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Chapter 8 



Residual Background 



This chapter describes the calculation of the residual e background. Residual pho- 
tonic come from the decays of mesons more massive than pions, which may produce 
e~e + pairs with M inv (e~e + ) > 150MeV/c 2 . The yield of residual photonic e ± from 



light-flavor sources (denoted Blfp) is estimated in Section 8.1. The yield (-BlfaO 
of residual non-photonic e , which come predominantly from the decays of kaons, is 



estimated in Section 8.2 Decays of J/tp and T and the Drell-Yan process also con- 
tribute to the residual photonic e ± background; estimates of the e ± spectra from those 



sources (denoted Bj/^, Br, and Bd-y, respectively) are described in Section 8.3 



8.1 Residual Photonic e 

from Light-Flavor Sources 



The analysis method described in Chapter |3| identifies photonic e~e + pairs with in- 
variant mass less than 150MeV/c 2 , allowing the removal of photonic e ± from photon 
conversions and tt° Dalitz decays from the e ± sample. However, other sources pro- 
duce photonic e~e + pairs with M inv .(e~e + ) > 150MeV/c 2 . The p°, u>, or <fi mesons 
may undergo two-body decays to e~e + pairs. Those vector mesons, as well as the 
T], rf, and Kg mesons, may undergo three- and four-body decays which include e~e + 
pairs among the products; some of those pairs may have M ira ,.(e _ e + ) > 150MeV/c 2 . 
In this section, the yield of residual photonic e from light-flavor sources [B^pp) is 



estimated using PYTHIA 75 simulations. 

In these simulations, 10 7 of each particle (tt°, K s , t], p°, u, rj', and <fi) are gener- 
ated. These simulated particles are distributed uniformly in transverse momentum 
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Table 8.1: Simulated decays studied to estimate the residual photonic e ± background. 

Photonic e~e + pairs are removed 



17 



The values of the branching ratios are from 
from the e by the invariant mass cut; simulations of ir° are included in this chapter 
only for reference. There are two decays for which only an upper limit is given; the 
value of the branching ratio is assumed to be that upper limit, with a fractional 
uncertainty of ±100%. Regarding the rf — > 7p° and — > 7r°p° decays: the branching 

instead of p°. Those values 



17 include non-resonant modes with 7r 



7T 



ratios given m 

are used here, although they are likely overestimates of the branching ratios for decays 
to p°. 



Decay 


Branching Ratio 


7T° — » e~e + 7 


(1.174 ±0.035)% 


Kg e~e + n~n + 


(4.69 ±0.30) x 10~ 5 


r] — > e~e + 7 


(7.0 ±0.7) x 10" 3 


7] —7- e~e + 7r~7r + 


(2.68 ±0.11) x 10" 4 


p° ->■ e"e + 


(4.72 ± 0.05) x 10~ 5 


p° — >• e~e + 7r° 


< 1.2 x 10- 5 (CL=90%) 


a; — ?• e~e + 


(7.28 ±0.14) x 10" 5 


u —7- e~e + 7r° 


(7.7 ±0.6) x 10" 4 


r^/ —7- e~e + 7 


< 9 x 10" 4 (CL=90%) 


r^/ —7- e _ e + 7r _ 7r + 


(2.4+^) x 10" 4 


—7- e _ e + 


(2.954 ± 0.030) x 10" 4 


— > e~e + 7T° 


(1.12 ±0.28) x 10~ 5 


(ft —7- e~e + r] 


(1.15 ±0.10) x 10" 4 


r/f —7- 7Y 7r°rj 


(21.7 ±0.8)% 


1]/ —7- 7T~7T + r] 


(43.2 ± 0.7)% 


7]! —7- 7p° 


(29.3 ± 0.5)% 


—7- 7W 


(2.75 ± 0.22)% 


— t- 7?7 


(1.309 ± 0.024)% 


7T p° 


(15.32 ±0.32)% 



(OGeV/c < Pt < 30GeV/c), rapidity (\y\ < 1.5), and azimuth (0 < < 2n). These 
particles are allowed to decay via modes that produce e~e + pairs. Also simulated are 
the decays of r/ and mesons to r], p°, and u> mesons, with the subsequent decays of 



those particles to e e + pairs. Table 8.1 summarizes the decays simulated, along with 



the physical branching ratios 17 



As was the case in the calculation of the background rejection efficiency (see Chap- 
ter [7]), the spectra produced in these simulations must be weighted to correct for the 
fact that the decaying particles are generated uniformly in transverse momentum. The 
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Table 8.2: The values of R(X/n°) used to scale the meson weighting functions. The 
Kg ratio is estimated from |124| . The other ratios are the same as those used in [55] 
and [80]. 



Ratio 


Value 


K s /ir° 


0.6 ±0.3 


77/7T 


0.48 ±0.10 


p°/vr 


1.0 ±0.3 




0.90 ±0.27 


r)f/ir° 


0.250 ± 0.075 




0.4 ±0.12 



7T° and 77-meson weighting functions that were used in Chapter [7] (functions F[7r°,2] 
and F[i], 2], respectively) are used as weighting functions here. For the other particle 
species, the shapes of the weighting functions are estimated from the 7r° weighting 
function by making the substitution pr(^°) ^/pt(X) 2 + m(X) 2 c 2 — m(ir°) 2 c 2 (my 
scaling), where X denotes the various particle species. In addition, the weighting 
function for each particle species is scaled by constant factors according to Equa- 
tion 8.1 to give the invariant yield (Ed 3 N/dp 3 ) of particle species X in each of the 



three Cu + Cu centrality classes. 



F[X, Cu + Cu] = E d3N °l +Cu = R* uCu x x R(X/n°) x F[X\. (8.1) 

Here, F[X] denotes the m^-scaled ir° weighting function and R(X/ir°) denotes the 
ratio of the measured yield of particle species X to the ir° yield. The values of 



R(X/tt°) are given in Table 8.2 the PHENIX and STAR collaborations use the same 
ratios in their calculations of photonic e spectra in 200-GeV p + p and Au + Au 
collisions. (Those analyses do not consider the K s contribution; the choice of 0.6 ±0.3 



as the kaon-to-pion ratio is explained in the next section.) In Equation 8.1, -RJuCu 
indicates the high-p^ charged-pion nuclear modification factor measured by the STAR 
collaboration in 200-GeV Cu ± Cu collisions [111] . It is assumed that neutral pions 
will have the same suppression as the charged pions. The values of -Rc uCu and (Nun) 
(the number of binary nucleon-nucleon collisions) for the three Cu + Cu centrality 



classes are given in Table 8.3 The inelastic p + p cross-section, a ine i, is 42 mb |17|. 
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Table 8.3: High-p^ charged-pion nuclear modification factors and values of (A^ n ) [ill 
for various 200-GeV Cu + Cu centrality classes. 



Centrality 


D7T 

-"-CuCu 


(N Hn ) 


- 20% 


0.66 ±0.10 


86.84+l^ 


20 - 60% 


0.99 ±0.14 


OO.U±_ 48 


- 60% 


0.82 ±0.12 


51 3n+ u - YS 



Each of the simulated photonic e is weighted by the branching ratio for the 
process that produced it. Each e ± is also weighted according to the transverse mo- 
mentum of its parent particle by a factor of px(X) x F[X, Cu ± Cu] . For each particle 
species X, the factor Pt(X) is needed because 

d 3 N 1 d 2 N 
F[X, Cu ± Cu] = E— oc — x — — . (8.2) 
dp 6 pt dprdy 



M inv (e"e + ) Distributions, 3.0<p_<4.0 
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Figure 8.1: Invariant mass of photonic e e + pairs from various meson-decay processes 
(3GeV/c<p T <4GeV/c). 



Figure 8.1 shows invariant- mass distributions (for 3GeV/c < p? < 4GeV/c) for 
the e~e + pairs produced in the simulated decays. Peaks in the distribution due to 
the two-body decays of the vector mesons are visible, as are wide distributions with 
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Photonic (e +e + )/2 from PYTHIA Simulations (cent. 0-60%) 






Pllj.Pl! NPP 




Ixg — ?G e 71 71 




t| — >© e + y 




1 1 7C C Jv Jl* 


o 


— >6"6 + 




. . 
n — ^P 0*77 




CO — >6 6 + 




- * 

fil :»P P TT 


A 


r|'->e"e + Y 


X 


r|'— >e"e + 7t"jt + 


* 


if^X^e"e + X 


o 


(])->e"e + 


A 


(])->e"e + 7t 


□ 


(|)^e"e + T| 


* 


^X^e"e + X 


ft 


Sum 



8 9 10 
P T [GeV/c] 



Figure 8.2: Spectrum of photonic e from pairs that fail the Mi nv _(e e + ) < 
150MeV/c 2 cut, scaled for the 0-60% centrality class of 200-GeV Cu + Cu collisions. 



peaks near due to Dalitz decays. Figure |8.2 shows the simulated invariant yield of 
photonic |(e~ + e + ) for the 0-60% centrality class in 200-GeV Cu + Cu collisions. 
For the spectra shown in this figure, the cut on pair invariant mass has been applied. 
Also shown is the efficiency-corrected invariant yield for non-photonic found in 
this analysis (prior to the subtraction of the residual background). 

The systematic uncertainties in Blfp, the total residual light-flavor photonic e 
yield, are calculated by adding in quadrature the uncertainties in the total spectrum 
due to the uncertainties in the branching ratios, the R(X/ir°) values, Rq uCu , an d 
the uncertainty due to the choice of pion weighting function. As described above, 
function F[ir°, 2] was chosen as the primary 7r° weighting function and the weighting 
functions for the other particle species were derived from F[tt°, 2] by ttit scaling. Four 
other sets of weighting functions were also derived, each set based on a different 7r° 
weighting function (see Chapter 17]). Blfp was calculated using each set of weighting 
functions. In each bin, the difference between the primary calculation and the 
most extreme calculations above (below) was used as an upper (lower) systematic 
uncertainty. These weighting-function systematic uncertainties were then added in 
quadrature to the uncertainties due to the scaling factors. 
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M inv (eV) Distributions, PYTHIA co^e"e f 7?V2.0<p T <3.0 
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Figure 8.3: The effect on the e~e + pair invariant mass distribution of using different 
momentum distributions in the simulation of the decay u — > e~e + ir°. The blue distri- 
bution was generated by allowing the momenta of the decay products to be uniformly 
distributed in the available phase space. The more realistic black distribution was 
generated according to ll 161 (Kroll-Wada). 



Note that the three-body decays of p° and <p may not have been simulated using 
the correct invariant-mass distributions: the products of those decays have momenta 
distributed uniformly in the available phase space. The u — > e~e + n° decay was 



simulated using the Kroll-Wada distribution 116 , as well as the "phase-space" dis- 
tribution; the M im) .(e~e + ) distributions for these cases are shown in Figure 8.3. For 
this decay, the incorrect phase-space distribution has a greater fraction of pairs with 
M inv Xe~e + ) > 150 MeV/c 2 ; using the phase-space distribution has given an overesti- 
mate of the residual background. A similar effect is observed when the phase-space 
distributions are used for the pseudoscalar-meson Dalitz decays. It is assumed that 
the same will hold true for the three-body p° and cf) decays and that the spectra of 
e from those decays shown in Figure 8.2 are overestimates. Note that those decays 
make negligible contributions to the residual background. 

Similarly, the four-body decays of Kg, 7], and rj (none of which are studied in the 
other RHIC non-photonic e ± analyses 55 80j) were also simulated by distributing 
their products uniformly in the available phase space. However, these decays make 



154 



only small contributions to the residual background. Even the contribution from the 
four-body r( decay is only about 5% of Blfp, a contribution that is smaller than the 
uncertainties in the background due to the uncertainty in -RcuCu an d the uncertainty 
due to the choice of weighting function. Note also that for the calculations in this 
section (though not in the next section) the Kg were assumed to have lifetimes similar 
to the other mesons and to decay at the collision vertex. If this has any effect, the 
(negligible) K s contribution to the residual photonic background is an overestimate 
(some of the decay e ± tracks may have failed the GDC A cut if the decays occurred 
far from the vertex). 
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Table 8.4: Simulated decays studied to estimate the residual non-photonic e back- 
ground. The values of the branching ratios come from 17 . 



Decay 


Branching Ratio 




1 


7r + — > e + v e 


(1.230 ± 0.004) x 10~ 4 


K + ->• e+z/ e 


(1.584 ±0.020) x 10~ 5 


K + ->• e + z/ e 7r° 


(5.07 ±0.04)% 


Kg — > e~z/ e 7r + (e + t' e 7r _ ) 


(7.04 ±0.08) x 10" 4 


->• e-^ e vr + (e + i/ e 7r") 


(40.55 ±0.12)% 


A — > e~9 e p 


(8.32 ±0.14) x lO" 4 


S+ e+z/ e A° 


(2.0 ±0.5) x 10" 5 


£~ —7- e~v e n 


(1.017 ±0.034) x 10" 3 


S- e"z/ e A 


(5.73 ±0.27) x 10- 5 


H° e-z/ e S+ 


(2.53 ±0.08) x 10" 4 


2~ -)• e~z/ e A° 


(5.63 ±0.31) x 10" 4 


5" e-p e S° 


(8.7 ± 1.7) x 10" 5 


fl~ — >■ e~z/ e S° 


(5.6 ±2.8) x 10" 3 



8.2 Residual Non-Photonic e 

The non-photonic e sample includes contributions from sources other than heavy- 
flavor decays. The largest such sources are the decays of kaons; there are also contribu- 
tions from charged pions, direct muons, and baryons. The total yield of non-photonic 



e ± from these sources (denoted B LFN ) is estimated in this section using PYTHIA 75 
simulations. In these simulations, 10 7 of each particle (/x - , 7r + , K + , K s , K^, A , X + , 
E~, 2°, E~, and Q~) are generated. These simulated particles are generated uni- 
formly in transverse momentum (OGeV/c < px < 30GeV/c), rapidity (\y\ < 1.5), 
and azimuth (0 < < 2tt). These particles are allowed to decay via modes that pro- 
duce e . Table 8.4 summarizes the decays studied, along with the physical branching 



ratios (17) . 

As was the case in the previous section, the spectra produced in these simulations 
must be weighted to correct for the fact that the decaying particles are uniformly 
distributed in transverse momentum. The same weighting function used for the 7r° 
in the previous section {F[ir°, 2]) is used to describe the charged-pion spectrum here. 
The shapes of the kaon weighting functions were estimated from the pion weighting 
function using my scaling. An m^-scaled pion spectrum was also used as a rough 
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guess of the muon spectrum. However, for reasons described below, the muon-decay 
contribution to the residual background is believed to be negligible, so the use of a 
questionable muon weighting function will have no effect on the background. 



STAR Au+Au proton spectra 



* 0-12% (x10) 

• 10-20% 




8 9 
p r [GeV/c] 

Figure 8.4: The STAR Collaboration's measurements of h(p+p) spectra in 200-GeV 
Au + Au collisions 125 . The spectra are corrected for hyperon (A and S + ) feed- 



down. The error bars indicate combined statistical and systematic uncertainties. The 
fit to the spectrum for the 40-60% centrality class (blue) is used as the basis for the 
baryon weighting functions in this section. 



The STAR collaboration has measured 125 the \{p+p) yield in various centrality 



classes in 200-GeV Au + Au collisions (see Figure 8.4). The 40-60% centrality class 



(for which (N part ) and (Nu n ) have similar values to the Cu + Cu 0-60% centrality 
class) is fit with function F[p,4] of the form A(l + p\/B) n + Cexp(—pT/D) (see 



Section D.4 for the fit parameters). The shapes of the baryon weighting functions are 
estimated from this proton fit using my scaling. 

The weighting function for each particle species is scaled by constant factors ac- 



cording to Equation |8.3| for muons and mesons and according to Equation |8.4| for 
baryons. This gives the invariant yield (Ed 3 N/dp 3 ) of particle species X in each of 
the three Cu + Cu centrality classes. 



F[X, Cu + Cu] = R, 



Oinel. 



[l + Ranti) X F[X\. 



(8.3) 
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F[X, Cu + Cu] = ( ^a u) x R ( X /P) x C 1 + x W ( 8 - 4 ) 

-Rant? is the antiparticle-to-particle ratio for species X. For and K^, which are 
their own antiparticles, R an ti is set to 0. For Q, only the combined (particle plus 
antiparticle) yield is given in |126| ; R an u is set to and the Vt + contribution is 
accounted for in the fi-to-proton ratio. {N^ +An } = 90 is the number of binary 
nucleon-nucleon collisions 127 for the 40-60% centrality class of 200-GeV Au + Au 
collisions (the baryon weighting functions are based on a fit to the \{p + p) spectrum 
for this centrality class). 

The values used for R(tt^/7t + ) and R(K~/K + ) are based on the PHENIX col- 



laboration's measurements of those quantities in 200-GeV Au + Au collisions 124 



The values quoted in Table 8.5 are the means of the measurements for the 40-50% 
and 50-60% centrality classes. These measurements are for p? < 2GeV/c, but no 
transverse-momentum dependence is observed and it is assumed that the ratios will 
retain those values at high px- The value of R(K + /vr + ) is a rough estimate from Fig- 
ure 17 in 124 . The desired value is the K + /tt + ratio at high-p^, for which no data 
are available. The value R(K + /tt + ) = 0.6 is an estimate based on the apparent trend 
of the ratio as pt increases; large uncertainties are therefore assigned. It is assumed 
that R(fi + /fi~) = 1 and that the muon-to-pion ratio is also 1. The latter assumption 
is not intended to be a realistic estimate of the muon yield, but rather a "worst-case" 
scenario. The resulting spectrum of e from muon decays will be discussed below. 

The baryon ratios are the ratios of the STAR collaboration's measurement of the 
integrated rapidity densities (dN/ dy) for each particle species for the 40-60% central- 
ity class in 200-GeV Au + Au collisions. The value of dN/dy for protons is taken 
from |127| , while the values of dN/dy for A, Cl~ , and their anti-particles are taken 
from |126| . For the E baryons, it is assumed that R(T l + /p) = R(E~ /p) = R(A/p) and 
#(£-/£+) = = R(A/A). For H° it is assumed that R(Z°/p) = R(E~/p) 

and R(E°/E°) = R(E+/E~). 

As in the previous section, each of the simulated non-photonic e is weighted 
by the branching ratio for the process that produced it and by a factor related to 
the transverse momentum of its parent particle, pr{X) x F[X, Cu + Cu]. Unlike the 
previous section, the parent particles here may have long enough lifetimes that a 
significant number of decays may occur farther than 1.5 cm from the primary vertex. 
It is possible that an produced in such a decay will fail the GDC A cut and not 
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Table 8.5: Particle ratios used to scale the weighting functions for the calculation of 
the residual non-photonic background. Meson ratios are from 
are from 



127 



124 



proton yields 



and other baryon yields are from |126| . Values in parentheses are 
assumptions and not based on independent measurements. See the text for further 
explanation. 



R{X/tt+) 
or R(X/p) 


Value 


Ranti 


Value 






7T~ /tT + 


0.97 ±0.06 


K+/TI + 


0.6 ±0.3 


K-/K+ 


0.92 ±0.05 


A/p 


0.206 


A/A 


0.792 


E+/p 


(0.206) 


E~/E+ 


(0.792) 


E"/p 


(0.206) 


E+/E" 


(0.792) 


E°/p 


(0.259) 




(0.885) 


E lp 


0.259 




0.885 


(Q- + Q+)/p 


0.00627 







be included in the e sample. For each parent particle X, a decay point along its 
track is randomly generated. Using that decay point a helical track is constructed 
for the produced e , the DCA to the primary vertex is calculated, and the cut of 
GDC A < 1.5 cm is applied. 

To generate a random decay point, a decay length L is generated with a uniform 
distribution in the range < L < A^, with A^=3 m for mesons and muons and 
A^=4.7 m for baryons. Those values were chosen so that a K + with p? = 

lGeV/c traveling the full distance Al along a helical path in the STAR magnetic 
field (0.497952 T) would always decay a radial distance > 102 cm from the z-axis, 
even when the particle had the extremal value of rapidity, y = 1.5. However, the 
decays are still constrained to be within a few meters of the STAR detector, rather 
than the large distances possible if the physically correct exponential distribution of 
decay lengths was used. For each decay, where the decaying particle had mass m, 
lifetime r, and momentum p, a decay probability D was calculated and used to scale 
the uniform decay-length distribution to match the physical exponential distribution. 

„ . , mc 2 / mc 2 L \ 

This is demonstrated in Figure [8~5j In black is an exponential distribution of decay 
lengths (with decay constant m/rp =0.47 m _1 ), while in blue is a uniform distribution 
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Simulated Decay-Length Distributions 



^ 0.5 
E 

j 0.45 
| 0.4 
0.35 
0.3 
0.25 
0.2 
0.15 
0.1 
0.05 
0, 



0.5 



-| — Exponential Dist. 
□ Uniform Dist. 
° Scaled Uniform Dist. 



1 1 M 1 1 n 1 1 1 n i m 1 1 1 n 1 1 1 1 1 1 ( t ( 1 1 1 



1.5 



2.5 



3.5 



4.5 5 
L[m] 



Figure 8.5: Demonstration of how a uniform distribution is scaled to match an expo- 



nential distribution using Equation 8.5 See the text for explanation 



of decay lengths (with = 3m). Both distributions are normalized by the number 
of entries (5 x 10 6 ) and the bin width. The red distribution is the blue uniform 



distribution scaled by the decay probability D calculated using Equation |8.5| The 
uniform distribution scaled by D matches the exponential distribution for L < A^. 
Each simulated e ± is weighted by the decay probability D for the decay that produced 
it. 

The residual non-photonic e background (Blfn) calculated using the method 



described above is shown in Figures 8J3 and 8/7 Figure 8J3 shows the background 
spectra for e from meson and muon decays in the three Cu + Cu centrality classes. 
This spectrum is dominated by the K + — > e + i , e Tr decay (called K^ 3 ) and the K® — > 
e + v e n~ decay (called K® 3 ) and their charge conjugates. 

The systematic uncertainty in the total background non-photonic spectrum 
from meson decays is the sum in quadrature of the uncertainty in each bin due 
to the uncertainty in R(K + /vr + ) and the uncertainty due to the choice of weighting 
function shape. As in the previous section, four additional sets of weighting functions 
were generated, each set based on a different n° weighting function. Residual non- 
photonic e spectra were calculated using each set of weighting functions. In each pt 
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Non-Photonic (e +e*)/2 from Non-Heavy-Flavor Mesons (cent. 0-60%) 



10" 1 




* 


Cu 


+Cu NPE 


• 


7I + 


-> e + v e + c.c. 




K + 


-> e + v e + c.c. 




K + 




— > e v e 7i + c.c. 






-> e"v e 3i + + c.c. 


□ 


K° L 


-> e'v e Ti* + c.c. 




Mesons 




H" - 


-> e"V e v„ + c.c. 


* 


Baryons 



9 10 
[GeV/c] 



Figure 8.6: Spectrum of non-photonic e ± from decays of mesons and muons, scaled 
for the 0-60% centrality class of 200-GeV Cu + Cu collisions. 



bin, the difference between the primary calculation and the most extreme calculations 
above (below) was used as an upper (lower) systematic uncertainty. These weighting- 
function systematic uncertainties were then added in quadrature to the uncertainties 
due to the R(K+/n + ). 



The PHENIX collaboration has measured 49 the low-p T direct-muon yield at 
forward rapidity (1.5 < rj < 1.8) in 200-GeV p + p collisions; those values are < 1% 



of the mid-rapidity pion yields 120 . The main sources of feed-down muons should 
be 7r , K , and K\ decays but those particles have long enough lifetimes that (with 
p > lGeV/c) most of them will decay outside of the TPC. Even in the extreme 
scenario in which the spectrum of direct muons is the same shape and magnitude as 
the charged-pion spectrum, the yield of e from the decays of those muons would 
amount to less than 10% of the total yield of e ± from pion and kaon decays and less 
than 1% of the efficiency-corrected non-photonic e yield. The simulated spectrum 
of e ± from muon decays is therefore not included in the estimates of B^p^. 



Figure 8/7 shows the background spectra for e from baryon decays for the 0-60% 
centrality class in 200-GeV Cu + Cu collisions. The baryon contribution is many 



161 



orders of magnitude below the contribution from the K eS decays and the efficiency- 
corrected non-photonic e signal. The baryon contribution is therefore disregarded. 



Non-Photonic (e +e*)/2 from Non-Heavy-Flavor Baryons (cent. 0-60%) 
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Figure 8.7: Spectrum of non-photonic e ± from decays of baryons, scaled for the 0-60% 
centrality class of 200-GeV Cu + Cu collisions. 
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8.3 Other Sources 



This section describes the calculation of the residual photonic e ± background due 
to J/ip and T decays and the Drell-Yan process. This analysis estimates the yields 
of e ± from those sources (denoted Bj/^, By, and B D _ Y ) using the same estimated 
spectra that the STAR collaboration uses in its measurement [80] of non-photonic e ± 
in p + p collisions. Those spectra are scaled by a nuclear modification factor -RcuCu 
(for the quarkonia) and by (N bin ) / a ineL to convert them to an expected yield in Cu + 
Cu collisions. The spectra of photonic e from J/ ip and T decays and the Drell-Yan 



process will be shown in Figure 9.2| (page 172); the calculations that produced those 
spectra are described below. 



8.3.1 Photonic e from J/i{> Decays 



JAj/ in Cu+Cu Collisions 



STAR 0-20% 
PHENIX 0-20% 
Fit 0-20% 
PHENIX 20-40% 
PHENIX 40-60% 
Combined 20-60% 
Fit 20-60% 




9 10 
p r [GeV/c] 

Figure 8.8: Spectra of J /ip in various centrality classes of 200-GeV Cu + Cu collisions 



measured by the STAR [47] and PHENIX 128 collaborations. Also shown are fits 
to these spectra. 



The decay J/tp —¥ e e + , which has a branching ratio of (5.94 ± 0.06)%, [TF] 
produces e _ e + pairs that fail the invariant- mass cut. While these e do come from 
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JA|/ Nuclear Modification Factor in Cu+Cu Collisions 
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Figure 8.9: The nuclear modification factor measured 128 by the PHENIX collab- 
oration for J /ip in Cu + Cu collisions for the 0-20% (circles) and 20-60% (squares) 
centrality classes. Also shown in the STAR collaboration's measurement [47] of the 
nuclear modification factor for > 5 GeV/ c in the 0-60% centrality class. The curves 
indicate the parametrized nuclear modification factor for the 20-60% centrality class; 
the upper dashed curve is used to estimate the upper limit of the J/tp yield in that 
centrality class. 



a heavy-flavor source, they are photonic and must be subtracted from the efficiency- 
corrected non-photonic e spectrum. The yield of from J/ip decays is estimated 
using PYTHIA simulations of J/ ip production in 140 million p + p collisions, with 
the subsequent decays of the J /ip to e~e + pairs. The shape of the simulated J /ip 
spectrum as a function of transverse momentum is not the same as the shape of 
the spectrum of J/ip in real Cu + Cu collisions; the simulated J/ij) are weighted to 



account for this. Figure 8.8 shows measurements 47 128 of the J/ip yields for Cu + 



Cu collisions in several centrality classes. For the 0-20% centrality class, the combined 
STAR and PHENIX data are fit with a function F[J/ip 7 - 20%], which has the form 



F(p T ) = A[exp(-ap T - bp 2 T ) +p T /Po 



(8.6) 
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J/y Spectra: Real Data and PYTHIA 
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Figure 8.10: Upper panel: The curves represent a fit (with systematic uncertainties) 
to the spectrum of J/ip in the 0-20% centrality class of Cu + Cu collisions. The 
black histogram represents the cross-section of J/ip in p + p collisions simulated using 
PYTHIA. Both quantities have been multiplied by the branching ratio for the J/ip — > 
e~e + decay. Lower panel: The ratio (Rpythia) of the fit to the simulated PYTHIA 
spectrum; the gray lines indicate the ratio of the dashed curves to the simulated 
spectrum. 



The fit parameters are given in Section D.5 The dashed curves, which indicate 
systematic uncertainties of the fit, are derived from the calculated covariance matrix]]] 

lr rhe calculated covariance matrix is used to evaluate the upper and lower limits of the 68% 
confidence interval at many values of pt', the results are then fit to obtain the curves describing the 
lower and upper limits. 
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The J/ip spectrum for the 20-60% centrality class is parametrized by a function 



F[J/tp, 20 — 60%], which has the form given by Equation 8.6 The parameters of 
F[J/ip, 20 — 60%] are found as follows. The J/ip spectrum for the 20-40% centrality 
class is fit with a function F[J/ip, 20—40%] (which also has the form given by Equation 



8.6). The values of the parameters a, 6, p , and n in FfJ/^,20 — 60%] are set equal 
to the values those parameters have in F[J/ip, 20 — 40%] (i.e., F[J/ip, 20 — 60%] is 
constrained to have the same shape as F[J/ip, 20 — 40%]); the value of the scale 
parameter A for F[J/ip, 20 — 60%] is found by x 2 minimization. The curve describing 
the lower limit of F[J/ijj, 20 — 60%] is found using the method described previously. 
A different method is used to find the curve describing the upper limit of F[J/ip, 20 — 
60%]: the J/ip nuclear modification factor is parametrized and then multiplied by 



a fit of the J/ip yield in p + p collisions. Figure 8.9 shows the measured values of 



the -RcuCu' the nuclear modification factor for J/ip in Cu + Cu collisions. Sigmoid 
functions of the form 



>j/V • B — A 



^CuCu,20-60%(^) ~ A 



erf ^ + 1 



\.7) 



are used to parametrize the central value and the lower and upper limits ofi^u 

over 

the range of transverse momentum shown. At low p?, the functions are constrained 
to the values of Rq^cu 20-60% measure d by the PHENIX collaboration, while at high 
pr, the functions are constrained to -RcuCuo-60% = •"••^ ^ 0-45, measured^ by the 
STAR collaboration for p? > 5 GeV/ c. The parameters of the function f up used to 
describe the upper limit of -RcuCu 20-60% are ^ = 0-87, B = 1.85, a = 4.2GeV/c and 
b = lGeV/c. The upper limit of the J/ip yield in the 20-60% centrality class of Cu 
+ Cu collisions is taken to be 

F[J/ij;, 20 - 60%+}(p T ) = f up (p T ) ■ (Tju) ■ F[J/^,pp}(p T ), (8.8) 
where the function F[J/tp,pp\ is a fit to the combined STAR (47] and PHENIX (46] 



measurements of the J/ip yield in p + p collisions (see Section D.5 for the fit param- 
eters) . 

For each centrality class, the ratio of the J/ ip fit functions for Cu + Cu collisions 



to the simulated PYTHIA J/ip spectrum is calculated, as illustrated in Figure 8.10 



The upper panel shows the J/ip cross-section simulated using PYTHIA, as well as the 



2 A separate measurement for the 20-60% centrality class is not available. 
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fit function F[J/ip, — 20%] and its systematic uncertainties. The lower panel shows 
RpythiAi the fit divided by the PYTHIA spectrum, with the gray lines indicating 
the the lower and upper limits of the fit divided by the simulated spectrum. For each 
centrality class, the spectrum of e from J/ip decays was calculated by weighting each 
e ± by the value of Rpythia (evaluated at the transverse momentum of the parent 
J/ip). The results of these calculations will be shown in the next chapter. 

8.3.2 Photonic e ± from Y Decays 




There are currently no RHIC measurements of the T cross-section as a function 
of transverse momentum at a/saT/v = 200 GeV. Instead, the shape of the T spectrum 
is estimated using the Color Evaporation Model (CEM). In this model, 129 131 the 
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bb production cross-section is calculated perturbatively at the Next-to-Leading Order 
("NLO", 0(a^)) level for values of Vs, the parton-parton CM collision energy, in the 
range 2 x (mass of b quark) < v3 < 2 x (mass of B ± meson) . It is assumed that any net 
color charge in the bb pair will be neutralized through non-perturbative interactions 
with the color field in the collision region (color evaporation). The b quark and b 
antiquark may then combine with each other to form a bb bound state, or may combine 
with light quarks to form B mesons. However, since V§ is less than the threshold for 
the production of two B mesons, any energy needed to form these mesons is assumed 
to come from non-perturbative interactions with the color field in the collision region. 
The T cross-section^ will be a fraction Fy of the calculated bb cross- sect ion; Fx has 
been found |132| to be independent of y/s, rapidity, and pr- Therefore, Fx need 
not be calculated at 200 GeV, but can be determined from measurements at lower 



collision energies. Figure 8.11 shows a CEM calculation 129 131 of the T cross- 
section (including the T(IS'), T(2S), and T(3S) states) at mid-rapidity as a function 
of transverse momentum. The px-integrated T cross-section estimated using the 



CEM is found to be consistent with the STAR collaboration's measurement 48 of 
that quantity: o r x(is+2S , +3S) = [114±38(stat.)+24(sys.)]pb. Adding the statistical and 
systematic uncertainties in quadrature gives a 39% uncertainty in this measurement. 

PYTHIA was used to simulate proton-proton collisions, producing 4.8 million 
T, which were allowed to decay to e~e + pairs. These PYTHIA calculations do not 
reproduce the shape of the T spectrum given by the CEM, so the PYTHIA data 
are weighted by the ratio of the CEM spectrum to the PYTHIA spectrum. This 
ratio is found to vary between 4 and 0.5, generally decreasing with pt- For each 
Cu + Cu centrality class, the spectrum of e ± from T decays is assumed to be the 
spectrum calculated for p+p collisions, scaled by {Taa) -^cuCu- Since no measurement 
of -Rq uCu is available, it is assumed that the nuclear modification factor is 1, with large 
uncertainties: -Rq uCu = 1 ± 0.7. The results of this calculation are shown in the next 
chapter; the contribution of T decays to the spectrum is small enough that the 
large uncertainties in the nuclear modification factor do not have a large effect on the 
final corrected spectrum. 



3 The CEM can also be used to estimate the production cross-sections for other heavy quarkonium 
states, including J/tp. 
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8.3.3 Photonic eF from the Drell-Yan Process 



A Leading-Order calculation 114 of e~e + pair production via the Drell-Yan process 
(qq 7* e~e + ) is used by the PHENIX collaboration as a correction to their 
measurement |55] of non-photonic e ± in 200-GeVp+p collisions. This same spectrum, 



scaled by (N bin ) /a ineL , is used in this dissertation. 
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Chapter 9 



Results 



9.1 Correction of Spectra 

The normalized efficiency-corrected yields of inclusive^ photonic, and non-photonic 
e± (Iec, Pec, an d N EC respectively) are found as follows. 



K I 



^ EC 



EC 



(Abemc)£r ■ £t ■ £dE/dx 

p 

{AbEMc) £ R • e T ■ £ dE/dx ■ £ B 



(9.1) 
(9.2) 



N EC = Iec ~ Pec = Vy- A ^ ( K inc I - —p) . (9.3) 

\ABEMC)£R ■ £ T ■ £ dE/dx \ £ B J 

Here, I is the uncorrected yield of inclusive e ± and P = n unUke — 2y/n n ++ is the 

uncorrected yield of photonic e ± . The correction factors^] are summarized in Table 3.3 



(page 79). The normalization factor v is 



„ = ^ (94) 

where A<ft = 2tt is the azimuthal range of the measurement, Ay = 0.8 is the rapidity 
range (for ultra- relativistic e ± , Ay is very well approximated by the pseudorapidity 
range), and Ap? is the width of each transverse momentum bin (usually 1 GeV/c). 
The initial factor of 2 accounts for the fact that two particle types, electrons and 



1 all e , the sum of the photonic and non-photonic contributions 

2 The spectra of e^- in the minimum-bias data set are not corrected by the high-tower trigger 
efficiency St- 
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positrons, are included in the measurement. (For the sake of brevity in the following 
discussion, a normalized \{e~ + e + ) yield will continue to be referred to as an e 
yield.) For each centrality class, the number of events (N evt ) is given by Table [3^2 



(page 57), with N^ t B used for the minimum-bias data set and N^J (normalized) 



used for the high-tower-triggered data set. The weighting factor {pj 1 ) is the mean 
value of 1/pt for the inclusive e ± spectrum in each px bin. The factor of (p T 1 )c is 



necessary for the conversion to an invariant yield of the form Ed 3 N/dp 3 . Figure 9.1 
shows the efficiency-corrected e yields as functions of transverse momentum for the 
0-20% most central 200-GeV Cu + Cu collisions. These spectra are combinations of 
the e spectra measured in minimum-bias triggered events (for pt < 4 GeV/ c) and 
high-tower triggered events (for pt > 4GeV/c). 



200-GeV Cu+Cu, corrected, normalized (e"+e*)/2 spectra, cent. 0-20% 
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Figure 9.1: The efficiency-corrected e yields in the 0-20% most central Cu 
collisions: inclusive (black), photonic (red), non-photonic (blue). 



Cu 



The efficiency-corrected non-photonic e yield, Nec, includes some contamina- 
tion. Chapter [8] describes the estimation of the various components of the residual 
background, which is subtracted from ./Vso The yield Y HF of e ± from the decays of 
hadrons containing heavy quarks is 

Y HF = N EC — B LFP — B LFN — B D _ Y , (9.5) 
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(e"+e + )/2 Yields from Various Sources, cent. 0-20% Efficiency -corrected view 
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Figure 9.2: The efficiency-corrected non-photonic e yield {Nec), various components 
of the residual background (see the text in this section and Chapter [8]), and the yields 
of e from heavy-flavor (Yhf) and open-heavy-flavor sources (Yd,b)- The spectra 
shown are for the 0-20% centrality class of 200-GeV Cu + Cu collisions. 



where -Blfp is the residual yield of photonic e from light-flavor sources (rj, p°, u, rj', 



and Kg\ see Section 8.1), Blfn is the yield of non-photonic e from light-flavor 



sources (predominantly the semileptonic decays of kaons, see Section 8.2), and Bd-y 
is the yield of produced via the Drell-Yan process. The yield Y HF includes photonic 
from the decays of mesons made up of cc and bb pairs, predominantly J/ip and 
T. The yields of e ± from those sources are subtracted from Yhf to give Yd,b the 
yield of non-photonic e ± from the decays of D and B mesons, which are said to have 
"open heavy flavor" because they have non-zero net charm or bottomness. The yield 
of non-photonic e from open-heavy-flavor sources is 
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(9.6) 



where Bj/^ (B r ) is the yield of photonic e ± from the decays of J/xjj (T). Figure 9.2 
shows the corrected e yields Nec, Yhf, arid Ye>,b for the 0-20% most central Cu 
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+ Cu collisions, as well as the various components of the residual background. The 
non-photonic e yields Yhf and Yd,b are presented for three centrality classes in 
Section 19.31 



9.2 Uncertainties 

For each transverse-momentum bin, the statistical uncertainty of the uncorrected 

photonic e ± yield, P, is calculated based on the assumptions that n and n ++ are 

statistically independent, and that the statistical uncertainty of the combinatorial 

background n ++ is uncorrelated with the uncertainty in n un uk e - For n un uk e , 

n , and n ++ , the uncertainty of each quantity is assumed to be the square root of 

that quantity. The statistical uncertainty of P is 



oP = y/n un i ike + n__ + n ++ . 
The statistical uncertainty of Iec is 



<tI. 



al 



EC 



■ I 



EC 



Iec 

7f 



(9.7) 



(9. 



and the statistical uncertainty of Pec is 



aP, 



oP 



EC 



■P, 



EC- 



(9.9) 



As shown in Section B.2, the statistical uncertainty of Nec is calculated based on 
the assumption that the inclusive e yield is statistically independent of the ratio 
r = P/Ki nc JeB, the fraction of e that are photonic. The statistical uncertainty of 
N E c is found to be 



EC 



N 



EC \1 JT 



oP 

~p 



(9.10) 



The absolute statistical uncertainties of Yhf and Yd,b are also assumed to be a Nec] 
the uncertainties in the components of the residual background are accounted for in 
the systematic uncertainties. 

The systematic uncertainties in I EC , Pec, and N E c are calculated by assuming 
that the uncertainties in each of the correction factors are uncorrelated. See Sec- 



tion B.3 for the equations used. The systematic uncertainties of Yhf and Yd b are 
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the quadrature sums of the systematic uncertainty of N EC and the systematic uncer- 
tainties of the components of the residual background that were subtracted. 



9.3 Bin-Shift Correction 

In the next section, the non-photonic e yield (Xd,b) in Cu + Cu collisions will be 
compared to the non-photonic e yield in p + p collisions measured by the STAR |80] 



and PHENIX 51 collaborations, and the nuclear modification factor will be calcu- 



lated. In Reference [80] one more correction is applied to Yd,b'- the bin-shift correc- 
tion factor Cbs- The calculation of Cb s for non-photonic e in Cu + Cu collisions is 
described in this section. 

In Reference (80] the non-photonic e ± cross-section in p + p collisions is presented 
as the value of Ed 3 a jdp 3 at the center of each bin. In contrast, the non-photonic e 
yield in Cu + Cu collisions has been presented in this dissertation as the average value 
of the invariant yield Ed 3 N/dp 3 for each bin (i.e., the integral of Ed 3 N/dp 3 over the 
bin, divided by the bin width). Unless Ed 3 N/dp 3 is constant or varies linearly across 
the bin, the average value will not be the same as the value of Ed 3 N/dp 3 at the bin 
center. The bin-shift correction factor Cb s accounts for this difference. 

For each Cu + Cu centrality class, the non-photonic e ± yield Yd,b is fit with a 
function of the form / = A(l + pT/B) n . In the calculation of x 2 ; the fit intergal in 
each bin is compared to the measured value of the spectrum. For each transverse- 
momentum bin, Cbs is the ratio of the value of / at the bin center to the average 



value of /. This method is illustrated in Figure 9.3 for the 0-20% centrality class. 



In the upper panel, the gray points show the values of Yd,b- The green curve is a 



fit to that spectrum (the fit parameters are given in Section D.6) and the black lines 
represent the fit integral in each bin. For each bin, the bin-shift correction is the 
ratio of the green curve evaluated at the bin center to the fit integral. The values 
of Cbs are shown in the lower panel. The dashed curves are used to calculate the 
systematic uncertainties of C& s . The "soft fit" ("hard fit") curve is generated from 
the fit by decreasing (increasing) the parameter n by its uncertainty]^] The bin-shift 
correction is calculated using these two extremal curves; the uncertainties in Cb s are 
the difference between the extreme lower or upper calculations and the central value. 



3 The dashed curves are also multiplied by scale factors, but this has no effect on the calculation 
of the bin-shift correction. 
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For each centrality class, the non-photonic e ± yield Yd,b is multiplied by Cb s , giving 
the final corrected spectrum]^] These spectra are shown for three centrality classes in 



Figures 9A\ (for Y H f) and WJ> (for Y DjB )- 



Calculation of Bin-Shift Correction 
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Figure 9.3: Upper panel: illustration of the calculation of the bin-shift correction Cb s 
for the 0-20% centrality class of Cu + Cu collisions. See the text for a description of 
the quantities plotted. Lower panel: the bin-shift correction. 



4 An alternative approach, described in 133 , is to shift the abscissa of each data point away from 
the bin center to the point where the fit function is equal to the fit integral (leaving the ordinate 
unchanged). 
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Non-Photonic e ± Spectra, 200-GeV Cu+Cu Collisions 
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Figure 9.4: The yields (Yhf) of e 1 * 1 from heavy-flavor decays as functions of transverse 
momentum for three centrality classes of 200-GeV Cu + Cu collisions. 
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Figure 9.5: The yields (Yd,b) of ^ from open-heavy-flavor decays (D and B meson 
decays) as functions of transverse momentum for three centrality classes of 200-GeV 
Cu + Cu collisions. 
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9.4 Nuclear Modification Factor 



The nuclear modification factor, which was defined in Section 1.4.3 (page [l7j) , is used 
to quantify the differences in particle yields between nucleus-nucleus collisions and 
proton-proton collisions. This factor is the ratio of the yield of a particle species in 
nucleus-nucleus collisions to the yield of the same particle species in p + p collisions. 
The ratio is scaled to account for the fact that a nucleus-nucleus collision consists of 
many nucleon-nucleon collisions. For a given p? bin, the nuclear modification factor 
for particle species X is 



[Ed 3 N x 



AA 



\ dp 3 



r x aa = — • (9-ii; 



X 

1 (N bi n) 



dp 3 

Here, Ed 3 N AA /dp 3 is the invariant yield of particle species X in nucleus-nucleus 
collisions and Ed 3 a v x p j dp 3 is the cross-section for particle species X in proton-proton 
collisions. The inelastic nucleon-nucleon cross-section is <Jj ne ;. = 42 mb, and (N^n) is 
the mean number of binary nucleon-nucleon collisions for the A + A collision system. 
A value of R* A = 1 indicates that the yield of particle species X is not suppressed in 
nucleus-nucleus collisions relative to proton-proton collisions. The symbol u Raa is 
generic; the subscripts may be used to indicate a specific collision system, which will 
be the practice in this chapter. 



Equation 9.11 is used to calculate RquOu^ the nuclear modification factor for non- 
photonic e ± in 200-GeV Cu + Cu collisions (in comparison to 200-GeV p+p collisions). 



The STAR 80 and PHENIX [51] collaborations have measured the cross-sections of 
non-photonic e in 200-GeV p + p collisions. In both of those measurements, the 
contributions from J/ip and T decays have been removed; those spectra are therefore 
measurements in p + p collisions of the quantity called Yd # in this dissertation. 



The upper panel of Figure 9.6 shows the STAR measurement of Yd.d in V + P an d 



the 0-20% most-central Cu + Cu collisions; the Cu + Cu data have been scaled by 



1/(Taa) — Vinei.J (Nbi n )- The lower panel of Figure 9.6 shows the nuclear modification 
factor i?cuCu as a function of transverse momentum. The p + p spectrum is fit with 
a function of the form A(l + pj</B) n . The parameters of this fit function, called 



F[ejv,2], are given in Section D.6. The lower (upper) systematic uncertainty in the 



fit is estimated by moving all p + p data points down (up) by the quadrature sum 
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of their statistical and systematic uncertainties, then fitting the resulting spectrum 
with another function of the same form. The statistical uncertainties of -RcuCu are 
due to the statistical uncertainties of the Cu + Cu spectrum only; the systematic 
uncertainties of the Cu + Cu spectrum and the p+p fit are assumed to be uncorrelated. 



Figure 9.7 is similar to Figure 9.6, but the PHENIX p + p measurement is used in 
place of the STAR p + p measurement. The nuclear modification factor i?cuCu * s 
calculated by dividing the scaled Cu + Cu spectrum by a fit (called F[eN, 1], of the 
same functional form as F[e^, 2]) to the PHENIX p + p measurement of Yd.b- 
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• STAR p+p 
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Figure 9.6: Illustration of the calculation of -R^uOu' the miclear modification factor 
for non-photonic e in the 0-20% centrality class of 200-GeV Cu + Cu collisions. 
The upper panel shows the non-photonic e ± yield {Yd,b) found in this analysis as 
a function of pt along with the STAR collaboration's measurement 80 of the same 
quantity in p + p collisions (and a fit to those data). The lower panel shows the 
nuclear modification factor, the scaled ratio of the Cu + Cu measurement to the p + p 
measurement, also ClS db function of pr- 
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Figure 9.7: Illustration of the calculation of -RquCu' the nuclear modification factor for 
non-photonic e ± in the 0-20% centrality class of 200-GeV Cu + Cu collisions (see also 
Figure 9.6). The upper panel shows the non-photonic e ± yield {Yd,b) found in this 



analysis as a function of px along with the PHENIX collaboration's measurement 51 
of the same quantity in p + p collisions (and a fit to those data). The lower panel 
shows the nuclear modification factor, the scaled ratio of the Cu + Cu measurement 
to the p + p fit, also as a function of px- 
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9.5 Discussion 
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Figure 9.8: Nuclear modification factors of non-photonic e as functions of in Cu 
+ Cu and Au + Au collisions for centrality classes with similar values of (N part ). In 
the upper panel, -RcuCu f° r the 0-20% centrality class is compared to the PHENIX 
collaboration's measurements 52 of R^aL f° r the 20-40% and 40-60% centrality 



classes. In the lower panel, -RcuCu f° r the 20-60% centrality class is compared to the 
PHENIX collaboration's measurements of R^aL f° r the 40-60% and 60-92% centrality 
classes. The gray shaded boxes represent the systematic uncertainties of R^aL f° r 
the 40-60% centrality class, while the brown shaded boxes represent the systematic 
uncertainties of R^aL f° r the other two centrality classes. 



In this section, Rq^q u , the nuclear modification factor for non-photonic e in 
200-GeV Cu + Cu collisions, will be compared with related measurements of the 
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nuclear modification factor. Figure 9.8 shows the nuclear modification factors for non- 
photonic e in 200-GeV Cu + Cu and Au + Au collisions 52 with similar values of 
(Npart)- The upper panel shows -RcuCuo-20%' ^ e nuc lear modification factor for non- 
photonic e ± in the 0-20% most central Cu + Cu collisions (for which {N part ) = 86); the 
lower panel shows -R^uOu 20-60%' the nuc lear modification factor for non-photonic e ± 
in the 20-60% centrality class of Cu + Cu collisions ((N part ) = 33) j^] Measurements of 
the nuclear modification factor for non-photonic e in three centrality classes of Au + 
Au collisions are also shown: R^aL 20-40% (impart) = 140) is shown in the upper panel, 
^^,40-60% (( N part) = 60) is shown in both panels, and #^,60-92% (( N P art) = 14) 
is shown in the lower panel. In the upper panel, the values of R^f u 20-40% tend to 
be lower than the values of R^aI 4o-60% : non-photonic e ± appear to be suppressed 
more when (N part ) (and therefore the volume of the QGP) is greater. This behavior 
is typical and is observed in several collision systems for various particle species. 
However, the measured values of i?AuAu 60-92% deviate from this trend: the values of 
-^AuAu 60-92% a PP ear f° be below the measured values of ^AuAu4o-60%- F° r both Cu 
+ Cu centrality classes, the measured values of i?cuCu are generally consistent with 
the measured values of RjIuaL f° r collision systems with similar values of (N part ). 

Measurements of the nuclear modification factor in Au + Au collisions indicate 
that, at high p?, the suppression of non-photonic e ± |5lJ|52j|55j[70) is similar in mag- 
nitude to the suppression of pions 30 125] j^] Figure |9.9| shows measurements of 
i?CuCu for non-photonic e ± and charged pions |lll in two centrality classes]^] The 
upper panel shows -RcuCuo-20% plotted with -RcuCuo-20%' while the lower panel shows 
-^CuCu 20-60% plotted with Rq uCu 20 _60% • At high transverse momenta, the measured 
values of -RcuCuo-20% anc l -^CuCuo-20% are consistent with each other. However, for 
the 20-60% centrality class, non-photonic e ± appear to be more suppressed at high 
Pt than pions. 



Figure 9.10 shows measurements (as functions of {N part )) of the nuclear modifica- 
tion factor at high transverse momenta for non-photonic e ± and ^ in 200-GeV Cu 
+ Cu and Au + Au collisions. These Raa measurements are the scaled ratios of the 
integrated particle yields at high pr- The measured values of -RauAu are consistent 
with the values of -RauAu; except for the most peripheral centrality class (with the 



5 (Np ar t) = 51 for the 0-60% most centr al Cu + Cu collisions. 
6 See also the discussion of Figure 



9.10 



below. 



7 Note that -RJuCu increases as collisions become more periphera l. T his trend is expected, and 
differs from the behavior of R^uaL seen m the lower panel of Figure 9.8 
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Figure 9.9: Nuclear modification factors of non-photonic e and ^ 111 in Cu + Cu 
collisions as functions of pt- 



lowest value of (N part )), where non-photonic e ± may be more suppressed than pions. 
In Cu + Cu collisions, it appears that non-photonic may be more suppressed than 
pions. This is particularly visible for the peripheral 20-60% and combined 0-60% cen- 
trality classes ((N part ) = 33 and 51, respectively). For the 0-20% centrality class, the 
measured value of i?cuCu i s ^ ess than Rq uCu , though the measurements are consistent 
with uncertainties. 

A reason for measuring -RquCu was to determine if the shape of the QGP, not 
just its size, plays a role in the suppression of non-photonic e ± and heavy quarks. 
A Cu + Cu collision in the 0-20% centrality class and a Au + Au collision in the 
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Figure 9.10: High-pr nuclear modification factor of non-photonic e in Cu + Cu 
collisions (green) and Au + Au collisions (red) 152) as a function of (N part ) compared 
with the nuclear modification factors of 7r =l= in Cu + Cu collisions (blue) |111| and Au 
+ Au collisions (black) 11251. 



40-60% centrality class have similar values of (N part ) and should therefore produce 
similar volumes of quark-gluon plasma. A Cu + Cu collision in the 0-20% centrality 
class should have a roughly circular cross section (in the xy-plane), while a Au + Au 
collision in the 40-60% centrality class would have a more elliptical cross-section. On 
average, a heavy quark produced near the center of the medium would have a greater 
path length through a quark-gluon plasma produced in the Cu + Cu collision than 
in the Au + Au collision, which might be expected^ to result in more heavy-quark 
energy loss for Cu + Cu collisions, and more suppression of high-p^ non-photonic 



e . As shown in Figure |9.10[ the measurements of the nuclear modification factor 



for non-photonic in Cu + Cu collisions are consistent within uncertainties with 
measurements of the same quantity in Au + Au collisions at similar values of (N part ). 
It appears that the central value of -R^uOu 0-20% (Npart) = 86) may lie somewhat 
below the expected value based on the measurements in Au + Au (if R^uAu varied 
linearly with (N part ) between {N part ) = 60 and (N part ) = 140). However, the large 



3 assuming that the two collision systems have the same energy density 
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uncertainties on the R^a E measurements prevent a strong statement from being made 
regarding this issue. 
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Figure 9.11: The measured nuclear modification factor for non-photonic e in Cu 
+ Cu collisions compared with theoretical predictions. The range of Fff PL 



predicted by the Collisional Dissociation model [67 ,68 is shown as the dashed curves. 



^CuCu.O-lO^o 



Predictions of -RauAu by the Gossiaux/ Aichelin model 65 66 for three Au + Au 
centrality classes are shown as solid curves. 
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Figure 9. 1 1| shows the measured values of RquCu compared to predictions by two 
theoretical model^Jfor the nuclear modification factor of non-photonic e ± . The mea- 
sured values of Rq^cu 0-20% are consistent with the prediction of the Collisional Disso- 



ciation model 



67 



(38 



for -RcuCuo-10%- The measured values of -RcuCu f° r the 20-60% 
centrality class are generally consistent with the prediction of the Collisional Disso- 
ciation model for the 0-10% most central Cu + Cu collisions. The values of R^aL 



measured 55 by the PHENIX collaboration are well described for all but the most 
peripheral centrality class by the predictions of the Gossiaux/Aichelin model |65|66jp^ 



The measured values of -RcuCuo-20% ( u PPer panel) are closer to the Gossiaux/Aichelin 
prediction of R^f u 20-40% than R^f u 40-60%' even though the value of (N part ) for the 
20% most central Cu + Cu collisions is closer to the value of (N part ) for the 40-60% 
centrality class of Au + Au collisions. The measured values of Rq^q u 20-60% (l° wer 
panel) are consistent with the Gossiaux/Aichelin prediction for R^nAm 20 — 40% de- 
spite the factor of ~ 4 difference in (N part ) between the two centrality classes. 



9 These models are among those discussed in Section 1.5.2 

ln This model predicts that the values of R^uAu wm increase as collisions become more peripheral. 
This trend is supported by the PHENIX measurements, except for the measurement in the most 
peripheral (60-92%) centrality class. 
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In summary, this dissertation has presented the first measurement of the yields 
of non-photonic e in 200-GeV Cu + Cu collisions. The limited number of events 
analyzed, the large photonic e ± background, and limited knowledge of the high-pr 
J I ip spectrum in Cu + Cu collisions contribute to the large uncertainties in this mea- 
surement. Non-photonic e are suppressed at high transverse momentum in Cu + Cu 
collisions relative to p + p collisions. For Cu + Cu collisions, non-photonic e appear 
to be more suppressed at high px than charged pions, with the difference in suppres- 
sion becoming larger for more peripheral collisions. In contrast, measurements in Au 
+ Au collisions indicate that the suppression of non-photonic e ± is similar in magni- 
tude to the suppression of pions, although non-photonic e may be more suppressed 
than pions in collisions with low values of (N part ). Although the measured values of 
RquGu have large uncertainties, they are generally consistent with measurements of 
^AuAu f° r collision systems with similar values of (N part ). The suppression of non- 
photonic e measured in this dissertation is generally greater than the predictions of 
the Gossiaux/Aichelin model for similar values of (N part ). The Collisional Dissoci- 
ation model prediction is consistent with the observed suppression of non-photonic 
e ± for central Cu + Cu collisions. However, that model is also consistent with the 
measured suppression of non-photonic e for the 20-60% centrality class, despite a 
factor of ~ 3 difference in (N part ). 



9.6 Looking Forward 

The production of heavy quarks and their interactions with the quark-gluon plasma 
are active subjects of study at both RHIC and the Large Hadron Collider (LHC). It 
is intended that the non-photonic e ± spectra found in this dissertation will be com- 
bined with the D° yield measured in [42] to provide an updated measurement of the 
total charm cross-section in 200-GeV Cu + Cu collisions. The PHENIX collaboration 
has submitted for publication new measurements [55J of R^uaL an d the elliptic-flow 
parameter t>2 for non-photonic e . The STAR collaboration will soon submit a new 
measurement 80 of the non-photonic e yield in p + p collisions; a new measurement 
of the non-photonic e yield in Au + Au collisions is in progress. Efforts are also un- 
derway in the STAR collaboration to measure the flow of heavy quarks by measuring 
V2 for D° in Au + Au collisions. Measurements of the J / ip yield in p + p and d + Au 
collisions and of the T yield in d + Au collisions are also in progress. 
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The LHC has recently begun operations, opening a new energy regime of heavy-ion 



collisions for study. Pb + Pb collisions at ^snn = 2.76 TeV have been recorded 134 



along with p + p collisions at energies including 900 GeV, 2.36 TeV, and 7 TeV 135- 



137 . The total charm (bottom) cross-section is predicted [53] to increase by a factor 



of ~ 10 (~ 100) in collisions at A /s^ r /v — 5.5 TeV (the maximum energy for Pb + 
Pb collisions) in comparison to collisions at ^/snn = 200 GeV allowing heavy flavor 
to be measured at higher momenta and with increased statistical significance relative 
to RHIC. It may be possible to study the path- length dependence of heavy-flavor 
suppression by measuring heavy-flavor yields in the reaction plane and perpendicular 
to it (cf. the measurement of jet quenching in and out of the reaction plane shown 



in Figure 1.5 ). The ATLAS collaboration 138 has measured 139 the suppression of 



J/ip in Pb + Pb collisions. The ALICE experiment |140| is beginning a broad program 
of heavy-flavor measurements. The yields of D°(D°) and D^ 1 in heavy-ion collisions 
have been measured 141 by reconstructing their hadronic decays. Simulations |142| 
indicate that for ten million 5. 5- TeV Pb + Pb collisions in the 0-5% centrality class, 
the yield of D° can be measured at momenta as high as 20GeV/c with a relative 
statistical uncertainty of < 20%. The yields of heavy quarkonia may be measured by 
reconstructing their decays to e~e + pairs at mid-rapidity or pairs at forward 

rapidity (2.5 < y < 4). Measurements of single muons and non-photonic are 
also underway. In the ALICE experiment, e ± are identified through measurements 
of energy loss in the Time Projection Chamber, particle speed in the Time-of-Flight 
detector, the radiation yield in the Transition Radiation Detector, and/or energy 
deposition in the Electromagnetic Calorimeter. It is expected 142 that for ten million 
5.5-TeV Pb + Pb collisions in the 0-5% centrality class, the yield of non-photonic e 
from 5-meson decays may be measured at momenta as high as 20 GeV/ c with a 
relative statistical uncertainty of < 10%. 
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Appendix A 

Conversion of (dE/dx) to ncr e 



This appendix describes the method used to estimate the most probable value of the 
normalized TPC energy loss (na e ) for a given value of (dE/dx), the 70% truncated 
mean TPC energy loss. See Section 2.3.2 and Chapter [6] for more information on 
(dE/dx) and na e . The method described in this appendix is used to generate the 
curves shown in several figures, including Figure |3~To| (page 65). 

Figure A.l shows the relationship between \n(dE/dx) and na e for particles that 
pass all other e identification cuts. The relationship is not one-to-one because B e 
and a e depend upon the length of the track in the TPC. When fitting the inclusive e 
distributions with multiple Gaussians, it can be useful to know the expected locations 



of the means of the hadron peaks. Bichsel's parameterization 108 can be used to 
generate an expected value of (dE/dx) for hadron tracks with typical characteris- 
tics; this section presents a method to convert those values to an approximate most 
probable value of na e . 

The black data points shown in Figure A.l give an approximate most probable 
value of na e for each bin in Yn.(dE/dx). This value is found as follows. For each bin 
in ln(dE/dx), the mean value of na e and standard deviation are calculated. Then, 
the mean value of na e is recalculated including only those values of na e within one 
standard deviation of the original mean. This truncated mean (not related to the 70% 



truncated mean used to calculate (dE/dx)) is shown in Figure A.l The truncated 
mean seems to be better than the mean at reproducing the maximum of the n<j e 



distribution in each \n(dE/dx) bin. Figure A. 2 shows the mean na e , the standard 
deviation, and the truncated mean na e as a function of \n(dE/dx). 
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The truncated mean na e as a function of \n{dE/dx) can be fit with a linear 
function of the form 



na e = [ln(dE / dx) — lnf? 



eOj 



O'eO 



(a.i; 



The values of the fit parameters lni? e0 and l/cr e o are shown in Figure A. 3 The 
intercept In B e0 does not appear to depend on p for the momentum range of interest for 
this analysis; the slope l/a e o appears to increase slightly with increasing momentum. 
These fits are used to calculate an approximate most probable value of na e for a 
given input value of \n(dE/dx) obtained from Bichsel's parameterization. The curves 
shown in Figures |6.2 6.8, and 6.9 are generated using a single linear conversion 
function calculated for one large momentum bin. 
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200-GeV Cu+Cu, MB+HT, na 8 vs. ln<dE/dx>, cent. 0-60%, 2.0<p<11.5 
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Figure A.l: The normalized TPC energy loss na e vs. \a(dE/dx) for both minimum- 
bias and high-tower triggered events. The black circles indicate an estimate of the 
most probable value of na e (see the text for details). 

200-GeV Cu+Cu, MB+HT, no, vs. l^dE/dx), cent. 0-60%, 2.0<p<11.5 
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Figure A. 2: The mean (red, central values) value of na e as a function of ln(dE/dx), 
as well as the standard deviation (red boxes). Also shown (in black) is the estimated 
most probable (truncated mean) value of na e . 
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ln(dE/dx)-to-nc e Conversion Fit Parameters 
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Figure A. 3: The values of the intercept (lnS e0 ) and slope (l/er e0 ) of linear fits to the 
most probable value of na e vs. hx(dE/dx). These values are used to convert a value 
of (dE/dx) given by Bichsel's parametrization |108| to an approximate most probable 
value of na e . 
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Appendix B 



Treatment of Uncertainties 



This Appendix describes the calculation of uncertainties for three cases. Section |B.l 



describes the estimation of the statistical uncertainties for weighted efficiency calcu- 



lations. In Section B.2 the expression for the statistical uncertainty of the efficiency- 
corrected non-photonic e ± yield (Equation 9.10) is derived. Section B.3 describes the 



calculation of the systematic uncertainties in the efficiency-corrected e yields. 

B.l Statistical Uncertainties of 
Weighted Efficiencies 

In Chapters |4] and [7] the e ± reconstruction efficiency and the background rejection 
efficiency (eb) are calculated by finding the ratio (as a function of px) of two simulated 
e spectra. Those efficiencies are properly calculated as the ratio of two weighted 
spectra, although the unweighted ratios are sometimes shown for illustration purposes. 



The unweighted efficiencies are calculated using the "Bayes error formula" from 143 
which gives the variance of an efficiency when the numerator k is a subset of the 
denominator n and the two are highly correlated. The uncertainty in the efficiency 
e = k/n is 



k + l fk + 2 k + l\ 
aS Y n + 2\n + 3 n + 2J { ' ) 

For the weighted efficiencies, the statistical uncertainties in the efficiency e are 
estimated as follows. The uncertainty in e is calculated with the assumption that 
the numerator and the denominator both independently follow Poisson statistics; 
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Figure B.l: Ratios of Bayesian to Poissonian uncertainties as a function of the effi- 
ciency 6b = k/n for various values of n. This ratio (with n = 100) is used to convert 
a calculated Poissonian uncertainty in the efficiency calculation into an estimated 
Bayesian uncertainty. 



the resulting uncertainty is then multiplied by a conversion factor to generate an 
estimated error based on the Bayes error formula. This calculation is explained in 
greater detail below. It should be noted that the systematic uncertainties assigned 
to er and Eb are large enough to account for the statistical fluctuations in those 
efficiency calculations. 

For the sake of clarity, the calculation of the uncertainty of the background rejec- 
tion efficiency (eb) will be described. Only a single e ± transverse-momentum bin is 
considered and each e is weighted according to the pt of its parent photon. The nu- 
merator and denominator of the efficiency calculation are called k and n, respectively. 
The number of rejected e ± that come from a photon with transverse momentum pip 
is N re j.(p?p) , and the photon weighting function is W(p^). The numerator k is a 
weighted sum: 

k = J2lW(p y T )xN rej .(pl)} (B.2) 

P T 

If each measurement of N re j, follows Poisson statistics, the uncertainty of k would be 
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4 



-1 2 



W{jpl)^N rej .{jP T 

Pt ' Pt 

Similarly, if N Tec _ is the number of reconstructed e ± 



E [W\pD x N rej .{jP T )] 



(B.3) 



and 



n 



J2[W(pl) x N rec .(pl)} 



7 



[BA) 



°l = E k(PT)^W^)] 2 = E P*^) x iv re ,(p: 



7 



7 



(B.5) 



If k and n were statistically independent, the uncertainty of the ratio k/n would be 



2 /a \ 2 



W(t) +( 



(B.6) 



Of course, and n are not independent. In the unweighted case, when k and n are 



integers, the ratio of the uncertainty calculated using Equation B.l to the uncertainty 
calculated using the Poisson formula is 



R 



k+l fk+2 k+l 



ri 



n + 2\n + 3 n + 2J ky/l/k + l/n 



(B.7) 



n + 2 + 3 n + 2 J e By /l/(e B n) + ijn 



Figure B.l shows this ratio as a function of the efficiency eb = k/n for different 



values of n. For the range of efficiencies of interest for this analysis (0.4 < Eb < 0.8) 
this ratio exhibits less than a 10% change when n varies from 10 to 10 6 . For the 
weighted efficiency calculations shown in Chapter [7j equation B.7 (with n = 100) 
is used to convert the calculated Poissonian uncertainty to an estimated Bayesian 
uncertainty. The same arguments apply for the calculation of the uncertainties of the 
e ± reconstruction efficiency (sr)- The numerator is the number of reconstructed e , 
each weighted according its simulated transverse momentum; the denominator is the 
number of simulated e ± . 
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B.2 Statistical Uncertainties of 
the Non-Photonic e Yield 



In this section, Equation |9.10| for the statistical uncertainty of Nec, the efficiency- 
corrected non-photonic yield is derived, following 80 . The ratio (r) of the 



efficiency-corrected photonic e yield to the inclusive e yield is 

where P and I are the uncorrected photonic and inclusive e yields, respectively. The 
uncertainties in Ki nc , and Eb are accounted for in the systematic uncertainties of the 
spectra. It is assumed that the statistical uncertainties of I and r are uncorrelated. 
Therefore, 



aP 
~P 



al 
T 



'cry 
r / 



ar = r\ 



oP 
~P 



al 

T 



(B.9) 



where aX denotes the statistical uncertainty of quantity X. The efficiency-corrected 
non-photonic yield is 



N 



1 



-P 



(B.10) 



£b J Q 

where Q is shorthand for a combination of several correction factors, which do not 
influence the statistical uncertainty of Nec- The statistical uncertainty of Nec is 



aN 



EC 



N, 



EC 



arl 
T 



o-n — rl 



(B.11) 



Since a[l — r] = ar, 



aN, 



ec — N E c 



N, 



EC 




(B.12) 
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B.3 Systematic Uncertainties of e ± Spectra 



The systematic uncertainties of the inclusive, photonic, and non-photonic e spectra 
are due to the uncertainties in the various correction factors. In the calculation of the 
systematic uncertainties of the efficiency-corrected e ± yields {Iec, Pec, and Nec) 
it is assumed that the uncertainties in (A B emc), £ r, £ t, £ dE/dx, Nine, and sb are 
uncorrelated. In this analysis, the lower and upper systematic uncertainties are not 
required to be the same. In the following discussion, X[X] and v[X\ denote the lower 
and upper systematic uncertainties of quantity X, respectively. It is convenient to 
define the variable Q as the product of the four correction factors that are common 
to the calculations of Iec, Pec, and Nec- 

Q = (AbEMc) £ R " £ T ■ £ dE/dx- (B.13) 

Therefore, 



^-'/HWBW (R14) 

and 



v[Q] = Qi 



v[{A B emc)]\ A , ( v[ £ r]\ 2 , ( v[ £ t\\ , ( v[e dE /dx 



+ 



(Abemc) J V £ R 
The systematic uncertainties of Iec are 



+ I + 

£ T J \ £ dE/dx 



[B.15) 



x[iEc\=iEc^m 2 J m *> 



Q 



K 

- 1 - 7,r. 



and 



XOY (cK, 1 
v[Iec] = Iec\I I I + 



Q 



K, 



The systematic uncertainties of Pec are 



(B.16) 



(B.17) 



A [P EC ] = Pec\ 



V \Q]\ + ( v \ £ b\ 



Q 



£ B 



[B.18) 



and 
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^■w(f) ,+ (^)" 



The systematic uncertainties of N E c are 



k n ec\ = \ii ^^) 2 + (^mru,]f+ (^ A ^]) 2 ( R2 °) 



and 
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Appendix C 



Additional Material 



This appendix contains supplementary figures which were excluded from the main 
body of the dissertation for the sake of brevity. Additional figures related to the 
calculation of the reconstruction efficiency (er, see Chapter El) are presented in 



Section C.I Additional figures and tables related to the inclusive e purity (K^ 



see Section 6.3) are presented in Section C.2 



C.l e Reconstruction Efficiency 

This section contains additional figures relating to the calculation of the e recon- 
struction efficiency (er), which is described in detail in Chapter |4| 

C.l.l Calculated Efficiencies from Embedding 

This section contains additional plots of Er, calculated from embedding (simulation 
data set Ni), with the values of various e ± identification cuts changed from their 



standard values. These plots are an extension of the series presented in Section 4.3 
In these figures, the black circles indicate the calculation of Er with all cuts having 
their standard values. 



Figures |C.1| shows Er calculated with various values for the cut on the BSMD 



cluster displacement A$md = \/(^Vsmd) 2 + {^■4>smd) 2 - Changing the value of the 
A-smd cut above 0.01 has negligible effect on er. 



Figure |C.2| shows the e reconstruction efficiency calculated with various values 



for the cut on Ritpc- 
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Figure C.3 shows Er calculated with various values for the global DCA cut. No 
strong dependence of the efficiency on the GDC A cut is observed. 



Figures C.4 and C.5 show er calculated with various values for the cuts on Ntpcju 
and Rfp, respectively. 



e± Rec. Eff. (e R ), MB trigger, cent. 0-60% 
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Figure C.I: The e reconstruction efficiency as a function of px for various values of 
the cut on A$md- 
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e 1 Rec. Eff. (e R ), MB trigger, cent. 0-60% 
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Figure C.2: The reconstruction efficiency as a function of px for various values of 
the cut on R\t pe- 
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Figure C.3: The e reconstruction efficiency as a function of px for various values of 
the GDC A cut. No large change in Er is observed as the value of the GDC A cut is 
varied above 0.5 cm. 
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e 1 Rec. Eff. (e R ), MB trigger, cent. 0-60% 
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Figure C.4: The e ± reconstruction efficiency as a function of px for various values of 
the cut on N TPC fit- 
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>. 0.7 

u 

c 

0) 

£ 0.6 
«^ 

HI 

I 0.5 

u 



£ 0.4 

o 

o 

I 0-3 



0.2 
0.1 

0, 




• R FP > 0.52 

□ R FP >0.7 

i R FP > 0.8 

a R FP >0.9 



_L 



_L 



10 12 14 

p r [GeV/c] 



Figure C.5: The e reconstruction efficiency as a function of for various values of 
the cut on Rfp- 
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C.1.2 Comparison of Embedding data to Real Data 



The plots in this section are an extension of the series presented in Section 4.4 



Figure C.6 shows the distributions of Nxpcfit and Figure C.7 shows the partial 
efficiencies of the N TPC fit > 20 cut. 

Figure C.8 shows the distributions of Rfp'i the partial efficiencies of the Rpp > 
0.52 cut are all very close to unity and are not shown. 



Figure |C.9| shows the distributions of R\tpc and Figure |C.10| shows the partial 
efficiencies of the R\tpc < 102 cm cut. 
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MB Trigger, 1^. Distributions, cent. 0-60%, 2.0<p T <3.0 
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Figure C.6: Distributions of N TPC fit for real and simulated e ± . 
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Figure C.7: Partial efficiencies of the cut Nxpcfit > 20. For pt < 4GeV/c, the 
results are for minimum-bias events; for p? > 4 GeV/ c, the results are for high-tower 
triggered events. 
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MB Trigger, R pp Distributions, cent. 0-60%, 2.0<p T <3.0 
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Figure C.8: Distributions of the TPC-fit-points ratio (Rfp) for real and simulated 
e ± . Nearly all e ± pass the Rfp > 0.52 cut and the partial efficiency is very close to 
1. 
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MB Trigger, R Distributions, cent. 0-60%, 2.0<p T <3.0 
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Figure C.9: Distributions of R\tpc for rea l an d simulated e . 
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Figure CIO: Partial efficiencies of the cut R\tpc < 102 cm. 
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C.2 Inclusive e Purity 



This section contains additional figures and tables relating to the calculation of the 
inclusive e ± purity (K inc ), which is described in detail in Section 6.3 Figures C.ll 



through C.24 are plots of the the multi-Gaussian fits to the na e distributions for the 
0-60% centrality class. In each figure, the green peak accounts for the e contribution, 
while the blue and red peaks account for hadron contributions (with the blue peak 



intended to account for charged pions). Tables C.I and C.2 contain the calculated 
values of K inc , as a function of pt for two event trigger types and three centrality 
classes. 




Figure C.ll: Four-Gaussian fit of na e distribution: 0-60% centrality class of 
minimum-bias triggered events; 1.3GeV/c < p < 1.7GeV/c. The magenta peak 
accounts for the (anti)deuteron contribution, (identical to Figure 6.10, page 119) 
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Figure C.12: Three- Gaussian fit of na e distribution: 0-60% centrality class of 



minimum-bias triggered events; 2GeV/c < p < 2.5GeV/c. (identical to Figure 6.7 
page 



115) 



na e Distribution for 2.5<p<3.0, MB Trigger, cent. 0-60% 




Figure C.13: Three- Gaussian fit of na e distribution: 0-60% centrality class of 
minimum-bias triggered events; 2.5 GeV/c < p < 3 GeV/c. 
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Figure C.14: Three- Gaussian fit of na e distribution: 0-60% centrality class of 
minimum-bias triggered events; 3 GeV/c < p < 3.5 GeV/c. 



np e Distribution for 3.5<p<4.0, MB Trigger, cent. 0-60% | 




Figure C.15: Three- Gaussian fit of na e distribution: 0-60% centrality class of 
minimum-bias triggered events; 3.5 GeV/c < p < 4 GeV/c. 
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Figure C.16: Three- Gaussian fit of na e distribution: 0-60% centrality class of 
minimum-bias triggered events; 4GeV/c < p < 5 GeV/c. 
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Figure C.17: Three- Gaussian fit of na e distribution: 0-60% centrality class of 
minimum-bias triggered events; 5 GeV/c < p < 6.5 GeV/c. 
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Figure C.19: Three-Gaussian fit of na e distribution: 0-60% centrality class of high- 
tower triggered events; 4GeV/c < p < 5 GeV/c. 
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Figure C.20: Three-Gaussian fit of na e distribution: 0-60% centrality class of high- 
tower triggered events; 5 GeV/c < p < 6 GeV/c. 



no e Distribution for 6.0<p<7.0, HT Trigger, cent. 0-60% 




Figure C.21: Three-Gaussian fit of na e distribution: 0-60% centrality class of high- 
tower triggered events; 6GeV/c < p < 7GeV/c. 
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Figure C.22: Three-Gaussian fit of na e distribution: 0-60% centrality class of high- 
tower triggered events; 7GeV/c < p < 8 GeV/c. 



no e Distribution for 8.0<p<9.0, HT Trigger, cent. 0-60% 




Figure C.23: Three-Gaussian fit of na e distribution: 0-60% centrality class of high- 
tower triggered events; 8 GeV/c < p < 9 GeV/c. 



213 




Figure C.24: Three-Gaussian fit of na e distribution: 0-60% centrality class of high- 
tower triggered events; 9GeV/c < p < 11.3GeV/c. 
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Table C.l: Inclusive e ± purity (K inc ) as a function of pt for the minimum-bias 
triggered data set. 



FT 


^max 
FT 




lower 


upper 


[GeV/c] 


[GeV/c] 


purity 


uncertainty 


uncertainty 


MB trigger, centrality 0-20% 


2.0 


2.5 


0.97009 


0.01152 


0.01092 


2.5 


3.0 


0.93983 


0.01851 


0.01632 


3.0 


3.5 


0.88753 


0.04563 


0.04009 


3.5 


4.0 


0.84976 


0.10105 


0.08891 


MB trigger, centrality 20-60% 


2.0 


2.5 


0.97647 


0.00977 


0.00945 


2.5 


3.0 


0.96144 


0.01309 


0.01115 


3.0 


3.5 


0.93288 


0.03145 


0.02550 


3.5 


4.0 


0.87917 


0.06763 


0.05778 


4.0 


4.5 


0.83746 


0.10990 


0.08943 


MB trigger, centrality 0-60% 


2.0 


2.5 


0.97192 


0.01088 


0.01042 


2.5 


3.0 


0.95003 


0.01716 


0.01476 


3.0 


3.5 


0.90798 


0.04234 


0.03586 


3.5 


4.0 


0.87411 


0.08126 


0.06815 


4.0 


4.5 


0.82238 


0.13420 


0.11626 
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Table C.2: Inclusive e ± purity {K inc ) as a function of pt for the high-tower triggered 
data set. 



FT 


^.max 
FT 




lower 


upper 


[GeV/c] 


[GeV/c] 


purity 


uncertainty 


uncertainty 


HT trigger, centrality 0-20% 


3.5 


4.0 


0.98305 


0.01097 


0.00777 


4.0 


4.5 


0.97306 


0.00807 


0.00715 


4.5 


5.0 


0.95449 


0.01319 


0.01195 


5.0 


5.5 


0.92818 


0.02272 


0.02054 


5.5 


6.0 


0.91508 


0.02787 


0.02516 


6.0 


6.5 


0.88956 


0.03436 


0.03138 


6.5 


7.0 


0.84440 


0.03850 


0.03681 


7.0 


7.5 


0.80547 


0.04324 


0.04287 


7.5 


8.0 


0.75479 


0.06429 


0.06755 


8.0 


8.5 


0.73205 


0.06838 


0.07179 


8.5 


9.0 


0.66771 


0.10601 


0.11702 


HT trigger, centrality 20-60% 


3.5 


4.0 


0.98731 


0.00734 


0.00485 


4.0 


4.5 


0.97583 


0.00754 


0.00658 


4.5 


5.0 


0.96111 


0.01094 


0.00982 


5.0 


5.5 


0.92848 


0.01859 


0.01725 


5.5 


6.0 


0.90131 


0.02533 


0.02391 


6.0 


6.5 


0.86248 


0.03510 


0.03356 


6.5 


7.0 


0.82938 


0.04072 


0.03925 


7.0 


7.5 


0.78388 


0.04833 


0.04759 


7.5 


8.0 


0.67611 


0.06183 


0.06770 


HT trigger, centrality 0-60% 


3.5 


4.0 


0.98477 


0.00896 


0.00637 


4.0 


4.5 


0.97136 


0.00837 


0.00746 


4.5 


5.0 


0.95180 


0.01261 


0.01155 


5.0 


5.5 


0.92555 


0.02076 


0.01912 


5.5 


6.0 


0.91411 


0.02615 


0.02403 


6.0 


6.5 


0.89065 


0.03350 


0.03105 


6.5 


7.0 


0.84732 


0.03840 


0.03677 


7.0 


7.5 


0.80261 


0.04532 


0.04410 


7.5 


8.0 


0.71752 


0.07533 


0.07068 


8.0 


8.5 


0.67326 


0.07989 


0.07604 


8.5 


9.0 


0.56928 


0.13911 


0.12581 
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Appendix D 
Catalog of Functions 



This Appendix is a catalog of functions used to describe the spectra of pions, rj 
mesons, direct photons, (anti)protons, J/ip, and e . Most of the functions take the 
form of a power-law term and/or an exponential term. Unless otherwise noted, the 
parameters A and C have units of mb • GeV~ 2 c 3 and the parameters B, D, and T 
have units of GeV/c or (GeV/c) 2 . 



D.l Pions 

The following is a list of functions which are used to describe pion spectra. They are 
used as weighting functions in calculating the background rejection efficiency Eb (see 
Chapter [7|. The function F[7T°,2] is used to construct (via rrix scaling) weighting 
functions for other meson species in calculating (see Chapter [8]) the yields of residual 
photonic e ± and non-photonic e ± from light-flavor sources (Y LFP ) and (Y LFN ). These 



functions are plotted in Figure 7.1 (page 124) 

Formula A(l+p T /B) n 
Parameters 393, 1.212, -9.97 



Data PHENIX tt° in p + p (2003) [119 
iV,2] 

Formula A(l + p 2 T /B) n + Ce~ PT l D 

Parameters 6.79322, 1.4428, -4.21776, 274.899, 0.141229 
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Data PHENIX vr° in p + p (2007) (l20 

Formula A(l+p T /B) n 
Parameters 1690, 0.723, -8.61 
Data STAR vr° in p + p [m) 

• F[vr°,4] 

Formula A(l +p T /B) n 

Parameters 13131.6, 0.558664, -8.65379 

Data PHENIX vr° in Cu + Cu 0-10% [m 

Note The data were scaled by 1/{{T AA )R AA ) = 42mb/(182.7 x 0.515) before 
fitting. 

• F[*±l] 

Formula A(l + p 2 T /B) n + Ce- pT ' D 

Parameters 2.07378, 2.12146, -4.32962, 128.454, 0.188331 
Data STAR vr° in Cu + Cu 0-60% (111) 

Note The data were scaled by 1/{(T AA )R AA ) = 42mb/(80.41 x 0.7) before 
fitting. 



D.2 rj Mesons 

The following is a list of functions which are used to describe 77-meson spectra. They 
are used as weighting functions in calculating the background rejection efficiency Eb 
(see Chapter [8]). These functions are plotted in Figure 7.2 (page 126). 



F[ V ,1] 

Formula A(l +p T /B) n 

Parameters 94.2433, 1.23777, -9.63473 



Data PHENIX r\ mesons in p + p 32 
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F[7],2] 

Formula A(l + q 2 / B) n + Ce~^ D , q = ^Jp 2 T + m(r]) 2 - m(n ) 2 
Parameters 0.48 x 6.79322, 1.4428, -4.21776, 0.48 x 274.899, 0.141229 
Note This is not a fit, but rather an expected spectrum derived from F[ir°, 2] 



(see Section 7.2). 



Data derived from PHENIX 7r° in p + p [120 
• F[ V ,3] 

Formula A(l +p T /B) n 
Parameters 70, 1.33, -9.83 



Data STAR rj mesons in p + p 117 



D.3 Photons 

The following is a list of functions which are used to describe photon spectra. They 
are used as weighting functions in calculating the background rejection efficiency Eb 
(see Chapter [8]). In the function names, 7^ indicates photons from 77-meson decays, 7^ 
indicates the non-thermal component of a direct photon spectrum, and 77- indicates 
thermal direct photons. 

Formula A(l +p T /B) n 
Parameters 47.426, 1.04317, -9.69906 

Note This is a fit to a simulated photon spectrum, for which F[r), 1] was used 



as a weighting function (see Section 7.2). 



Data derived from PHENIX 77 mesons in p + p 32 

Formula A(l +p T /B) n 

Parameters 418.887, 0.675643, -9.07911 

Note This is a fit to a simulated photon spectrum, for which F[rj, 2] was used 



as a weighting function (see Section 7.2). 
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Data derived from PHENIX tt° in p + p (l20 

Formula A(l+p T /B) n 

Parameters 35.4076, 1.11683, -9.88181 

Note This is a fit to a simulated photon spectrum, for which F[q, 3] was used 



as a weighting function (see Section 7.2). 



Data derived from STAR r\ mesons in p + p 117 

Formula A(l +p T /B) n 
Parameters 163.226, 0.279184, -6.3762 



Data PHENIX direct photons in p + p 122, 123 



Note F[jd, 1] and F[jd, 2] are fits to the same direct-photon spectrum. F[jd, 1] 
was generated by the author of this dissertation. 

F[id,2] 

Formula A(l+p 2 T /B) n 

Parameters 0.00919245, 1.798, -3.28487 



Data PHENIX direct photons in p + p 122, 123 



Note F[jd, 1] and F[y D , 2] are fits to the same direct-photon spectrum. F[y D , 2] 
is the published fit generated by the PHENIX Collaboration. 

Formula A(l +p T /B) n 
Parameters 0.090482, 1.36984, -7.35996 
Data PHENIX direct photons in Au + Au 0-92% [33]p2 



Note The data were fit with (T AA ) x F[j D ,3], with (T AA ) = 6.14mb _1 . 
Formula A(l+p T /B) n 
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Parameters 0.248478, 1.2668, -7.60368 



Data PHENIX direct photons in Au + Au 0-20% p 122 



Note The data were fit with (T AA ) x F[y D ,A], with (T AA ) = 18.6 nib -1 . 

F[y D ,5] 

Formula A(l +p T /B) n 

Parameters 2.14531 ■ 10 14 , 0.00217448, -6.05396 



Data PHENIX direct photons in Au + Au 20-40% m,122 



Note The data were fit with (T AA ) x F[j D ,5), with (T AA ) = 7.07 nib -1 . 
F[itA] 

Formula Ce~ PT/T 
Parameters 3.24934, 0.233 

Data PHENIX direct photons in Au + Au 0-92% p|p2] 

Note The data were fit with (T AA ) x (F[y D ,2] + F[y T , 1]), with (T AA ) 
6.14mb _1 . 

F[itA 

Formula Ce~ PT/T 
Parameters 4.38509, 0.221 

Data PHENIX direct photons in Au + Au 0-20% [33|p2] 

Note The data were fit with (T AA ) x (F[y D ,2] + F[^ T , 1]), with (T AA ) 
18.6 mb" 1 . 

F[ 7T ,3] 

Formula Ce~ Pr/T 
Parameters 5.62872, 0.217 



Data PHENIX direct photons in Au + Au 20-40% m,122 



Note The data were fit with (T AA ) x (F[j D ,2] + F[jr, 1]), with (T AA ) 
7.07 mb" 1 . 
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D.4 Protons 



The following is a list of functions which are used to describe proton spectra. The 
function F[p, 4] is used as a weighting function in calculating the yield of non-photonic 



e from baryon decays (which is found to be negligible, see Section 8.2). These 
functions are plotted in Figure 8.4 (page |i"57 ). 



• F[p,l] 

Formula A(l + p 2 T /B) n + Ce~ PT ' D 

Parameters 0.953737, 3.27513, -4.43296, 44.4136, 0.343277 



Data STAR (anti)protons in Au + Au 0-12% 125 



Note Parameters A and C are measured in GeV 2 c 3 . 
F\p,2] 

Formula A(l + p 2 T /B) n + Ce- pT ' D 

Parameters 0.703249, 3.34638, -4.42059, 32.2110, 0.344818 



Data STAR (anti)protons in Au + Au 10-20% 125 



Note Parameters A and C are measured in GeV 2 c 3 . 
F[p,3] 

Formula A(l + p 2 T /B) n + Ce~ PT ' D 

Parameters 0.122992, 3.62065, -4.07714, 18.3999, 0.347876 



Data STAR (anti)protons in Au + Au 20-40% 125 



Note Parameters A and C are measured in GeV 2 c 3 . 
Ftp, 4] 

Formula A(l + p 2 T /B) n + Ce- pT ' D 

Parameters 0.307271, 3.68843, -4.74017, 9.59262, 0.310598 



Data STAR (anti)protons in Au + Au 40-60% 125 



Note Parameters A and C are measured in GeV 2 c 3 . 



F[p,5] 
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Formula A{1 + p 2 T /B) n + Ce~ PT ^ D 

Parameters 1.26394, 1.49505, -4.29711, 0.205822, 0.338809 



Data STAR (anti)protons in Au + Au 60-80% 125 



Note Parameters A and C are measured in GeV 2 c 3 . 



D.5 J/tp 

These J/ip weighting functions are used calculating the yield of photonic e from 



J/ip decays (see Section 8.3.1). The fit functions for Cu + Cu collisions are plotted 
in Figure |8.8[ 

• F[J/ip,pp], F[J/ip,pp-], F[J/ip,pp+] 
Formula A[exp(—ap T — bp^) + Pt/PgY ti 

Parameters (5.23757, 0.323357, 0.0573737, 2.59389, 8.44212), (3.74232, 0.354190, 
0.0320993, 2.74964, 8.71965), (6.86912, 0.297896, 0.0749802, 2.50527, 8.27808) 



Data PHENIX (46] and STAR (47] J/ip in p + p 

Note Three fits were performed to obtain a "central" fit (F[J/ip,pp]), and 
lower (F[J/i/;,pp—]) and upper (F[J/ijj,pp+}) limits. Units of parameters: 
A in nb GeV~ 2 c 3 ; a in (GeV/c)^ 1 ; b in (GeV/c)~ 2 ; po in GeV/c; n unit-less 

F[J/if>, - 20%], F[J/ip, - 20%-], F[J/V», - 20%+] 

Formula A[exp(— ap? — bp^) + Pt/Po} " - 

Parameters (9.10922, 0.221385, -0.0100289 ,4.06224, 10.0477), (5.35754, 0.449842, 
-0.0632743, 2.87785, 5.89605), (11.5138, 0.278625, 0.00748965, 3.07351, 
8.44259) 



Data PHENIX (128] and STAR |47| J/ip in Cu + Cu 0-20% 

Note Three fits were performed to obtain a "central" fit (F[J/i[),Q — 20%]), 
and lower (F[J/i/>, 0-20%-]) and upper [F[J/if>, 0- 20%+]) limits. Units 
of parameters: A in nbGeV _2 c 3 ; a in (GeV/c) -1 ; b in (GeV/c)~ 2 ; p in 
GeV/c; n unit-less 

F[J/4>, 20 - 60%] and F[J/^, 20 - 60%-] 
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Formula A[exp(— ap? — bpj) +Pt/Po] n 

Parameters (3.71419, 0.113357, -0.0073416, 7.67173, 18.2388), (3.236, 0.111221, 
-0.0090455, 7.87763, 17.3856) 



Data PHENIX 1281 J /if) in Cu + Cu 20-60% 



Note Two fits were performed to obtain a "central" fit (F[J/ij), 20 — 60%]) and 



a lower limit (F[J/ip, 20 — 60%—]). See Section 8.3.1 for a description of the 
calculation of the upper limit (F[J/ip, 20 — 60%+]). Units of parameters: 
A in nbGeV" 2 c 3 ; a in (GeV/c) _1 ; b in (GeV/c)~ 2 ; p in GeV/c; n unit-less 

D.6 Non-Photonic e 

The following is a list of functions which are used to describe non-photonic e spectra. 

• F[e N , 1], F[e N , 1-], and F[e N , 1+] 
Formula A(l+p T /B) n 

Parameters (A, B, n) = (0. 0757745, 1.31713, -9.20756), (0.0691622, 1.27749, - 
9.05677), (0.167812, 1.02855, -8.54378) 

Data PHENIX NPE in p + p |5l] 

Note Three fits were performed to obtain a "central" fit {F[en, 1]), and lower 
(F[ejy, 1—]) and upper (F[e^f, 1+]) limits. 

• F[e N ,2], F[e N ,2-}, and F[e N ,2+] 
Formula A(l+p T /B) n 

Parameters (A, B, n) = (0. 00228122, 3.27957,-11.8724), (1.85453-10- 4 , 15.6626, 
-31.0754), (0.0499159, 1.23543, -8.62345) 



Data STAR NPE in p + p 80 



Note Three fits were performed to obtain a "central" fit (F[sn, 2]), and lower 
(F[epf, 2— ]) and upper (F[ejv, 2+]) limits. 

Formula A(l+p T / J B) n 

Parameters 3488.44, 0.208018, -7.80847 
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Data STAR NPE in Cu + Cu [this dissertation] 

Note fit to iteration of the efficiency-corrected non-photonic e ± yield (N EC ); 
parameter A in GeV _2 c 3 

• F[e N ,h] 

Formula A(l+p T /B) n 

Parameters 3866.20, 0.207808, -7.83930 

Data STAR NPE in Cu + Cu [this dissertation] 

Note fit to iteration 1 of the efficiency-corrected non-photonic e ± yield (Nec)', 
parameter A in GeV~ 2 c 3 
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